Next Article in Journal
Bipolar Fuzzy Relations
Previous Article in Journal
Moisture Estimation in Cabinet Dryers with Thin-Layer Relationships Using a Genetic Algorithm and Neural Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Nanolayer and Second Order Slip on Unsteady Nanofluid Flow Past a Wedge

Department of applied Mathematics, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(11), 1043; https://doi.org/10.3390/math7111043
Submission received: 29 August 2019 / Revised: 22 October 2019 / Accepted: 28 October 2019 / Published: 3 November 2019
(This article belongs to the Special Issue Nonlinear Waves and Applications )

Abstract

:
This paper presents the study of unsteady nanofluids flow and heat transfer past a wedge with second order velocity slip and temperature jump. The model is modified by considering the existence of a nanolayer together with the effects of thermophoresis and Brownian motion. The fundamental equations were transformed into ordinary differential equations by a new set of similarity transformations and solved by using the homotopy analysis method (HAM). We determined that the error reached 10−6 and the effectiveness of HAM was attained. The influence of second-order slip on the fluid skin-friction coefficient was analyzed and we determined that the Nusselt number decreases and skin friction coefficient rises with an increase in the thickness of the nanolayer.

1. Introduction

Nanofluids have been used in electronic information induction, scientific experiments recently. Nanofluids, which consist of nanometer-sized solid particle or tubes suspended in base fluids, are solid–liquid composite materials [1]. Nanofluids can flow smoothly with such ultrafine nanoparticles without clogging. Researchers [2] have proven that nanofluids can enhance the thermal transport property of the original fluid. On the other hand, some people have also noticed the issue that layered solid-like structures will be formed between the liquid and solid molecules. Nanolayers form a thermal bridge between a bulk fluid and a solid nanoparticle [3]. Some experiments [4,5] have proven that the nanolayer has a major effect on the thermal conductivity of nanofluid, while the particle diameter is no more than 10 nm. The thickness of the nanolayer can improve the heat transfer capacity of the nanofluid. However, we know little about the relationship between the thermal properties of the fluid and the nanolayer. It is very necessary to take the nanolayer into account and it indicates the direction for the development of new coolant [6].
In many practical situations, the velocity of free flow varies with time, which leads to instability of flow. For many engineering problems, such as periodic fluid motion, the unsteady boundary layer plays a vital role [7]. Unsteady magnetohydrodynamic (MHD) Falkner–Skan flow of Casson nanofluid was studied by Imran Ullah et al. [8]. Hydro-magnetic unsteady channel flow of nanofluid with heat transfer was studied by Awan et al. [9]. Unsteady flow and heat transfer past a shrink sheet in nanofluids was studied by Jahan et al. [10]. Some scholars also studied non-Newtonian nanofluid [11,12,13].
With further research on nanofluids, scholars have shown greater attention to the slip boundaries than the nonslip boundary conditions. Babu et al. [14] studied three-dimensional MHD slip flow of nanofluids over a stretching sheet with thermophoresis and Brownian motion effects; Usman et al. [15] studied thermal and velocity slip effects on Casson nanofluid flow; Abbas et al. [16] studied slip effects on stagnation point flow of a micro polar nanofluid; Muhammad Ramzan et al. [17] studied the partial slip effect of MHD micropolar nanofluid flow on rotating disks; Seth et al. [18] studied the transient flow of MHD nanofluid, considering Navier’s slip boundary condition. However, the literature only discusses the first order slip model. The degree of slip at the boundary depends on a number of interfacial parameters [19]. The Navier’s slip condition breaks down at higher shear rates. The values calculated by using second order slip boundary conditions are closer to the experimental value [20].
The unsteady fluid slip-flow issue has gained extensive interest from researchers, and some mathematical models have been solved by numerical methods. The main numerical methods used are the Runge–Kutta method, the shooting method, and the finite difference method [21]. Sobamowo et al. [22] studied the flow of nanofluid under the effect of slip conditions by numerical methods. Job et al. [23] studied unsteady MHD free convection nanofluid flows by the mixed finite element method. Because most of the fundamental governing equations have a strongly nonlinear and non-conventional nature, obtaining good approximate solutions through traditional methods is difficult. In order to solve this problem, the homotopy analysis method [24] was found to be a worthy and applicable method.
In the paper, we implement the homotopy analysis method to investigate unsteady nanofluids flow and heat transfer based on high order velocity slip and temperature jump. The equations are transformed by a new different similarity transformation and analyzed by HAM. The effects of nanolayer, velocity slip parameter, temperature jump on the Nusselt number and skin-friction coefficient are discussed. The effects of Biot number, solid volume fraction, thickness of nanolayer, and wedge angle on velocity and temperature are found as well.

2. Mathematical Formulation of the Problem

When adding nanoparticles to the base liquid, a layer of thin film will be covered around the nanoparticles. The nanolayer can improve the thermal conductivity of particles in the liquid suspension. Assuming that each particle layer can be combined with nanoparticles, the result is that the nanolayer can increase the volume concentration ϕ e , which can be calculated as [25]:
ϕ e = ( 1 + δ r ) 3 ϕ ,
where r is the radius of the nanoparticles and δ is the thickness of the nanolayer. Nanofluids are not a simple mixture of liquid and solid particles. The physical properties characterizing the nanofluids are strongly affected by the characteristics of the suspended particle and base fluids, such as particle morphology, chemical structure of particles and base fluids. Furthermore, the effective density of nanofluid ρ n f , effective heat capacity ( ρ C p ) n f , thermal conductivity κ n f , thermal diffusivity α n f , effective dynamic viscosity μ n f , electrical conductivity σ n f , respectively, depend on the nanoparticle volume fraction ϕ , δ , r as follows [26]:
ρ n f = ( 1 ϕ e ) ρ b f + ϕ e ρ s p ,   κ b f κ n f = κ s p + 2 κ b f + ϕ e ( κ b f κ s p ) κ s p + 2 κ b f 2 ϕ e ( κ b f κ s p ) ,
( ρ C p ) n f = ( 1 ϕ e ) ( ρ C p ) b f + ϕ e ( ρ C p ) s p , α n f = κ n f ( ρ C p ) n f , μ n f = μ b f ( 1 ϕ e ) 2.5 ,
σ n f = σ b f [ 1 + 3 ϕ e ( σ s p σ b f 1 ) ( σ s p σ b f + 2 ) ϕ e ( σ s p σ b f 1 ) ] ,
where subscripts s p , b f , n f correspond to solid particle, base fluid, and nanofluid, respectively.

2.1. Mathematical Modeling Analysis

Assuming that the flow is two-dimensional, laminar flow, and the research of nanofluid is a volatile and incompressible fluid, we also hypothesize that the influence of fluid flow in the process of evaporation and surface tension can be ignored. Considering the effect of the nanolayer and Brownian motion and high order slip, the corresponding mathematical model around the wedge is constructed as follows:
u x + v y = 0 ,
u t + u u x + v u y = U d U d x + μ n f ρ n f 2 u y 2 σ n f B 0 2 ρ n f [ u U ( x ) ] + [ ϕ e ( ρ β ) s p + ( 1 ϕ e ) ( ρ β ) b f ] g ρ n f ( T T ) × cos ( β π 2 ) ,
T t + u T x + v T y = α n f 2 T y 2 + τ D T T ( ( T x ) 2 + ( T y ) 2 ) + Q 0 ( ρ C p ) n f ( T T ) 1 ( ρ C p ) n f q r y ,
where u and v are the velocity components along the x and y directions, respectively, U ( x ) = a x λ t is the uniform velocity, a denotes the initial stretching rate and is a positive constant, λ is the stretching parameter, the angle of the wedge is denoted by βπ, and T indicates the temperature inside the boundary layer, B 0 is the magnetic field along the y-axis, T denotes the temperature far from the surface and is a constant, Q 0 is the heat generation coefficient, τ = ( ρ C p ) s p ( ρ C p ) b f is the ratio of effective heat capacity. D T is the thermophoresis diffusion coefficient, β is the thermal expansion coefficients, g is the gravitational force due to acceleration [27].
In this paper, Rosseland’s diffusion approximation was adopted for the radiative heat flux q r , and the expression was written as [28]:
q r = 4 σ * 3 k * T 4 y ,
where σ * is the Stefan–Boltzmann constant, k * denotes the Rosseland mean spectral absorption coefficient. Then, we expanded T4 in a Taylor series about T and neglected higher terms to gain:
T 4 4 T 3 T 3 T 4 .
Thus, we obtained:
q r = 16 σ * T 3 3 k * T y .
Based on the first order and second order velocity slip and temperature jump, the boundary conditions based on Kandasamy et al. [29] were modified as follows.
When t = 0 ,
u = 0 , v = 0 , T = T 0 ,
when t > 0,
u = L 1 u y + L 2 2 u y 2 , v = 0 , h [ T T 0 ] = K 1 T y + K 2 2 T y 2 ,   at   y = 0 ,
u = U ( x ) , T = T   as     y ,
where T 0 is a hot fluid at temperature, L 1 and L 2 denote the slip parameters of velocity, K 1 and K 2 are the jump parameters of temperature, h is a heat transfer coefficient.
For the sake of deriving a simplified model by converting partial differential equations into ordinary differential equations, a stream function ψ was introduced in this paper such that u = ψ y , v = ψ x . First, we scaled a set of points ( x , y , ψ , u , v ) :
Γ : x * = x e ε α 1 , y * = y e ε α 2 , y * = y e ε α 2 , u * = u e ε α 4 , v * = v e ε α 5 .
Then, we introduced Lie-group transformations to obtain the relationships between the parameters α i , i = 1 , , 5 . Based on the above Lie-group transformations, the stream function and similar parameter can be prescribed as follows:
η = a ν b f λ t y , ψ = a ν b f λ t x f ( η ) , θ = T T T 0 T , U = a x λ t .
For the convenience of calculation, we defined the following parameters:
γ = G r b f Re 2 b f , ν b f = μ b f ρ b f , G r b f = g β * b f ( T 0 T ) x 3 ν 2 b f , Re b f = U x ν b f ,
H a = σ n f B 0 2 x ρ b f U ,   Pr = ν b f α b f ,   Q = Q 0 x ( ρ C p ) b f U ,   N = 16 σ * T 3 3 κ * ,
N t = D T υ T ( T 0 T ) ,   B i = h x κ b f Re b f 1 2 ,   K = K n x Re b f 1 2 ,   K n x = L x .
The governing fundamental equations were transformed into the following third order ordinary differential equations after transformation:
f + A ( D f + f f + 1 f 2 ) + A B γ θ cos ( β π 2 ) H a C ( f 1 ) = 0 ,
( 1 + N F ) θ + 1 2 d Pr E θ η + Pr E θ f + N t Pr E θ 2 + Q Pr E θ = 0 .
Accordingly, boundary conditions were transferred into the following form:
f ( 0 ) = K f ( 0 ) + K 2 L f ( 0 ) ,   f ( 0 ) = 0 ,   f ( ) = 1 ,
θ ( 0 ) = 1 B i G θ ( 0 ) + K L B i G θ ( 0 ) + 1 ,   θ ( ) = 0 ,
where:
A = ( 1 ϕ e ) 2.5 [ 1 ϕ e + ϕ e ( ρ s p ρ b f ) ] ,   B = [ 1 ϕ e 1 ϕ e + ϕ e ( ρ s p ρ b f ) + ( β * s p β * b f ) ϕ e ( 1 ϕ e ) ( ρ b f ρ s p ) + ϕ e ] ,   C = ( 1 ϕ e ) 2.5 , D = λ a ,   E = ( ρ C p ) n f ( ρ C p ) b f κ b f κ n f ,   F = 1 κ n f ,   G = κ b f κ n f ,   d = λ 2 a .
When the fluid flows over a wedge, it will create skin friction drag C f . Nu is defined as the dimensionless temperature gradient of the fluid near the wall. It determines the strength of the convective heat transfer fluid. The local Nusselt number and skin friction coefficient are important physical parameters in the study of fluid flow and heat transfer, and they are defined as follows:
C f = 2 μ n f ρ b f U 2 u y | y = 0 ,   C f Re b f 1 2 = 2 ( 1 ϕ e ) 2.5 f ( 0 ) ,
N u x = x q w κ b f ( T 0 T ) ,   q w = κ n f T y | y = 0 ,   N u x Re b f 1 2 = 1 G θ ( 0 ) .

2.2. Application of HAM

Due to the strongly nonlinear and uncommon feature of the third order ordinary differential equations, we used the homotopy analysis method to get an approximate analytical solution.
Firstly, we defined that: λ 1 = K , λ 2 = K 2 L , ω 1 = 1 B i G , ω 2 = K L B i G , assuming the initial approximations were:
f 0 ( η ) = a 0 + a 1 η + a 2 η e η ,   θ 0 ( η ) = b 0 + b 1 e η .
Substituting the initial guess into the boundary conditions, we could get:
f 0 ( η ) = x + 1 3 λ 2 2 λ 1 1 x e x ,   θ 0 ( η ) = 1 ω 1 ω 2 + 1   e x .
Then, selecting the linear operator:
L f = d 3 f d η 3 + d 2 f d η 2 ,   L θ = d 2 θ d η 2 + d θ d η ,
and they both satisfied the following conditions:
L f [ C 1 + C 2 η + C 3 η e η ] = 0 ,   L θ [ C 4 + C 5 e η ] = 0 .
After that, the zero-order deformation equations could be constructed as follows:
( 1 p ) L f [ F ( η ; p ) f 0 ( η ) ] = p h f N f [ F ( η ; p ) ] ,
( 1 p ) θ f [ Θ ( η ; p ) θ 0 ( η ) ] = p h θ N θ [ Θ ( η ; p ) , F ( η ; p ) ] ,
where:
N f [ F ( η ; p ) , Θ ( η ; p ) ] = 3 F η 3 + A ( D F η + F η 2 F η 2 + 1 ( F η ) 2 ) + H A C ( 1 F η ) + A B γ Θ cos ( β π 2 )
N θ [ F ( η ; p ) , Θ ( η ; p ) ] = ( 1 + N F 1 ) 2 Θ η 2 + Pr E F Θ η + Q Pr E Θ + Pr E N t ( Θ η ) 2 + 1 2 d Pr η E Θ η .
Differentiating the zero-order deformation Equations (30)–(31) m-times with respect to p at p = 0 and then dividing the resulting expression by m!, we could attain high-order equations:
L f [ f m ( η ) χ m f m 1 ( η ) ] = h f R f m ( η ) ,
L θ [ θ m ( η ) χ m θ m 1 ( η ) ] = h θ R θ m ( η ) .
High order equations satisfied the boundary conditions:
f m ( 0 ) = f m ( + ) = 0 ,   f m ( 0 ) = λ 1 f m ( 0 ) + λ 2 f m ( 0 ) ,
θ m ( + ) = 0 , θ m ( 0 ) = ω 1 θ m ( 0 ) + ω 2 θ m ( 0 ) ,
where:
f m ( η ) = 1 m ! m F [ η ; p ] p m | p = 0 , θ m ( η ) = 1 m ! m Θ [ η ; p ] p m | p = 0 , F = f 0 ( η ) + m = 1 + f m ( η ) p m , Θ = θ 0 ( η ) + m = 1 + θ m ( η ) p m .

3. Results and Discussion

First of all, we determined the value of auxiliary parameter hm through the BVPh2.0 (Shanghai Jiao Tong University, shanghai, China; http://numericaltank.sjtu.edu.cn/BVPh2_0.htm) procedure package. Table 1 shows the error values of same auxiliary parameters and different auxiliary parameters.
Through the above form, we can obtain that the error reached 10−3 in the sixth order solution, which reached the standard of the engineering calculation. We defined the error as the following:
0 f + A ( D f + f f + 1 f 2 ) + A B γ θ cos ( β π 2 ) H a C ( f 1 )   d η .
With the increase of the order, we can determine that the error reached 10 6 and the rate of convergence of the two methods was similar, so we chose the same parameters in this paper.
The first step of HAM is to acquire the value of auxiliary parameter hm. Figure 1 and Figure 2 are plotted by the case that when other parameters are constant ( β = 0.15 , γ = 0.1 , H = 2.0 , Pr = 6.5 , N t = 0.5 , Q = 0.3 ), the range of hm was [−0.3, 0.3]. This paper selected hm = −0.1 for the next calculation.
After determining the value of hm, this paper analyzed the effects of various physical parameters on CfRe1/2bf and NuxRe−1/2bf.
Table 2 reflects the effect of radius on CfRe1/2bf and NuxRe−1/2bf, and this trend is consistent with its definition. The increase in radius led to the decrease in the solid volume fraction of nano-particles, which reduced heat transfer and the viscosity together with the resistance of the nanofluids.
In other words, the increase in radius caused the reduction of NuxRe−1/2bf and CfRe1/2bf.
Table 3 points to the effect of the thickness of the nanolayer on CfRe1/2bf and NuxRe−1/2bf, and this trend is also accordance with their definitions. The increase in the thickness of nanolayer led to an increase in the solid volume fraction of nano-particles, which added to heat transfer and the viscosity with the resistance of the nanofluids. Obviously, this gave rise to an increase in NuxRe−1/2bf and CfRe1/2bf.
Figure 3 and Figure 4 illustrate the effects of various solid particles and solid volume fraction on CfRe1/2bf and NuxRe−1/2bf, and this change conforms to its mathematical calculation formulas. Nanoparticles cause a change in the thermal conductivity of the mixture, which affects the flow of the fluid. The specific analysis is shown below. The increase in solid volume fraction causes a rise in the heat transfer and the viscosity as well as the resistance of the nanofluids. From the physical meaning of the Nusselt number we could determine that the increase in thermal conductivity led to the decrease in NuxRe−1/2bf. The increase in fluid resistance caused the increase in surface skin friction coefficient in the meantime. In addition, the effect of various solid particles on NuxRe−1/2bf was greater than that of CfRe1/2bf.
Figure 5 and Figure 6 show the effects of thermophoresis parameter Nt on CfRe1/2bf and NuxRe−1/2bf. Thermophoresis parameter Nt of nano-fluid in the fluid was used to measure the temperature distribution. When the thermophoresis parameter increased, the thermophoretic force helped to migrate the nanoparticles from the region of high temperature to low temperature area, which caused the temperature and the thermal efficiency to rise, resulting in an increase in NuxRe−1/2bf. Nt principally affected NuxRe−1/2bf, and the influence of Nt on CfRe1/2bf only presented a little increase trend. Researchers have shown that velocity slip and temperature jump have a crucial influence on the performance of the fluid at micro or nano scales.
Table 4 demonstrates the effect of the first order and the second order velocity slip and temperature jump on CfRe1/2bf and NuxRe−1/2bf. Velocity slip parameters were used to measure the sliding loss in fluid flow, and the sliding loss affected the velocity by way of the change of fluid flow. With the increase in the first order velocity slip, the wall velocity decreased, which caused the decrease in CfRe1/2bf and NuxRe−1/2bf. With the increase in the second order velocity slip, the result was the the opposite. The temperature jump parameters mainly affected the dimensionless temperature. With the increase in the first order temperature jump, the wall temperature increased. Therefore CfRe1/2bf and NuxRe−1/2bf both decreased. When the temperature of the second order jump is zero, it is a dividing line. When ω 2 0 , CfRe1/2bf decreased and NuxRe−1/2bf increased with the rise of the second order temperature jump; but while ω 2 > 0 , CfRe1/2bf rose and NuxRe−1/2bf fell.
Figure 7 demonstrates the effect of solid volume fraction on dimensionless temperature. When the volume fraction of nanoparticles rose, the thickness of the thermal boundary layer and the thermal conductivity increased, leading to the rise of the wall temperature of the nanofluid. The volume fraction of fluid mainly affected the heat transfer resistance, so the change in speed curve was not very obvious. Figure 8 demonstrates the effect of the thickness of the nanolayer on dimensionless temperature. The thickness of nanolayer primarily affected the solid volume fraction. The solid volume fraction rose with the increase in the thickness of the nanolayer, which gave rise to the decline in wall temperature.
Figure 9 shows the effect of the angle β of the wedge on dimensionless velocity. With the increase in the wedge angle, the thermal diffusion reduced, which caused the drop in the temperature of the fluid. This change increased the impact of buoyancy, resulting in the decrease in the driving force of the fluid, which caused the decrease in the fluid velocity. Figure 10 reflects B i on the dimensionless temperature. The Biot number reflects the distribution of temperature field in the object under unsteady heat conduction. The smaller Bi is, the smaller the internal thermal resistance and the bigger the external thermal resistance. With an increase in Bi, the ratio of the heat resistance increased, which gave rise to the reduction in wall temperature.

4. Conclusions

Considering the effects of the thickness of the nanolayer, temperature jump, and velocity slip, the model of unsteady nanofluid flow and heat transfer was modified in this paper. Meanwhile, both thermophoresis and Brownian motion were considered, and the differential equation was solved by HAM. Some of the main finds of this paper can be summarized as follows:
An increase in the thickness of the nanolayer causes an increase in the local Nussels number NuxRe−1/2bf and the skin friction coefficients CfRe1/2bf, together with the decrease in wall temperature. An enhancement of the thermophoresis parameter gives rise to an increase in NuxRe−1/2bf and CfRe1/2bf.
With an increase in λ 1 and ω 1 , CfRe1/2bf and NuxRe−1/2bf reduce. With an increase in the second order velocity slip parameter λ 2 , the result is the opposite. When the second order temperature jump parameter ω 2 0 , CfRe1/2bf decreases and NuxRe−1/2bf increases, with an enhancement of ω 2 ; but while ω 2 > 0 , CfRe1/2bf rises and NuxRe−1/2bf reduces.
An increase in solid volume fraction ϕ causes a decrease in the wall temperature. With an increase in the Biot number the dimensionless wall temperature drops. An increase in the angle of the wedge β leads to a decrease in wall velocity.

Author Contributions

J.Z. conducted the original research, modified the model, and contributed analysis tools. J.C. analyzed the data, simulated the modified model, and prepared original draft. J.Z. revised the manuscript.

Funding

This paper was supported by the National Natural Science Foundation of China (No. 11772046; No. 81870345).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. ASME Publ. Fed 1995, 231, 99–106. [Google Scholar]
  2. Eastman, J.A.; Choi, S.U.S.; Li, S. Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles. Appl. Phys. Lett. 2001, 78, 718–720. [Google Scholar] [CrossRef]
  3. Li, Y.J.; Zhou, J.E.; Tung, S.; Schneider, E.; Xi, S.Q. A review on development of nanofluid preparation and characterization. Powder Technol. 2009, 196, 89–101. [Google Scholar] [CrossRef]
  4. Yu, W.; Choi, S.U.S. The role of interfacial layers in the enhanced thermal conductivity of nanofluids: A renovated Hamilton–Crosser model. J. Nanopart. Res. 2004, 5, 167–171. [Google Scholar] [CrossRef]
  5. Azizian, R.; Doroodchi, E.; Moghtaderi, B. Influence of controlled aggregation on thermal conductivity of nanofluids. J. Heat Transf. 2016, 138, 43–55. [Google Scholar] [CrossRef]
  6. Bigdeli, M.B.; Fasano, M.; Cardellini, A.; Chiavazzo, E.; Asinari, P. A review on the heat and mass transfer phenomena in nanofluid coolants with special focus on automotive applications. Renew. Sustain. Energy Rev. 2016, 60, 1615–1633. [Google Scholar] [CrossRef] [Green Version]
  7. Fang, T.G. A note on the unsteady boundary layers over a flat plate. Int. J. Non Linear Mech. 2008, 43, 1007–1011. [Google Scholar] [CrossRef]
  8. Ullah, I.; Shafie, S.; Makinde, O.D.; Khan, I. Unsteady MHD Falkner-Skan flow of Casson nanofluid with generative/destructive chemical reaction. Chem. Eng. Sci. 2017, 172, 694–706. [Google Scholar] [CrossRef]
  9. Awan, S.E.; Khan, Z.A.; Awais, M.; Rehman, S.U.; Raja, M.A.Z. Numerical treatment for hydro-magnetic unsteady channel flow of nanofluid with heat transfer. Results Phys. 2018, 9, 1543–1554. [Google Scholar] [CrossRef]
  10. Jahan, S.; Sakidin, H.; Nazar, R.; Pop, I. Unsteady flow and heat transfer past a permeable stretching/shrinking sheet in a nanofluid: A revised model with stability and regression analyses. J. Mol. Liq. 2018, 261, 550–564. [Google Scholar] [CrossRef]
  11. Hayat, T.; Sajjad, R.; Muhammad, T. On MHD nonlinear stretching flow of Powell–Eyring nanomaterial. Results Phys. 2017, 7, 535–543. [Google Scholar] [CrossRef]
  12. Tian, X.Y.; Li, B.W.; Hu, Z.M. Convective stagnation point flow of a MHD non-Newtonian nanofluid towards a stretching plate. Int. J. Heat Mass Transf. 2018, 127, 768–780. [Google Scholar] [CrossRef]
  13. Hsiao, K.L. Combined Electrical MHD Heat Transfer Thermal Extrusion System Using Maxwell Fluid with Radiative and Viscous Dissipation Effects. Appl. Therm. Eng. 2017, 112, 1281–1288. [Google Scholar] [CrossRef]
  14. Babu, M.J.; Sandeep, N. Three-dimensional MHD slip flow of nanofluids over a slendering stretching sheet with thermophoresis and Brownian motion effects. Adv. Powder Technol. 2016, 27, 2039–2050. [Google Scholar] [CrossRef]
  15. Usman, M.; Soomro, F.A.; Haq, R.U.; Wang, W.; Defterli, O. Thermal and velocity slip effects on Casson nanofluid flow over an Inclined permeable stretching cylinder via collocation method. Int. J. Heat Mass Transf. 2018, 122, 1255–1263. [Google Scholar] [CrossRef]
  16. Abbas, N.; Saleem, S.; Nadeem, S.; Alderremy, A.A.; Khan, A.U. On stagnation point flow of a micro polar nanofluid past a circular cylinder with velocity and thermal slip. Results Phys. 2018, 9, 1224–1232. [Google Scholar] [CrossRef]
  17. Ramzan, M.; Chung, J.D.; Ullah, N. Partial slip effect in the flow of MHD micropolar nanofluid flow due to a rotating disk—A numerical approach. Results Phys. 2017, 7, 3557–3566. [Google Scholar] [CrossRef]
  18. Seth, G.S.; Mishra, M.K. Analysis of transient flow of MHD nanofluid past a non-linear stretching sheet considering Navier’s slip boundary condition. Adv. Powder Technol. 2017, 28, 375–384. [Google Scholar] [CrossRef]
  19. Thompson, P.A.; Troian, S.M. A general boundary condition for liquid flow at solid surfaces. Nature 1997, 389, 360–362. [Google Scholar] [CrossRef]
  20. Zhu, J.; Yang, D.; Zheng, L.C.; Zhang, X.X. Effects of second order velocity slip and nanoparticles migration on flow of Buongiorno nanofluid. Appl. Math. Lett. 2016, 52, 183–191. [Google Scholar] [CrossRef]
  21. Liu, Q.F.; Zhu, J.; Bandar, B.M.; Zheng, L.C. Hydromagnetic flow and heat transfer with various nanoparticles additives past a wedge with high order velocity slip and temperature jump. Can. J. Phys. 2017, 95, 440–449. [Google Scholar] [CrossRef]
  22. Sobamowo, M.G.; Jayesimi, L.O.; Waheed, M.A. Magnetohydrodynamic squeezing flow analysis of nanofluid under the effect of slip boundary conditions using variation of parameter method. Karbala Int. J. Mod. Sci. 2018, 4, 107–118. [Google Scholar] [CrossRef]
  23. Job, V.M.; Gunakala, S.R. Unsteady MHD free convection nanofluid flows within a wavy trapezoidal enclosure with viscous and Joule dissipation effects. Numer. Heat Transf. Appl. 2015, 69, 1–23. [Google Scholar] [CrossRef]
  24. Liao, S.J. Beyond perturbation: The basic concepts of the homotopy analysis method and its applications. Adv. Mech. 2008, 38, 1–33. [Google Scholar]
  25. Xie, H.Q.; Fujii, M.; Zhang, X. Effect of interfacial nanolayer on the effective thermal conductivity of nanoparticle-fluid mixture. Int. J. Heat Mass Transf. 2005, 48, 2926–2932. [Google Scholar] [CrossRef]
  26. Mahalakshmi, T.; Nithyadevi, N.; Oztop, H.F.; Abu-Hamdeh, N. MHD mixed convective heat transfer in a lid-driven enclosure filled with Ag-water nanofluid with center heater. Int. J. Mech. Sci. 2018, 142, 407–419. [Google Scholar] [CrossRef]
  27. Abdul Hakeem, A.K.; Raja, K.; Ganga, B.; Ganesh, N.V. Nanofluid slip flow through porous medium with elastic deformation and uniform heat source/sink effects. Comput. Therm. Sci. Int. J. 2019, 11, 269–283. [Google Scholar] [CrossRef]
  28. Rosseland, S. Astrophysik und Atom-Theoretische Grundlagen; Springer: Berlin, Germany, 1931; pp. 41–44. [Google Scholar]
  29. Kandasamy, R.; Muhaimin, I.; Khamis, A.B. Unsteady Hiemenz flow of Cu-nanofluid over a porous wedge in the presence of thermal stratification due to solar energy radiation: Lie group transformation. Int. J. Therm. Sci. 2013, 65, 196–205. [Google Scholar] [CrossRef]
Figure 1. hm- f ( 0 ) .
Figure 1. hm- f ( 0 ) .
Mathematics 07 01043 g001
Figure 2. hm- θ ( 0 ) .
Figure 2. hm- θ ( 0 ) .
Mathematics 07 01043 g002
Figure 3. The effect of ϕ on Nu.
Figure 3. The effect of ϕ on Nu.
Mathematics 07 01043 g003
Figure 4. The effect of ϕ on Cf.
Figure 4. The effect of ϕ on Cf.
Mathematics 07 01043 g004
Figure 5. Effect of Nt on Cf.
Figure 5. Effect of Nt on Cf.
Mathematics 07 01043 g005
Figure 6. The effect of Nt on Nu.
Figure 6. The effect of Nt on Nu.
Mathematics 07 01043 g006
Figure 7. The effect of ϕ on temperature.
Figure 7. The effect of ϕ on temperature.
Mathematics 07 01043 g007
Figure 8. The effect of δ on temperature.
Figure 8. The effect of δ on temperature.
Mathematics 07 01043 g008
Figure 9. The effect of β on velocity.
Figure 9. The effect of β on velocity.
Mathematics 07 01043 g009
Figure 10. The effect of B i on temperature.
Figure 10. The effect of B i on temperature.
Mathematics 07 01043 g010
Table 1. The error values of various orders.
Table 1. The error values of various orders.
OrderSame hmDifferent hm
20.098045973639846490.10971050115341624
40.016952284557340580.01987321691368215
60.001497186582817360.00111284953763120
80.000920062004606960.00103228583202702
100.000166529681125830.00023106953667965
120.000289378172321420.00018686850247371
140.000046253120773640.00003407739648355
160.000044111865158860.00003247720532823
180.000041383488059060.00003498110747688
200.000009133142038490.00000687359503864
Table 2. The effect of radius on CfRe1/2bf and NuxRe−1/2bf with β = 0.25 , λ 1 = 0.6 , λ 2 = 0.5 .
Table 2. The effect of radius on CfRe1/2bf and NuxRe−1/2bf with β = 0.25 , λ 1 = 0.6 , λ 2 = 0.5 .
RadiusNuxRe−1/2bfCfRe1/2bf
5−0.8915021.26755
6−0.8549121.24286
7−0.8308581.22730
8−0.8139501.21633
9−0.8014661.20772
10−0.7918841.20219
11−0.7843131.19735
12−0.7781921.19350
Table 3. The effect of the thickness of nanolayer on CfRe1/2bf and NuxRe−1/2bf.
Table 3. The effect of the thickness of nanolayer on CfRe1/2bf and NuxRe−1/2bf.
ThicknessNuxRe−1/2bfCfRe1/2bf
0−0.7212041.15796
1−0.7473871.17419
2−0.7781921.1935
3−0.813951.21633
4−0.8549121.24286
5−0.9011761.27453
6−0.9526441.31082
7−1.008931.35419
8−1.069351.40348
9−1.132841.46137
10−1.197961.52895
Table 4. The effect of particle and slip parameters on CfRe1/2bf and NuxRe−1/2bf.
Table 4. The effect of particle and slip parameters on CfRe1/2bf and NuxRe−1/2bf.
ParameterValueNuxRe−1/2bfCfRe1/2bf
λ 1 0−0.7397261.279140
0.5−0.7378900.987631
1−0.7367060.810298
1.5−0.7359530.678731
2−0.7353820.586453
λ 2 −1−0.7368510.829446
−0.5−0.7392021.184720
0−0.7444041.996120
0.3−0.7525063.308360
0.5−0.7675275.848600
ω 1 0−1.0602600.944702
1−0.4711350.944367
2−0.4078750.944307
3−0.3026040.944202
ω 2 −2−0.2105570.918785
−1−0.3720540.918769
−0.5−0.6039140.918768
0−1.8398100.916633
0.52.4058300.914590
10.6813840.918776
20.2941310.918839

Share and Cite

MDPI and ACS Style

Zhu, J.; Cao, J. Effects of Nanolayer and Second Order Slip on Unsteady Nanofluid Flow Past a Wedge. Mathematics 2019, 7, 1043. https://doi.org/10.3390/math7111043

AMA Style

Zhu J, Cao J. Effects of Nanolayer and Second Order Slip on Unsteady Nanofluid Flow Past a Wedge. Mathematics. 2019; 7(11):1043. https://doi.org/10.3390/math7111043

Chicago/Turabian Style

Zhu, Jing, and Jiahui Cao. 2019. "Effects of Nanolayer and Second Order Slip on Unsteady Nanofluid Flow Past a Wedge" Mathematics 7, no. 11: 1043. https://doi.org/10.3390/math7111043

APA Style

Zhu, J., & Cao, J. (2019). Effects of Nanolayer and Second Order Slip on Unsteady Nanofluid Flow Past a Wedge. Mathematics, 7(11), 1043. https://doi.org/10.3390/math7111043

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