Next Article in Journal
Development of a Staggered PCD End Mill for Carbon Fiber Reinforced Plastic
Next Article in Special Issue
Magnetohydrodynamic Nanoliquid Thin Film Sprayed on a Stretching Cylinder with Heat Transfer
Previous Article in Journal
Reconstruction to Sensor Measurements Based on a Correlation Model of Monitoring Data
Previous Article in Special Issue
Slip Flow and Heat Transfer of Nanofluids over a Porous Plate Embedded in a Porous Medium with Temperature Dependent Viscosity and Thermal Conductivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural Convection Flow of Fractional Nanofluids Over an Isothermal Vertical Plate with Thermal Radiation

1
Academy of Romanian Scientists, Bucharest 050094, Romania
2
Department of Theoretical Mechanics, Technical University of Iasi, Iasi 700050, Romania
3
Abdus Salam School of Mathematical Sciences, Government College University, Lahore 54000, Pakistan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2017, 7(3), 247; https://doi.org/10.3390/app7030247
Submission received: 27 November 2016 / Accepted: 1 March 2017 / Published: 3 March 2017
(This article belongs to the Special Issue Recent Developments of Nanofluids)

Abstract

:
The studies of classical nanofluids are restricted to models described by partial differential equations of integer order, and the memory effects are ignored. Fractional nanofluids, modeled by differential equations with Caputo time derivatives, are able to describe the influence of memory on the nanofluid behavior. In the present paper, heat and mass transfer characteristics of two water-based fractional nanofluids, containing nanoparticles of CuO and Ag, over an infinite vertical plate with a uniform temperature and thermal radiation, are analytically and graphically studied. Closed form solutions are determined for the dimensionless temperature and velocity fields, and the corresponding Nusselt number and skin friction coefficient. These solutions, presented in equivalent forms in terms of the Wright function or its fractional derivatives, have also been reduced to the known solutions of ordinary nanofluids. The influence of the fractional parameter on the temperature, velocity, Nusselt number, and skin friction coefficient, is graphically underlined and discussed. The enhancement of heat transfer in the natural convection flows is lower for fractional nanofluids, in comparison to ordinary nanofluids. In both cases, the fluid temperature increases for increasing values of the nanoparticle volume fraction.

1. Introduction

Natural convection flows have been extensively studied due to their multiple engineering applications. Such flows over an infinite plate are usually met in different engineering processes, including petroleum resource gas production, and geothermal reservoirs, thermal insulation, etc. (see [1,2,3]). The effects of thermal radiation are also important in geophysics, and geothermic, chemical, and ceramics processing, and they have been investigated by many researchers. A short presentation of the main results, up until 2007, is given by Ghosh and Beg [4], who studied the convective radiative heat transfer over a hot vertical surface in porous medium. Moreover, the effects of thermal radiation on nanofluid flows have been studied by many scholars. Mondal et al. [5] considered the unsteady magneto-hydrodynamic axi-symmetric stagnation-point flow over a shrinking sheet with Navier slip, and the temperature-dependent thermal conductivity. Magneto-hydrodynamic (MHD) flows of nanofluids, with radiation heat transfer over a flat plate with a variable heat flux and first-order chemical reaction, were studied by Zhang et al. [6]. A numerical study on Cu-water and Ag-water nanofluids, focusing on the radiation effects over a stretching sheet, was made by Abd Elazem et al. [7].
Nowadays, it is well known that one way of enhancing the thermal conductivity of fluids is by suspending the metallic particles, such as alumina, gold, copper, iron, or titanium, in fluids [8]. These particles, also called nanoparticles, have a diameter of less than 100 nm and the obtained solution is named nanofluid. The concept of a nanofluid seems to have been introduced by Choi [9], and based on his work, many researchers have focused their attention on the heat transfer in natural convection flows of nanofluids. Khan and Aziz [10] have studied the natural convection flow of a nanofluid over a vertical plate with a uniform surface heat flux. Turkyilmazoglu [11] provided exact solutions for the MHD slip flow of a nanofluid over a stretching/shrinking sheet, while Bachok et al. [12] studied the heat transfer characteristics on a moving plate in a nanofluid. Interesting exact solutions have also been obtained by Turkyilmazoglu and Pop [13], for the velocity and temperature fields corresponding to the natural convection flow of some nanofluids, past an infinite vertical plate with radiation effects. Radiation and magnetic effects on the natural convection flow of a nanofluid, past an infinite vertical plate with a heat source, have been studied by Mohankrishna et al. [14]. Ellahi [15] studied the effects of MHD and temperature-dependent viscosity on the flow of a non-Newtonian nanofluid in a pipe, while an analysis of the flow and heat transfer of water and ethylene glycol-based Cu nanoparticles between two parallel disks with suction/injection effects, has been provided by Rizwan Ul Haq et al. [16]. Of course, the list of such studies can continue, but we close it with some of the most interesting analytical and numerical results that have been obtained in [17,18,19,20,21,22,23,24,25,26,27].
However, none of these papers took into consideration the fractional derivatives in their governing equations, although the fractional models have been found to be quite flexible in describing the complex behavior of many materials. More recently, it seems that fractional partial differential equations may be used to describe some physical phenomena more accurately, when compared to the corresponding partial differential equations. Our interest here is to provide exact solutions for the temperature and velocity fields corresponding to the radiative natural convection flow of fractional nanofluids over an infinite vertical plate with heat and mass transfer, and to investigate the enhancement of heat transfer in such a flow, utilizing the fractional model. The associated skin friction coefficient and Nusselt number will be also determined. These solutions, which are presented in equivalent forms in terms of the Wright functions or its fractional derivatives, are reduced to similar solutions, corresponding to ordinary nanofluids [13]. Finally, the influence of the fractional parameter on the thermal and hydrodynamic response of physical interest, is graphically underlined and discussed.

2. Statement of the Problem

Let us consider the unsteady free convection flow and heat transfer of a nanofluid, modeled by the Caputo time-fractional derivative, past an infinite vertical plate situated in the ( x 1 , z 1 ) -plane of a fixed Cartesian coordinate system O x 1 y 1 z 1 . At the initial moment t 1 = 0 , the fluid and the plate are at rest, with a constant ambient temperature T . We also consider the radiation effect and assume the radiative heat flux to be applied, perpendicular to the plate. Since the plate is infinite, all of the physical quantities describing the fluid motion are functions of y 1 and t 1 . The fluid is a water-based nanofluid containing nano particles of CuO or Ag , whose thermo-physical properties are given in Table 1 [13,28].
Assuming a small difference between the fluid temperature T ( y 1 , t 1 ) and the stream temperature T , and adopting the Rosseland approximation [5], the radiative heat flux q r ( y 1 , t 1 ) can be linearized to:
q r ( y 1 , t 1 ) = 16 σ T 3 3 k T ( y 1 , t 1 ) y 1 ,
where σ is the Stefan-Boltzman constant and k is the mean absorption coefficient.
In the following, we consider the nanofluid model proposed by Tiwari and Das [29], and take into consideration the usual Boussinesq’s approximation. In this case, the governing equations can be written as [13]:
ρ n f u 1 ( y 1 , t 1 ) t 1 = μ n f 2 u 1 ( y 1 , t 1 ) y 1 2 + g ( ρ β ) n f [ T ( y 1 , t 1 ) T ] ; y 1 , t 1 > 0 ,
( ρ c p ) n f T ( y 1 , t 1 ) t 1 = k n f ( 1 + 16 σ T 3 3 k n f k ) 2 T ( y 1 , t 1 ) y 1 2 ; y 1 , t 1 > 0.
If no slipping exists between the fluid and the plate, the appropriate initial and boundary conditions are:
u 1 ( y 1 , 0 ) = 0 , T ( y 1 , 0 ) = T ; y 1 0 ,
u 1 ( 0 , t 1 ) = 0 , T ( 0 , t 1 ) = T w ; t 1 > 0 ,
u 1 ( y 1 , t 1 ) 0 , T ( y 1 , t 1 ) 0 as y 1 .
In the above relations, u 1 ( y 1 , t 1 ) is the fluid velocity in the x 1 -vertical direction, T w is the constant plate temperature ( T w > T or T w < T corresponds to the heated or cooled plate, respectively), g is the acceleration due to gravity, ρ n f is the density of the nanofluid, μ n f is the dynamic viscosity of the nanofluid, and β n f is the thermal expansion coefficient of the nanofluid. Their expressions, as well as the expression of ( ρ c p ) n f , are given by:
ρ n f = ( 1 φ ) ρ f + φ ρ s , μ n f = μ f ( 1 φ ) 2.5 ,
( ρ β ) n f = ( 1 φ ) ( ρ β ) f + φ ( ρ β ) s , ( ρ c p ) n f = ( 1 φ ) ( ρ c p ) f + φ ( ρ c p ) s
where φ is the nanoparticle volume fraction, ρ f is the density of the base fluid, ρ s is the density of the solid particle, and c p is the specific heat at constant pressure. The effective thermal conductivity of the nanofluid, corresponding to the Hamilton and Crosser model, is given by [28,30]:
k n f k f = k s + 2 k f 2 φ ( k f k s ) k s + 2 k f + φ ( k f k s ) ,
where k n f , k f , and k s are the thermal conductivities of the nanofluid, the fluid, and the solid particles, respectively.
Next, the non dimensional variables and functions are introduced as:
t = ν f L 2 t 1 , y = y 1 L , u = L ν f u 1 , θ = T T T w T , L = [ ν f 2 g β f ( T w T ) ] 1 / 3 ,
Equations (2)–(6) take simplified dimensionless forms, as follows:
u ( y , t ) t = 1 a 1 2 u ( y , t ) y 2 + a 2 θ ( y , t ) ; y , t > 0
θ ( y , t ) t = 1 a 3 2 θ ( y , t ) y 2 ; y , t > 0 ,
u ( y , 0 ) = 0 , θ ( y , 0 ) = 0 ; y 0 ,
u ( 0 , t ) = 0 , θ ( 0 , t ) = 1 ; t > 0 ,
u ( y , t ) 0 , θ ( y , t ) 0 as y ,
where:
a 1 = ( 1 φ ) 2.5 [ ( 1 φ ) + φ ρ s ρ f ] , a 2 = 1 φ + φ ( ρ β ) s ( ρ β ) f 1 φ + φ ρ s ρ f , a 3 = P r 1 φ + φ ( ρ c p ) s ( ρ c p ) f k n f k f + N r , P r = μ f c p f k f , N r = 16 σ T 3 3 k k f ,
where P r is the Prandtl number and N r is the radiation parameter.
The fractional model of the nanofluid is described by the fractional differential equations:
D t α c u ( y , t ) = 1 a 1 2 u ( y , t ) y 2 + a 2 θ ( y , t ) ; y , t > 0 ,
D t α c θ ( y , t ) = 1 a 3 2 θ ( y , t ) y 2 ; y , t > 0 ,
together with the initial and boundary conditions given by Equations (13)–(15). The operator D t α c represents the Caputo time-fractional derivative, defined as [31]:
D t α c u ( y , t ) = 1 Γ ( 1 α ) 0 t ( t s ) α u ( y , s ) s d s ; 0 α < 1.
The Caputo derivative is:
L { D t α c u ( y , t ) } ( q ) = q α u ¯ ( y , q ) q α 1 u ( y , 0 ) i f u ¯ ( y , q ) = L { u ( y , t ) } ( q )
and:
lim α 1 D t α c u ( y , t ) = u ( y , t ) t

3. Solution of the Problem

In order to determine the solution of the fractional partial differential Equations (17) and (18), with the initial and boundary conditions (13)–(15), the Laplace transform technique will be used. Equation (18) is not coupled to the momentum equation. Consequently, we shall firstly determine the temperature field.

3.1. Determination of the Temperature Field

Applying the Laplace transform to Equation (18), and bearing in mind the corresponding initial and boundary conditions, we find that:
a 3 q α θ ¯ ( y , q ) = 2 θ ¯ ( y , q ) y 2 ; y > 0 ,
where q is the transform parameter and the Laplace transform θ ¯ ( y , q ) of θ ( y , t ) has to satisfy the following conditions:
θ ¯ ( 0 , q ) = 1 q ; θ ¯ ( y , q ) 0   as   y .
The solution to the problems (22) and (23) is:
θ ¯ ( y , q ) = 1 q exp ( y a 3 q α ) ; y > 0.
Applying the inverse Laplace transform to Equation (24), and using Equation (A1) from the Appendix A, we find that:
θ ( y , t ) = Ψ ( 1 , α 2 ; y a 3 t α 2 ) for 0 < α 1 ,
where:
Ψ ( a , b ; z ) = n = 1 z n n ! Γ ( a n b ) ; b ( 0 , 1 ) ,
is the Wright function [32]. For α = 1 , Equation (25) becomes:
θ ( y , t ) = e r f c ( y a 3 2 t ) ; y , t > 0.
Of course, a simple analysis clearly shows that this result is in accordance with that obtained by Turkyilmazoglu and Pop [13], Equation (3.19).

3.2. Calculation of the Velocity Field

Applying the Laplace transform to Equation (17), and taking into consideration the corresponding initial and boundary conditions, we find that:
2 u ¯ ( y , q ) y 2 a 1 q α u ¯ ( y , q ) = a 1 a 2 1 q exp ( y a 3 q α ) ; y > 0
where the Laplace transform u ¯ ( y , q ) of u ( y , t ) has to satisfy the conditions:
u ¯ ( 0 , q ) = 0 ; u ¯ ( y , q ) 0   as y .
A particular solution of Equation (28) is:
u ¯ p ( y , q ) = a 1 a 2 a 1 a 3 1 q α + 1 exp ( y a 3 q α ) ,
while its general solution is:
u ¯ ( y , q ) = A exp ( y a 1 q α ) + B exp ( y a 1 q α ) + a 1 a 2 a 1 a 3 1 q α + 1 exp ( y a 3 q α ) .
Considering the conditions of (29), it results that:
u ¯ ( y , q ) = a 1 a 2 a 1 a 3 1 q α + 1 [ exp ( y a 3 q α ) exp ( y a 1 q α ) ] ; 0 < α < 1.
Applying the inverse Laplace transform to Equation (32), using the convolution theorem and the Equality (A1), we find that:
u ( y , t ) = a 1 a 2 a 1 a 3 1 Γ ( 1 β ) 0 t 1 ( t s ) β [ Ψ ( 1 , β 1 2 ; y a 3 s β 1 ) Ψ ( 1 , β 1 2 ; y a 1 s β 1 ) ] d s ,
where β = 1 α .
Now, using the identitiy (A2) from the Appendix A, we can present our solution in an interesting, but equivalent, form:
u ( y , t ) = a 1 a 2 a 1 a 3 { D t β c [ t Ψ ( 2 , β 1 2 ; y a 3 t β 1 ) ] D t β c [ t Ψ ( 2 , β 1 2 ; y a 1 t β 1 ) ] } ,
in terms of the Caputo derivative of the Wright functions.

3.3. Nusselt Number and Skin Friction

In order to determine the two entities of physical interest, namely the Nusselt number N u and the skin friction coefficient C f , we use the relations:
N u = L q w k f ( T w T ) = L k n f k f ( T w T ) T ( y 1 , t 1 ) y 1 | y 1 = 0 = k n f k f θ ( y , t ) y | y = 0 = k n f k f lim y 0 + L 1 { θ ¯ ( y , q ) y } .
C f = τ w ρ f ( ν f L ) 2 = μ n f ρ f ( ν f L ) 2 u 1 ( y 1 , t 1 ) y 1 | y 1 = 0 = μ n f μ f u ( y , t ) y | y = 0 = 1 ( 1 φ ) 2.5 lim y 0 + L 1 { u ¯ ( y , q ) y } ,
where q w is the constant heat flux from the surface of the plate and τ w is the skin friction or shear stress on the boundary.
Introducing Equations (24) and (32) into (35) and (36), respectively, we find that:
N u = k n f k f a 3 t α 2 Γ ( 1 α 2 ) , C f = a 1 a 2 a 1 + a 2 t α 2 Γ ( 1 + α 2 ) 1 ( 1 φ ) 2.5 .
Using the identity (A3) from the Appendix A, we also provide equivalent forms for N u and C f , namely:
N u = 2 k n f k f a 3 π D t 1 + α 2 c ( t 1 / 2 ) , C f = 1 ( 1 φ ) 2.5 a 1 a 2 a 1 + a 2 D t 1 α 2 c ( t ) ,
in terms of the Caputo derivatives of t 1 / 2 and t .

4. Validation

In order to bring to light the accuracy of the results that have been obtained, it is suitable to show that they are in accordance with similar solutions from the existing literature. For that, let us use β = 0 in Equation (34). It corresponds to α = 1 , and the solution corresponding to the same unsteady natural convection flow of ordinary nanofluids, has to be obtained. When the Caputo derivative of zero order is the identity operator, Equation (34) becomes:
u ( y , t ) = a 1 a 2 a 1 a 3 { [ t Ψ ( 2 , 1 2 ; y a 3 t ) ] [ t Ψ ( 2 , 1 2 ; y a 1 t ) ] } .
On the other hand (see also Equation (A2)):
e r f c ( y a 2 t ) = Ψ ( 1 , 1 2 ; y a t ) = t [ t Ψ ( 2 , 1 2 ; y a t ) ]
and then:
t Ψ ( 2 , 1 2 ; y a t ) = 0 t e r f c ( y a 2 s ) d s ,
where e r f c ( ) is the complementary error fucntion of Gauss.
Now, by introducing Equation (41) in (39), and by using Equation (A4), we get the velocity field as the simple form:
u ( y , t ) = a 1 a 2 a 1 a 3 { ( t + y 2 a 3 2 ) e r f c ( y a 3 2 t ) ( t + y 2 a 1 2 ) e r f c ( y a 1 2 t ) + y a 1 t π exp ( y 2 a 1 4 t ) y a 3 t π exp ( y 2 a 3 4 t ) } .
Finally, bearing in mind the notations of Turkyilmazoglu and Pop [13], as well as their rescaling relation from equality (2.11), it is easy to show that our solution (42) is identical to Equation (3.20), from [13].
With regards to the Nusselt number N u and the skin friction coefficient C f , we use α = 1 in Equation (38) and use Equation (A5). The expressions corresponding to the ordinary nanofluid are:
N u = k n f k f a 3 π t , C f = 1 ( 1 φ ) 2.5 2 a 1 a 2 a 1 + a 3 t π
As expected, by changing a 1 , a 3 , and t , by 1 a 1 , 1 a 3 , and τ 1 , respectively, we recover the solutions (3.21), from [13].

5. Numerical Results and Discussion

The natural convection flow of water-based fractional nanofluids over an infinite vertical plate with thermal radiation and a uniform temperature on the boundary, is analytically studied. Closed form solutions for the dimensionless temperature and velocity fields, and the two entities of physical interest, the Nusselt number and skin friction coefficient, are determined in equivalent forms, in terms of the Wright function or its fractional derivatives. It is worth pointing out that all of these solutions have been immediately reduced to the known solutions, based on the literature for ordinary nanofluids. A table containing the thermo-physical properties of copper oxide (CuO) and silver (Ag) is also included for later use.
In order to bring to light the influence of the fractional parameter on the heat and mass transfer in the natural convection flow of the above-mentioned fractional nanofluids, and therefore to obtain some physical insight into the present results, some numerical calculations have been carried out for different values of the fractional parameter α, radiation parameter Nr, and the nanoparticle value fraction ϕ. For comparison, the diagrams of dimensionless temperature and velocity fields, and the Nusselt number and skin friction coefficient corresponding to fractional nanofluids (for different values of the fractional parameter α ( 0 , 1 ) ) and those of ordinary nanofluids (when α = 1 ), are depicted in Figure 1, Figure 2, Figure 3 and Figure 4. As was expected, in all cases, the diagrams corresponding to the fractional nanofluids tend to superpose over those of ordinary nanofluids, when α 1 .
Profiles of the dimensionless temperature θ ( y , t ) against y are presented in Figure 1a–c, for the different values of the fractional parameter α and the CuO nanoparticle volume fraction φ . The fluid temperature, as it results from these figures, increases with respect to α , up to a critical value of y (less than 0.5), and then decreases. Consequently, in terms of the plate proximity, the heat transfer is stronger when the thermal boundary layer is thinner for fractional nanofluids, in comparison to the ordinary ones. An opposite trend appears at a later point in time. With respect to the volume fraction φ , the temperature is an increasing function for both ordinary and fractional nanofluids, and it smoothly decreases from the maximum value of one on the boundary, to the zero value far away from the plate. In Table 2, the temperature at different values of y and of the fractional parameter α , is given for nanofluids containing nanoparticles of CuO or Ag, and these values are in full accordance with those resulting from Figure 1. Furthermore, as results from this table, the Ag nanoparticles induce larger temperature values and a lower heat transfer to ordinary or fractional nanofluids, in comparison to CuO nanoparticles.
The effect of enhancing the heat transfer rate against ϕ with the fractional parameter α, for CuO-water and Ag-water nanofluids, is presented in Figure 2, in the absence or presence of thermal radiation. For both fractional nanofluids, the Nusselt number Nu is an increasing function with respect to α. It is also an almost linearly increasing function of ϕ. Further, the heat transfer rate for the CuO-water fractional nanofluid is always a little higher than that corresponding to Ag-water fractional nanofluid. However, the difference between them increases with increasing values of ϕ.
The influence of the fractional parameter α and of the nanoparticle volume fraction φ , on the dimensionless velocity u ( y , t ) against y, is brought to light by Figure 3a–c; which is also seen for a nanofluid with CuO nanoparticles. Near the plate, the nanofluid velocity increases up to a maximum value and then asymptotically decreases to the zero value for y values greater than 2.5, but it is a decreasing function with respect to the two parameters α and φ . Consequently, the boundary layer thickness is lower for ordinary nanofields, in comparison to fractional nanofluids. From a physical point of view, it means that the nanofluid viscosity decreases for increasing values of α or φ , and the ordinary nanofluids exhibit a stronger capacity in flow. This implies that the viscoelasticity strengthens the flow resistance with a decrease in the fractional parameter [33]. The variation of the skin friction coefficient C f , against the nanoparticle volume fraction ϕ, is presented in Figure 4, for the same nanofluids at different values of α and Nr, equal to zero or one.
The variation of the skin friction coefficient C f , against the nanoparticle volume fraction ϕ, is presented in Figure 4, for the same nanofluids at different values of α and Nr, equal to zero or one. In both cases, i.e., in the absence or presence of radiation, the skin friction coefficient is a decreasing function with respect to the fractional parameter α. Considering its variation with respect to ϕ, two different situations appear. In the absence of radiation (Figure 4b), it is an increasing function with respect to ϕ, and its values are always greater for Ag-water fractional nanofluid compared to CuO-water fractional nanofluid. In the presence of radiation, when Nr = 1 (Figure 4a), the shear stress on the plate decreases up to a critical value of ϕ (about 0.1), and then increases. Furthermore, up to this value of ϕ, it is smaller for Ag-water fractional nanofluid in comparison to CuO-water fractional nanofluid, and a reverse situation is presented at a later point in time. Throughout this study, the value of Pr has been taken as 6.067.

6. Conclusions

An analytical study of the natural convection flow of some water-based fractional nanofluids over an infinite vertical plate with thermal radiation and a uniform temperature on the boundary, is developed using the Caputo time-fractional derivative. The closed forms of solutions for the dimensionless temperature and velocity fields, and the corresponding Nusselt number and skin friction coefficient, are established in equivalent forms, in terms of the Wright function and its fractional derivatives. They have been reduced to the solutions obtained in [13], corresponding to ordinary nanofluids, when the fractional parameter α 1 .
In order to get some physical insight into the results, which have been obtained for CuO-water and Ag-water fractional nanofluids, some numerical calculations and graphical representations have been presented in Table 2 and Figure 1, Figure 2, Figure 3 and Figure 4, for different values of the fractional parameter α , radiative parameter N r , and the nanoparticles volume fraction φ . The main findings are:
The enhancement of heat transfer in the natural convection flows is lower for fractional nanofluids, in comparison to ordinary nanofluids. The thermal boundary layer is thicker for fractional nanofluids. In both cases, the fluid temperature increases and the heat transfer declines for increasing values of the nanoparticle volume fraction φ .
The flows of water-based fractional nanofluids are faster than the ordinary nanofluids.
A decrease in the fractional derivative parameter increases the thickness of the velocity boundary layer. From a physical point of view, it means that the nanofluid viscosity increases with the decreasing of α .
The dimensionless velocity of water-based fractional nanofluids, as well as that of ordinary nanofluids, is a decreasing function with respect to the nanoparticle volume fraction φ .
The skin friction coefficient in the natural convection flow is higher for fractional nanofluids, in comparison to ordinary nanofluids.
The enhancement of heat transfer in natural convection is stronger for Cu-water fractional nanofluids, when compared with Ag-water fractional nanofluids.
In the presence of radiation, Ag-water fractional/ordinary nanofluids achieve a lower skin friction coefficient near the plate, relative to Cu-water fractional/ordinary nanofluids.
Our study shows that the nanofluids described by fractional derivatives have a significantly different behavior than ordinary nanofluids. The same conclusions were obtained in the references [33,34]. Even if the models studied in these articles are different from the model studied here, the obtained results showed that the fractional parameter has a strong influence on the heat transfer process.

Acknowledgments

Waqas Ali Azhar is grateful to Abdus Salam School of Mathematical Sciences, Government College University, Lahore, Pakistan for supporting this scientific research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

L 1 { 1 q a exp ( c q b ) } = t a 1 Ψ ( a , b ; c t b ) ; b ( 0 , 1 )
t [ t Ψ ( 2 , b ; c t b ) ] = Ψ ( 1 , b ; c t b ) ; b ( 0 , 1 )
D t α c ( t γ ) = Γ ( γ + 1 ) Γ ( γ + 1 α ) t γ α
0 t e r f c ( y a 2 s ) d s = ( t + a y 2 2 ) e r f c ( y a 2 t ) y a π exp ( a y 2 4 t )
D t 1 c ( t 1 / 2 ) = 1 2 t , D t 1 / 2 c ( t ) = 2 π t

References

  1. Raptis, A. Unsteady free convection flow through a porous medium. Int. J. Eng. Sci. 1983, 21, 345–348. [Google Scholar] [CrossRef]
  2. Kim, S.J.; Vafai, K. Analysis of natural convection about a vertical plate embedded in a porous medium. Int. J. Heat Mass Transfer 1989, 32, 665–677. [Google Scholar] [CrossRef]
  3. Narasimha, K.R.; Pop, I. Transient free convection in a fluid saturated porous media with temperature dependent viscosity. Int. Commun. Heat Mass Transf. 1994, 21, 573–581. [Google Scholar]
  4. Ghosh, S.K.; Beg, O.A. Theoretical analysis of radiative effects on transient free convection heat transfer past a hot vertical surface in porous media. Nonlinear Anal. Model. Control 2008, 13, 419–432. [Google Scholar]
  5. Mondal, S.; Haroun, N.A.H.; Sibanda, P. The effects of thermal radiation on an unsteady MHD axisymmetric stagnation-point flow over a shrinking sheet in presence of temperature dependent thermal conductivity with Navier slip. PLoS ONE 2015, 10, e0138355. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, C.; Zheng, L.; Zhang, X.; Chen, G. MHD flow and radiation heat transfer of nanofluids in porous media with variable surface heat flux and chemical reaction. Appl. Math. Model. 2015, 39, 165–181. [Google Scholar] [CrossRef]
  7. Abd Elazem, N.Y.; Ebaid, A.; Aly, H.E. Radiation effect of MHD on Cu-water and Ag-water nanofluids flow over a stretching sheet: Numerical study. J. Appl. Computat. Math. 2015, 4. [Google Scholar] [CrossRef]
  8. Bianco, V.; Manca, O.; Nardini, S.; Vafai, K. Heat Transfer Enhancement with Nanofluids; CRC Press: London, UK, 2015. [Google Scholar]
  9. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition 66, San Francisco, CA, USA, 12–17 November 1995; pp. 99–105.
  10. Khan, W.A.; Aziz, A. Natural convection flow of a nanofluid over a vertical plate with uniform surface heat flux. Int. J. Therm. Sci. 2011, 50, 1207–1214. [Google Scholar] [CrossRef]
  11. Turkyilmazoglu, M. Exact analytical solutions for heat and mass transfer of MHD slip flow in nanofluids. Chem. Eng. Sci. 2012, 84, 182–187. [Google Scholar] [CrossRef]
  12. Bachok, N.; Ishak, A.; Pop, I. Flow and heat transfer characteristics on a moving plate in a nanofluid. Int. J. Heat Mass Transf. 2012, 55, 642–648. [Google Scholar] [CrossRef]
  13. Turkyilmazoglu, M.; Pop, I. Heat and mass transfer of unsteady natural convection flow of some nanofluids past a vertical infinite flat plate with radiation effect. Int. J. Heat Mass Transf. 2013, 59, 167–171. [Google Scholar] [CrossRef]
  14. Mohankrishna, P.; Sugunamma, V.; Sandeep, N. Radiation and magnetic field effects on unsteady natural convection flow of a nanofluid past an infinite vertical plate with heat source. Chem. Process. Eng. Res. 2014, 25, 39–52. [Google Scholar]
  15. Ellahi, R. The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe: Analytical solutions. Appl. Math. Model. 2013, 37, 1451–1457. [Google Scholar] [CrossRef]
  16. Rizwan, U.H.; Khan, Z.H.; Hussain, S.T.; Hammouch, Z. Flow and heat transfer analysis of water and ethylene glycol based Cu nanoparticles between two parallel disks with suction/injection effects. J. Mol. Liq. 2016, 221, 298–304. [Google Scholar] [CrossRef]
  17. Sheikholeslami, M.; Zaigham, Q.M.; Ellahi, R. Influence of induced magnetic field on free convection of nanofluid considering Koo-Kleinstreuer-Li (KKL) correlation. Appl. Sci. 2016, 324, 1–13. [Google Scholar] [CrossRef]
  18. Sekrani, G.; Poncet, S. Further investigation on laminar forced convection of nanofluid flows in a uniformly heated pipe using direct numerical simulations. Appl. Sci. 2016, 332, 1–24. [Google Scholar] [CrossRef]
  19. Rizwan, U.H.; Prabhakar, B.; Bandari, S. Impact of inclined Lorentz forces on Tangent hyperbolic nanofluid flow with zero normal flux of nanoparticles at the stretching sheet. Neural Comput. Appl. 2016. [Google Scholar] [CrossRef]
  20. Ibrahim, W.; Ul Haq, R. Magnetohydrodynamic (MHD) stagnation point flow of nanofluid past a stretching sheet with convective boundary condition. J. Braz. Soc. Mech. Sci. Eng. 2016, 38, 1155–1164. [Google Scholar] [CrossRef]
  21. Waqar, A.; Khan, C.R.; Rizwan, U.H. Heat Transfer Analysis of MHD Water Functionalized Carbon Nanotube Flow over a Static/Moving Wedge. J. Nanomater. 2015, 16, 934367. [Google Scholar] [CrossRef]
  22. Khan, W.A.; Khan, Z.H.; Haq, R.U. Flow and heat transfer of ferrofluids over a flat plate with uniform heat flux. Eur. Phys. J. Plus 2015, 130. [Google Scholar] [CrossRef]
  23. Rizwan, U.H.; Rajotia, D.; Noor, N.F.M. Thermophysical effects of water driven copper nanoparticles on MHD axisymmetric permeable shrinking sheet: Dual-nature study. Eur. Phys. J. 2016, 39. [Google Scholar] [CrossRef]
  24. Sheikholeslami, M.; Ellahi, R. Electrohydrodynamic nanofluid hydrothermal treatment in an enclosure with sinusoidal upper wall. Appl. Sci. 2015, 5, 294–306. [Google Scholar] [CrossRef]
  25. Rahman, S.U.; Ellahi, R.; Nadeem, S.; Zia, Q.M. Simultaneous effects of nanoparticles and slip on Jeffrey fluid through tapered artery with mild stenosis. J. Mol. Liq. 2016, 218, 484–493. [Google Scholar] [CrossRef]
  26. Sheikholeslami, M.; Ellahi, R. Three dimensional mesoscopic simulation of magnetic field effect on natural convection of nanofluid. Int. J. Heat Mass Transf. 2015, 89, 799–808. [Google Scholar] [CrossRef]
  27. Akbar, N.S.; Raza, M.; Ellahi, R. Copper oxide nanoparticles analysis with water as base fluid for peristaltic flow in permeable tube with heat transfer. Comput. Meth. Prog. Biomed. 2016, 130, 22–30. [Google Scholar] [CrossRef] [PubMed]
  28. 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]
  29. Tiwari, R.K.; Das, M.K. Heat transfer argument in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  30. Kakac, S.; Pramuanjaroenkij, A. Review of convective heat transfer enhancement with nanofluids. Int. J. Heat Mass Transf. 2009, 52, 3187–3196. [Google Scholar] [CrossRef]
  31. Povstenko, Y. Linear Fractional Diffusion-Wave Equation for Scientists and Engineers; Springer: Basel, Switzerland, 2015. [Google Scholar]
  32. Stankovic, B. On the function of E. M. Wright. Publ. I’Inst. Math. 1970, 10, 113–124. [Google Scholar]
  33. Cao, Z.; Zhao, J.; Wang, Z.; Liu, F.; Zheng, L. MHD flow and heat transfer of fractional Maxwell viscoelastic nanofluid over a moving plate. J. Mol. Liq. 2016, 222, 1121–1127. [Google Scholar] [CrossRef]
  34. Pan, M.; Zheng, L.; Liu, F.; Zhang, X. Modeling heat transport in nanofluids with stagnation point flow using fractional calculus. Appl. Math. Model. 2016, 40, 8974–8984. [Google Scholar] [CrossRef]
Figure 1. Variation of the dimensionless temperature θ(y,t) with the fractional parameter α and the CuO nanoparticle volume fraction ϕ. (a) for φ = 0 ; (b) for φ = 0.1 ; (c) for φ = 0.2 .
Figure 1. Variation of the dimensionless temperature θ(y,t) with the fractional parameter α and the CuO nanoparticle volume fraction ϕ. (a) for φ = 0 ; (b) for φ = 0.1 ; (c) for φ = 0.2 .
Applsci 07 00247 g001
Figure 2. Values of the Nusselt number against the nanoparticle volume fraction ϕ for the two water-based nanofluids. (a) for Nr = 1; (b) for Nr = 0.
Figure 2. Values of the Nusselt number against the nanoparticle volume fraction ϕ for the two water-based nanofluids. (a) for Nr = 1; (b) for Nr = 0.
Applsci 07 00247 g002
Figure 3. Variation of the velocity u(y,t) with the fractional parameter α and the nanoparticle volume fraction ϕ. (a) for φ = 0 ; (b) φ = 0.1 ; (c) for φ = 0.2 .
Figure 3. Variation of the velocity u(y,t) with the fractional parameter α and the nanoparticle volume fraction ϕ. (a) for φ = 0 ; (b) φ = 0.1 ; (c) for φ = 0.2 .
Applsci 07 00247 g003
Figure 4. Values of the skin friction coefficient C f against the nanoparticle volume fraction ϕ for two water-based fractional nanofluids. (a) for Nr = 1; (b) for Nr =0.
Figure 4. Values of the skin friction coefficient C f against the nanoparticle volume fraction ϕ for two water-based fractional nanofluids. (a) for Nr = 1; (b) for Nr =0.
Applsci 07 00247 g004
Table 1. Thermo-physical properties of water and nanoparticles.
Table 1. Thermo-physical properties of water and nanoparticles.
Basic Fluid/Nanoparticles ρ (Kg/m3)Cp(J/Kg K)k (W/m·K)β × 105 (1/K)
Pure water997.141790.61321
Copper oxide (CuO)6320531.876.51.80
Silver (Ag)10,5002354291.89
Table 2. Values of the dimensionless temperature for ϕ = 0.2, t = 0.5, and different values of the fractional parameters α and y.
Table 2. Values of the dimensionless temperature for ϕ = 0.2, t = 0.5, and different values of the fractional parameters α and y.
y CuOAg
α = 0.1 α = 0.5 α = 0.85 α = 1 α = 0.1 α = 0.5 α = 0.85 α = 1
011111111
0.10.863410.843480.848160.876890.866750.846940.850130.87864
0.20.743630.725360.759250.770860.749370.730830.764480.7768
0.30.640570.628030.654750.66090.6480.635040.662730.6691
0.40.551870.543410.556640.558680.560410.551660.566280.56886
0.50.475510.468510.468870.464760.484720.477660.479610.47635
0.60.409760.402180.391110.380350.41930.41190.40250.39275
0.70.353130.343790.323020.306080.362740.353760.334610.31872
0.80.304360.292770.264150.242130.313830.302720.275560.25448
0.90.262340.24850.213880.188210.271540.258210.224810.19985
10.226130.21030.171480.143710.234960.21960.18170.15432

Share and Cite

MDPI and ACS Style

Fetecau, C.; Vieru, D.; Azhar, W.A. Natural Convection Flow of Fractional Nanofluids Over an Isothermal Vertical Plate with Thermal Radiation. Appl. Sci. 2017, 7, 247. https://doi.org/10.3390/app7030247

AMA Style

Fetecau C, Vieru D, Azhar WA. Natural Convection Flow of Fractional Nanofluids Over an Isothermal Vertical Plate with Thermal Radiation. Applied Sciences. 2017; 7(3):247. https://doi.org/10.3390/app7030247

Chicago/Turabian Style

Fetecau, Constantin, Dumitru Vieru, and Waqas Ali Azhar. 2017. "Natural Convection Flow of Fractional Nanofluids Over an Isothermal Vertical Plate with Thermal Radiation" Applied Sciences 7, no. 3: 247. https://doi.org/10.3390/app7030247

APA Style

Fetecau, C., Vieru, D., & Azhar, W. A. (2017). Natural Convection Flow of Fractional Nanofluids Over an Isothermal Vertical Plate with Thermal Radiation. Applied Sciences, 7(3), 247. https://doi.org/10.3390/app7030247

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