Next Article in Journal
Field Navigation Using Fuzzy Elevation Maps Built with Local 3D Laser Scans
Next Article in Special Issue
The Rotating Flow of Magneto Hydrodynamic Carbon Nanotubes over a Stretching Sheet with the Impact of Non-Linear Thermal Radiation and Heat Generation/Absorption
Previous Article in Journal
Numerical Simulation of Early Age Cracking of Reinforced Concrete Bridge Decks with a Full-3D Multiscale and Multi-Chemo-Physical Integrated Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Analysis of Nanofluid Flow over a Curved Stretching Surface Suspended by Carbon Nanotubes with Internal Heat Generation

1
Department of Mathematics, Faculty of Sciences, HITEC University, Taxila Cantt 47080, Pakistan
2
Department of Mathematics, COMSATS Institute of Information Technology, Abbottabad 22060, Pakistan
3
Center for Research (CFR), University of Islamabad (UoI), Islamabad 44000, Pakistan
4
School of Mathematical Sciences, Faculty of Sciences and Technology, Universiti Kebangsaan Malaysia, Bangi 43600, Selangor, Malaysia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2018, 8(3), 395; https://doi.org/10.3390/app8030395
Submission received: 19 December 2017 / Revised: 7 February 2018 / Accepted: 12 February 2018 / Published: 8 March 2018
(This article belongs to the Special Issue Multi-Walled Carbon Nanotubes)

Abstract

:
We have investigated a two-dimensional radiative flow of a boundary layer nature. The fluid under consideration is carbon nanotube (CNT)-based nanofluid and it flows over a curved surface. The heat transfer through the flow is analyzed under the influence of internal heat generation. Water (base fluid) along with single or multi-walled carbon nanotubes is taken to compose the nanofluid. After introducing the suitable similarity variables, the consequent equations are reduced to a system of nonlinear ordinary differential equations. The solution to the system is computed by using the shooting method accompanied by Runge–Kutta–Fehlberg algorithm. Various parameters, emerging in the governing equations, influences the flow and heat transfer distribution. These changes are captured and portrayed in the form of graphs. The changes in local rate of heat transfer and skin friction coefficient are also enlisted. To ensure the correctness of applied numerical scheme, the results are compared with some already existing studies.

1. Introduction

Nanomaterials have revolutionized many industrial and household appliances. Their use is gaining considerable popularity, and the fields of science associated with them enjoy a similar repute. Nanotechnology and its consequent products have many useful applications in various fields of science and engineering, from enhancing thermal properties of traditional fluids to construction of an effective drug delivery system. Other examples include electronics, life sciences, chemical synthesis, fuel cells, medical sciences, microsystems including mechanical and electrical components, etc. Nanofluids are formed by the addition of nano-sized particles in base fluids, which traditionally have poor thermal properties. These additives work as boosters and the resulting substance, i.e., nanofluid, bears some remarkable thermal properties [1,2,3,4,5]. Several mathematical models describe the properties of nanofluids. Choi [6] was the pioneer who witnessed and announced the existence of the tiny particles and their effectiveness. He examined through several experiments that the addition of nanometer sized particles of various nature exceptionally enhance the thermal properties of base fluids. Later, Choi et al. [7] presented a mathematical model which provides some new insights related to the field of nanotechnology. Following in his footsteps, many researchers carried out several research activities. Subsequently, several models have been proposed for different physical and industrial problems. One of these models is presented by Buongiorno [8]. His model considers both thermophoretic and Brownian motion effects. Hamilton and Crosser [9] formulated another model which describes the effects of shapes of nanoparticles. Some other famous models include Maxwell’s [10] and Xue’s [11] models. Khan and Pop utilized Buongiorno’s concept to model a boundary layer flow over a stretching sheet [12]. There has been several efforts to use these models in different situations and geometries. Some of the most recent can be found in [13,14,15,16,17,18,19,20].
In 1991, Iijima [21] introduced a cylindrical member of fullerenes family named as carbon nanotubes (CNTs). They possess unique electronic, mechanical, chemical, optical, catalytic, adsorption, and transport properties, which make them very valuable for a wide range of applications [22,23]. They are mainly categorized in two modules, one Single-walled carbon nanotubes (SWCNTs) and the other multi-walled carbon nanotubes (MWCNTs). The term single refers to a single layer or wall that holds all the particles together, whereas the term multi-walled refers to a bunch of nested tubes whose diameter is continuously growing. It has been observed that the enhancement of thermal properties of base fluid using CNTs as nanoparticles is remarkable. Murshed et al. [24] showed that CNTs produce six times better thermal conductivity, at room temperature, as compared to other substances. The role of CNTs in enhancement of thermal properties of the base fluid has been highlighted by many scientists, in different scenarios and geometries [25,26,27,28].
During the last few years, researchers are showing a keen interest in understanding the flows of boundary layer nature over a dilating or a squeezing curved sheet. These types of flows bear commendable uses in different engineering and industrial processes such as production of papers, polymer sheet, manufacturing of rubber and plastic sheets, glass fibers, wire coating, food manufacturing, etc. In 1961, Sakiadis [29] was the first to come up with the idea of boundary layer flows. He presented the flow over a stretching surface. Crane [30] extended his idea for both linear and exponentially stretching sheets. He derived a closed form solution for the viscid flow caused by the stretching of a sheet. McLeod and Rajagopal [31] examined the uniqueness of the exact solution [30]. Since the pioneering work, flow problems related to a stretching surface under different situations have been analyzed by many scientists. Some of them can be found in References [32,33,34,35,36,37].
The energy crisis has become a major cause of concern in modern world. To meet a growing demand for energy, a new wave of innovative models is required that can replace the traditional ones. The resources for renewable energy, including solar and wind power, are now well placed to contribute to energy requirements in both mature and developing economies. Thermal radiation gains a great significance for such processes where the systems are running at very high temperatures, e.g., heating and cooling chambers, solar power technology, electrical power generation, thermal energy storage, nuclear power plants and many other industrial areas. Due to these important applications, many authors have contributed in studying the radiative flows of different fluids in various geometrical configurations [38,39].
Besides conduction and radiation, another means of heat transfer, i.e., convection, also plays a pivotal role in many real-life problems. Heat transfer via convection has been a core subject of many scientific investigations. Its role is very prominent in a wide range of energy-related applications such as automobile radiators, lubricants and coolants in various mechanical processes, etc. An asymptotic solution, by involving homotopy analysis method, to the stagnation point nanofluid flow (with a heat generation source) over a dilating permeable sheet was explained by Malvandi et al. [40]. Tsai et al. [41] focused on the influence of flow and heat transfer coupled with non-uniform heat source over an unsteady stretching surface. Pal [42] examined the radiative flow over an unsteady expanding surface with variable heat generating/consuming mechanism.
Previously, studies on the flows over flat surfaces were carried out using Cartesian coordinate systems. Recently, Sajid et al. [43] came up with an idea of curved surface and modeled the flow phenomena using a curvilinear coordinate system. Abbas with his coworkers [44] followed their approach and applied it to present a thermal analysis of the flow over a bent sheet. In 2016, Abbas et al. [45] studied the hydromagnetic and radiative slip flow of a nanofluid over a curved stretching surface along with heat generation. Most recently, a considerable interest has been shown by many authors in studying the flow over curved surfaces [46,47,48]. These recent efforts by renowned scholars have motivated us to carry out this work. We have examined the heat transfer in a radiative boundary layer flow, heat generation, of a carbon nanotube (CNT)-based nanofluid over a curved stretching sheet. Shooting method followed by Runge–Kutta–Fehlberg scheme is employed to obtain a numerical solution of the problem. Consequently, a comprehensive graphical description (highlighting the effects of various parameters on the velocity, pressure and temperature profiles) of the flow behavior is presented.

2. Formulation of the Problem

We consider a laminar flow of an incompressible nanofluid over a curved stretching sheet at r * = * (see Figure 1). The flow is steady and two dimensional. The nanofluid under consideration is composed of CNT nanoparticles added to a base fluid (water). Since we are dealing with a curved sheet that is being stretched in opposite direction by two forces of same magnitude, we use curvilinear coordinates ( r * ,   s * ) . The linear velocity by which the sheet is being stretched is taken to be 𝓊 ´ = 𝒶 ´ s * , where 𝒶 ´ is a positive constant. The temperature of the sheet is represented by T * w   ( T * w ( s * ) = A ´ ( s * / ) ; where A ´ is a constant). The free stream temperature is denoted by   T * and it is assumed that T * < T * w .
In context of the aforementioned coordinate scheme and the assumptions along with the boundary layer approximation, the equations representing the conservative balance for mass, momentum and energy are given as follows:
r * ( Υ v ´ ) + * 𝓊 ´ s * = 0
p * r * ρ ` n f 𝓊 ´ 2 Υ = 0
v ´ 𝓊 ´ r * + 𝓊 ´ Υ 𝓊 ´ s * + 𝓊 ´ v ´ Υ = 1 ρ ` n f * Υ p * s * + μ n f ρ ` n f { 2 𝓊 ´ r * 2 + 1 Υ 𝓊 ´ r * 𝓊 ´ Υ 2 }
v ´ T * r * + * 𝓊 ´ Υ T * s * = k n f ( ρ ` C p ) n f [ 2 T * r * 2 + 1 Υ T * r * 1 k n f Υ r * ( Υ ) q ˜ r * + Q ´ k n f ( T * T * ) ]
where Υ = r * + * . In the above set of equations, the velocities along s * and r * are denoted by the symbols 𝓊 ´ and v ´ , respectively. Moreover, p * represents the dimensional pressure, q ˜ r * describes the radiative heat flux and T * is the base fluid temperature. The term Q ´ arises due to heat generation and it represents the volumetric rate of heat generation through a source. Furthermore, ( ρ ` C p ) n f denotes volumetric heat capacity of the nanofluid and μ n f is its dynamic viscosity. Effective density and the thermal conductivity of the nanofluid are respectively symbolized by ρ ` n f , and k n f [11,27,49]. Mathematically:
υ n f = ( μ ρ ` 1 ) n f ,     μ n f = μ f ( 1 Φ ) 5 / 2 ,     α n f = ( k ρ ` C p ) n f ,     ρ ` n f ρ ` f = [ 1 ( 1 ρ ` C N T ρ ` f ) Φ ]
( ρ ` C p ) n f ( ρ ` C p ) f = [ 1 ( 1 ( ρ ` C p ) C N T ( ρ ` C p ) f ) Φ ]
k n f k f = ( ( 1 Φ ) ( k C N T k f ) + 2 Φ k C N T ln ( k C N T + k f 2 k f ) ( 1 Φ ) ( k C N T k f ) + 2 Φ k f ln ( k C N T + k f 2 k f ) )
where the base fluid viscosity is represented by μ f , Φ denotes the nanoparticle volume fraction and α n f is the thermal diffusivity. Moreover, ( k f , k C N T ) and ( ρ ` f , ρ ` C N T ) represent the thermal conductivities and densities of base fluid and CNT, respectively. The following table highlights the thermo-physical properties of pure water (base fluid) and CNTs.
The associated boundary conditions are imposed as follows:
A t   r * = 0 ,     𝓊 ´ = 𝓊 ´ w = 𝒶 ´ s * ,         v ´ = 0 ,         T * = T * w , 𝓊 ´ 0 ,       𝓊 ´ r * 0 ,     T * T *   as   r *
The radiative heat flux is defined by the means of Rosseland approximation [50] as:
q ˜ r * = 16 σ ´ T * 3 3 a R T * 4 r *
where Stefan–Boltzmann constant and Rosseland mean absorption coefficient, are σ ´ and   a R , respectively. The small changes in temperature diffusion are assumed in such a way that we may expand the Taylor series expansion of T * 4 about T * and ignore the terms of higher order. Consequently, we get
T * 4 ( 3 T * 4 4 T * 3 T * )
By substituting Equations Equations (9) and (10) into Equation (4), we get
v ´ T * r * + * 𝓊 ´ Υ T * s * = k n f ( ρ ` C p ) n f ( 1 + 16 σ ´ T * 3 3 a R k f ( k n f / k f ) ) [ 2 T * r * 2 + 1 Υ T * r * ] + Q ´ ( ρ ` C p ) n f ( T * T * )
Now by assuming R d = 16 σ ´ T * 3 / 3 a R k f as a radiation parameter, following Magyari and Pantokratoras [51], Equation (11) becomes
v ´ T * r * + * 𝓊 ´ Υ T * s * = υ f Λ 2 1 P r k n f k f ( 1 + R d ( k n f / k f ) ) [ 2 T * r * 2 + 1 Υ T * r * ] + Q ´ ( ρ ` C p ) n f ( T * T * )
where Prandtl number is given as P r = υ f / α f .
The dimensionless transformation variables, for the simplification of the flow equations, are defined as:
𝓊 ´ = 𝒶 ´ s * f ( ζ ) ,     v ´ = * Υ 𝒶 ´ υ f f ( ζ ) ,       ζ = 𝒶 ´ υ f r * ,       p ´ = ρ 𝒶 ´ 2 s * 2 P ( ζ ) , T * = T * + A ´ s * Θ ( ζ ) ,       Θ ( ζ ) = T * T * T * w T *  
The substitution of above equation implies an automatic satisfaction of the continuity Equation (1), while Equations (2), (3) and (12) are given as:
P ζ = f 2 ζ + κ
2 κ ζ + κ P = υ n f υ f [ f + f ζ + κ f ( ζ + κ ) 2 ] κ ζ + κ f 2 + κ ζ + κ f f + κ ( ζ + κ ) 2 f f
( 1 + R d ( k n f / k f ) ) ( Θ + Θ ζ + κ ) P r ( k n f / k f ) [ Λ 2 κ ζ + κ ( f Θ f Θ ) λ Θ ] = 0
where dimensionless radius of curvature and heat generation parameter are κ = * 𝒶 ´ / υ f and   λ = Q ´ / 𝒶 ´ ( ρ ` C p ) f , respectively. The classical energy equation can be recovered from Equation (6) by ignoring the radiative   ( R d = 0 ) and heat generation   ( λ = 0 ) effects. Moreover, it is noteworthy that by taking κ and in the absence of pressure gradient, Equation (15) gets converted to the classical problem of flat stretching sheet as discussed by Crane [30],
f f 2 + f f = 0 .
The implementation of Equation (13) reduces the boundary conditions to a dimensionless form given as:
f ( 0 ) = 0 ,       f ( 0 ) = 1 ,       Θ ( 0 ) = 1 , f ( ) = 0 ,       f ( ) = 0 ,       Θ ( ) = 0
We can eliminate P from Equations (14) and (15). Consequently,
f i v + 2 f ζ + κ f ( ζ + κ ) 2 + f ( ζ + κ ) 3 + Λ 1 [ κ ζ + κ ( f f f f ) κ ( ζ + κ ) 2 ( f 2 f f ) κ ( ζ + κ ) 3 f f ] = 0
where
Λ 1 = υ f υ n f = ( 1 Φ ) 5 / 2 [ 1 ( 1 ρ ` C N T ρ ` f ) Φ ] , Λ 2 = [ 1 ( 1 ( ρ ` C p ) C N T ( ρ ` C p ) f ) Φ ]                                                   }
After obtaining   f ( ζ ) , P ( ζ ) can be easily determined from Equation (15) as:
P ( ζ ) = ζ + κ 2 κ [ 1 Λ 1 { f + f ζ + κ f ( ζ + κ ) 2 } κ ζ + κ f 2 + κ ζ + κ f f + κ ( ζ + κ ) 2 f f ]
The skin-friction coefficient C f and the local Nusselt number N u s * in s * -direction are given by:
C f = ( τ r * s * ´ ) r * = 0 ρ ` f 𝓊 ´ w 2 ,         N u s * = s * q ˜ w / k f ( T * w T * )
where τ r * s * represents the wall shear stress and q ˜ w symbolizes the wall heat flux in s * -direction. Their mathematical expression are as follows:
τ r * s * = μ n f ( 𝓊 ´ r * 𝓊 ´ Υ ) | r * = 0 ,           q ˜ w = k n f ( 1 + 16 σ ´ T * 3 3 a R k f ( k n f / k f ) ) T * r * | r * = 0  
Using Equations (13), (21) and (22) becomes:
R e s * 1 / 2 C f = 1 ( 1 Φ ) 2.5 ( f ( 0 ) f ( 0 ) κ ) , R e s * 1 / 2 N u s * = k n f k f ( 1 + R d ( k n f / k f ) ) Θ ( 0 )
where the local Reynolds number is symbolized by R e s * = 𝒶 ´ s * 2 / υ f t .

3. Solution Procedure

The shooting method, followed by Runge–Kutta–Fehlberg scheme, has been implemented to solve the problem. The shooting method aids in transforming the given system into an initial value problem which can then be solved by Runge–Kutta–Fehlberg (RKF) method. Mathematical software Mathematica (Version 10, Wolfram: Computation Meets Knowledge, Champaign, IL, USA) has been used to solve the present problem by assuming the step size ∆ζ = 0.001, and the convergence criteria is fixed at 10 6 .

4. Results and Discussions

The size and shape of nanocomposites matter greatly for enhaning the thermal properties of the base fluid. Carbon nanotubes are long nanotube wires with distinctive size and shape and also possess significant physical properties. They are supposed to be much better than steel and Kevlar due to the material’s exceptional tensile strength and stiffness. Thus, the flow and heat transfer characteristics of single- and multi-walled carbon nanotubes over a curved stretching sheet have been investigated. The mass flow rate is assumed to be 0.02 kg / s . The core objective of this section is to graphically demonstrate the influence of various parameters like dimensionless radius of curvature κ , nanoparticle volume fraction Φ (ranging from 0 Φ 0.2 ), heat generation parameter λ and radiation parameter Rd on the velocity, temperature and pressure profiles. Table 1 displays the thermo-physical properties of the base fluid and the nanoparticles. Since water has been utilized as a base fluid, therefore, a fixed value for Prandtl number ( P r = 6.2 ) has been used.
Figure 2 visualizes the effects of dimensionless curvature on the velocity, temperature and pressure profiles. Figure 2a depicts the graphical description of the velocity profile for both SWCNTs and MWCNTs based nanofluid under the influence of dimensionless radius of curvature κ . One can clearly notice an increase in the velocity as well as the momentum boundary layer thickness, with increasing dimensionless curvature (i.e., decreasing κ ). As the fluid particles trace the curved path along the surface of the sheet, the curvature of the sheet, under the influence of centrifugal force, enhances the secondary flow of the fluid. This flow is relatively minor, which is then superimposed on the primary flow which as a result enhances the fluid velocity.
The effect of curvature seems to be small in the s * -direction, while, in the r * -direction, the impact is quite prominent which is mainly due to the centrifugal force acting towards the origin. Figure 2a also exhibits the role of dimensionless radius of curvature κ on temperature profile Θ ( ζ ) . The increment of the dimensionless curvature enhances the temperature of the fluid, i.e., the maximum amount of heat is transfer when the radius of curvature decreases. Figure 2b illustrates the influence of dimensionless curvature on the pressure profile. Clearly, a decline in dimensionless radius of curvature enhances the magnitude of pressure inside the boundary layer. However, far away from the boundary, the pressure remains almost negligible because the behavior of stream lines away from the sheet follows the same pattern of flow as in the case of flat stretching sheet. Furthermore, it is pertinent to mention that the pressure variations, near and away from the surface, seem insignificant in the case of flat stretching sheet (i.e., κ = 1000 ). On the other hand, the pressure does not remain constant when the curvature is employed into the sheet and, therefore, a certain variation has been perceived, especially inside the boundary layer.
Figure 3 has been placed to show the influence of nanoparticle volume fraction Φ on the velocity, temperature and pressure profiles. A certain rise in velocity has been observed when the nanoparticle volume fraction increases. Besides, it has been observed that the SWCNTs have a slightly lower velocity than MWCNTs, which is due to the higher density values for SWCNTs with which the resistance within the fluid increases and consequently the dimensionless horizontal velocity component, i.e., f ( ζ ) , for SWCNTs decreases. Figure 3a also shows the variations in temperature Θ ( ζ ) with increasing Φ . The temperature of the fluid rises and consequently, the thermal boundary layer increases. Since the carbon nanotubes bear high thermal conductivity and low specific heat as compared to the base fluid, the inclusion of sufficient nanoparticles enhances the temperature of the fluid quite significantly. Moreover, the temperature for MWCNT-nanoparticles exhibits lower values as compared to SWCNT-nanoparticles. Figure 3b depicts the behavior of pressure distribution for distinct values of solid volume fraction Φ . A decline has been noticed in the magnitude of pressure distribution with an increasing Φ .
The results in Figure 4 highlight the variations in temperature for increasing values of radiative factor   R d . As expected, because of radiation, nanofluid temperature raises quite significantly. The reason behind is the decline in mean absorption coefficient and as a result the rate with which radiative heat transfers seems to be increasing at every point distant from the sheet. It has also been observed that the temperature of SWCNTs-nanoparticles increases more than MWCNTs nanoparticles because of the less thermal conductivity of the MWCNTs nanoparticles. Figure 5 depicts the behavior of heat generation parameter λ on the temperature profile. The increasing λ > 0 causes an upsurge in temperature as well as in thermal boundary layer thickness for both SWCNTs and MWCNTs nanoparticle. The heat generation parameter usually involves the heat generation coefficient and, thus, the increment of heat generation parameters leads to an increase in heat generation coefficient which means that the temperature at the surface is higher than the free stream temperature and, therefore, by transferring the heat from the sheet to the fluid, the rise in the temperature is quite obvious. Moreover, the increment of heat absorption parameter λ < 0 implies a decline in temperature in the region adjacent to the wall for both SWCNTs and MWCNTs based nanofluid.
Figure 6 describes the influential behavior of various parameters on skin friction coefficient. The increment in nanoparticle volume fraction considerably enhances the magnitude of skin friction coefficient, as seen in Figure 6a. The reason is the higher volume fraction of nanoparticles, in terms of composite carbon chains, which plays a prominent role within the fluid as well as on the surface and as a result increases the surface friction in both directions. Furthermore, it has been observed that the SWCNTs remain on the higher side because the higher density value of SWCNTs as compared to MWCNTs, therefore, the resistance within the fluid, seems to be more significant in the case of SWCNTs and thus exhibits the higher skin friction. Figure 6b depicts the impact of dimensionless radius of curvature κ on skin friction coefficient for both SW- and MW-carbon nanotubes. A decline in skin friction coefficient has been observed with the increasing values of dimensionless radius of curvature κ , which indicates that the increment in κ decreases the curvature of the surface and, therefore, offers a less resistance to the fluid particles which eventually decay the skin friction coefficient.
Figure 7 illustrates the effects of various parameters such as solid volume fraction Φ , heat generation parameter λ , radiation parameter R d and dimensionless radius of curvature κ . on the rate of heat transfer for both single walled and multi walled carbon nanotubes. In Figure 7a, one can clearly observe an enhancement in the local heat flux rate with an increasing nanoparticle volume fraction Φ for both SW- and MW-carbon nanotubes. Nusselt number is usually a product of thermal conductivity ratio and temperature gradient.The temperature gradient experiences a decline due to the inclusion of nanoparticles and it is much smaller than the thermal conductivity ratio, which as a result enhances the magnitude of Nusselt number. Furthermore, it has also been noticed in Figure 7a that the local heat flux experiences a decline with the increasing heat generation parameter λ > 0 . Since the increment of heat generation parameter significantly raises the temperature of the fluid, the decrement in the rate of heat transfer is obvious. The impact of radiation parameter Rd on the magnitude of the local Nusselt number is presented in Figure 7b. It has been noticed that the rate of heat flux is an increasing function of R d . The increment of radiation parameter R d drastically enhances the heat transfer coefficient which as a result shows an increment in heat transfer rate. Moreover, Figure 7b also describes the effect of dimensionless radius of curvature κ on the rate of heat flux. The increment in the magnitude of local Nusselt number has been perceived with an increasing dimensionless curvature (i.e., decreasing κ ). Since the temperature of the fluid rises with an increasing dimensionless curvature, rate of heat flux seems to be significant with decreasing κ .
Table 2 has been organized to see the variations in thermal and physical properties of base fluid with the dispersion of carbon nanotubes. It has been clearly noticed that the inclusion of both types of carbon nanotubes significantly enhances the thermal conductivity of the base fluid. The density of the fluid also increases while the specific heat experiences a clear decline. Table 3 and Table 4 has been organized to discover the impact of various embedded parameters on local skin friction and local Nusselt number for both SW- and MW-carbon nanotubes. Table 3 illustrates the influence of both solid volume fraction Φ and dimensionless radius of curvature κ on skin friction coefficient C f . In addition, comparison with the previously existing results of Abbas et al. [44] has been presented. It has been witnessed that the skin friction coefficient experiences a decline because of increasing dimensionless radius of curvature κ which indicates that a larger drag force is required to drag the fluid over the curved stretching sheet. Table 4 provides the numerical values for skin friction coefficient R e s * 1 / 2 C f and the magnitude of local Nusselt number R e s * 1 / 2 N u s *   for numerous values of solid volume fraction Φ and heat generation parameter λ under the assumption that the dimensionless radius of curvature κ i.e., κ = 1000 (for flat dilating sheet). It has been observed that the skin friction coefficient increases with increasing values of Φ , which can also be observed in Table 3. Moreover, less resistance has been offered by MWCNTs as compared to SWCNTs. In addition, the magnitude of local Nusselt number seems to be an increasing function of Φ for both SWCNTs and MWCNTs. The reason is the supremacy of CNTs over traditional base fluid, i.e., water, which is primarily due to the higher thermal conductivity and low specific heat. Thus, the ability to conduct heat significantly improves the rate of heat flux at the sheet. This proves the fact that the addition of carbon nanotubes enhances the local heat flux. On the other hand, the local rate of heat flux experiences a decrement with the increasing values of λ for CNTs based nanofluid as well as traditional base fluid. Finally, the validity of our results for a dimensionless radius of curvature κ , i.e., κ = 1000 (for flat dilating sheet), is shown in Table 5. An excellent agreement has been found between obtained results and previously published results of Ali [52], Grubka and Bobba [53] and Ishak et al. [54] in the absence of radiation parameter R d , heat generation parameter λ and carbon nanotubes Φ .

5. Conclusions

A time independent, steady, incompressible boundary layer flow of CNT (single-walled and multi-walled carbon nanotubes) based nanofluid over a curved dilating sheet has been considered. Impact of different physical parameters have enabled us to draw the following conclusions:
  • The increment in dimensionless curvature (i.e., decreasing κ ) and nanoparticle volume fraction Φ significantly enhances the fluid velocity. Moreover, velocity dominates for multi-walled carbon nanotubes.
  • An increment in temperature along with thermal boundary layer has been observed for increasing dimensionless curvature (i.e., decreasing κ ), nanoparticle volume fraction Φ and radiation parameter Rd. Furthermore, SWCNTs remains on the higher side as compared to MWCNTs.
  • An upsurge in temperature has been perceived for increasing heat generation parameter λ > 0 while the opposite behavior has been detected for increasing heat absorption parameter λ < 0 .
  • A decline is observed in the magnitude of pressure distribution inside the boundary layer as the dimensionless radius of curvature increases while a rise in nanoparticle volume fraction depicts a clear enhancement in it.
  • Skin friction enhances for higher values of solid volume fraction Φ while dimensionless radius of curvature κ implies a reduction in skin friction coefficient for both SWCNTs and MWCNTs.
  • The magnitude of local heat flux rate increases for growing values of dimensionless curvature (i.e., decreasing κ ), nanoparticle volume fraction Φ and radiation parameter Rd while the behavior is opposite for increasing heat generation parameter λ.

Acknowledgments

The work here is supported by grant: UKM-2017-064. Moreover, first author Syed Tauseef Mohyud-Din is also thankful to Chairman Bahria Town/ Patron and Chairman of FAIRE; Chief Executive FAIRE; Administration of University of Islamabad (a project of Bahria Town) for establishment of Center for Research (CFR) and the provision of conducive research environment.

Author Contributions

Fitnat Saba and Naveed Ahmed developed the problem and its MATHEMATICA code. Umar Khan in collaboration with Syed Tauseef Mohyud-Din, did the literature review, developed and implemented the computer code, and interpreted the subsequently obtained results. Saqib Hussain did the editing and removed the grammatical mistakes. Maslina Darus in consultation of rest of the Authors, re-confirmed the credibility of obtained solutions. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A ´ Constant
κ Dimensionless radius of curvature
k Thermal conductivity, W / mK
R d Radiation parameter
C p Specific heat at constant pressure,   J / kg · K
p ´ Pressure
P Dimensionless pressure
P r Prandtl number
Q ´ Volumetric rate of heat source
q ˜ r * Radiative heat flux, W / m 2
q ˜ w Wall heat flux,   W / m 2
* Radius of curvature
R e s * Local Reynold number
N u s * Local Nusselt number
T * Local fluid temperature,   K
T * w Surface temperature, K
T * Free stream temperature,   K
𝓊 ´ s * component of velocity, m / s
v ´ r * component of velocity, m / s
Greek Symbols
Φ Solid volume fraction of CNT
ζ Similarity variable
μ Dynamic viscosity,   N · s / m 2
υ Kinematic viscosity,   m 2 / s
ρ ` Density, kg / m 3
α Thermal diffusivity, m 2 / s
ρ ` C p Heat capacitance
Θ Dimensionless temperature
λ Heat generation parameter
τ r * s * Wall shear stress
σ ´ Stefan–Boltzmann constant
a R Mean absorption coefficient
Subscripts
n f Nanofluid
f Base fluid
CNT Carbon nanotube
SW Single walled
MW Multi walled

References

  1. Nnanna, A.G.A. Experimental model of temperature-driven nanofluids. J. Heat Transf. 2007, 129, 697–704. [Google Scholar] [CrossRef]
  2. Gu, B.; Hou, B.; Lu, Z.; Wang, Z.; Chen, S. Thermal conductivity of nanofluids containing high aspect ratio fillers. Int. J. Heat Mass Transf. 2013, 64, 108–114. [Google Scholar] [CrossRef]
  3. Milanese, M.; Iacobazzi, F.; Colangelo, G.; De Risi, A. An investigation of layering phenomenon at the liquid-solid interface in Cu and CuO based nanofluids. Int. J. Heat Mass Transf. 2016, 103, 564–571. [Google Scholar] [CrossRef]
  4. Iacobazzi, F.; Milanese, M.; Colangelo, G.; Lomascolo, M.; De Risi, A. An explanation of the Al2O3 nanofluid thermal conductivity based on the phonon theory of liquid. Energy 2016, 116, 786–794. [Google Scholar] [CrossRef]
  5. Colangelo, G.; Milanese, M.; De Risi, A. Numerical simulation of thermal efficiency of an innovative Al2O3 nanofluid solar thermal collector: Influence of nanoparticles concentration. Therm. Sci. 2016, 2016, 2769–2779. [Google Scholar] [CrossRef]
  6. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. Dev. Appl. Non Newton. Flows 1995, 231, 99–105. [Google Scholar]
  7. Choi, S.U.S.; Zhang, Z.G.; Yu, W.; Lockwood, F.E.; Grulke, E.A. Anomalous thermal conductivity enhancement in nanotube suspensions. Appl. Phys. Lett. 2001, 79, 2252–2254. [Google Scholar] [CrossRef]
  8. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  9. Hamilton, R.L.; Crosser, O.K. Thermal conductivity of heterogeneous two-component systems. Ind. Eng. Chem. Fundam. 1962, 1, 187–191. [Google Scholar] [CrossRef]
  10. Maxwell, J.C. A Treatise on Electricity and Magnetism, 3rd ed.; Oxford University Press: Oxford, UK, 1904. [Google Scholar]
  11. Xue, Q.Z. Model for thermal conductivity of carbon nanotube-based composites. Phys. B Condens. Matter 2005, 368, 302–307. [Google Scholar] [CrossRef]
  12. Khan, W.A.; Pop, I. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53, 2477–2483. [Google Scholar] [CrossRef]
  13. Kandelousi, M.S.; Ellahi, R. Simulation of Ferrofluid Flow for Magnetic Drug Targeting Using the Lattice Boltzmann Method. J. Z. Naturforschung A 2015, 70, 115–124. [Google Scholar] [CrossRef]
  14. Akbar, N.S.; Raza, M.; Ellahi, R. Influence of induced magnetic field and heat flux with the suspension of carbon nanotubes for the peristaltic flow in a permeable channel. J. Magn. Magn. Mater. 2015, 381, 405–415. [Google Scholar] [CrossRef]
  15. Sheikholeslami, M.; Ganji, D.D.; Younus Javed, M.; Ellahi, R. Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model. J. Magn. Magn. Mater. 2015, 374, 36–43. [Google Scholar] [CrossRef]
  16. Rashidi, S.; Dehghan, M.; Ellahi, R.; Riaz, M.; Jamal-Abad, M.T. Study of stream wise transverse magnetic fluid flow with heat transfer around an obstacle embedded in a porous medium. J. Magn. Magn. Mater. 2015, 378, 128–137. [Google Scholar] [CrossRef]
  17. Mabood, F.; Khan, W.A.; Ismail, A.I.M. 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]
  18. Mohyud-Din, S.T.; Zaidi, Z.A.; Khan, U.; Ahmed, N. On heat and mass transfer analysis for the flow of a nanofluid between rotating parallel plates. Aerosp. Sci. Technol. 2015, 46, 514–522. [Google Scholar] [CrossRef]
  19. Mohyud-Din, S.T.; Khan, U.; Ahmed, N.; Hassan, S. Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels. Appl. Sci. 2015, 5, 1639–1664. [Google Scholar] [CrossRef]
  20. Khan, U.; Ahmed, N.; Mohyud-Din, S.T. Thermo-diffusion, diffusion-thermo and chemical reaction effects on MHD flow of viscous fluid in divergent and convergent channels. Chem. Eng. Sci. 2016, 141, 17–27. [Google Scholar] [CrossRef]
  21. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  22. Saito, R.; Dresselhaus, G.; Dresselhaus, M.S. Physical Properties of Carbon Nanotubes; Imperial College Press: Singapore, 2001. [Google Scholar]
  23. Endo, M.; Hayashi, T.; Kim, Y.A.; Terrones, M.; Dresselhaus, M.S. Applications of carbon nanotubes in the twenty–first century. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2004, 362, 2223–2238. [Google Scholar] [CrossRef] [PubMed]
  24. Murshed, S.M.S.; Nieto de Castro, C.A.; Lourenço, M.J.V.; Lopes, M.L.M.; Santos, F.J.V. A review of boiling and convective heat transfer with nanofluids. Renew. Sustain. Energy Rev. 2011, 15, 2342–2354. [Google Scholar] [CrossRef]
  25. Akbar, N.S.; Butt, A.W. CNT suspended nanofluid analysis in a flexible tube with ciliated walls. Eur. Phys. J. Plus 2014, 129, 174. [Google Scholar] [CrossRef]
  26. Ul Haq, R.; Nadeem, S.; Khan, Z.H.; Noor, N.F.M. Convective heat transfer in MHD slip flow over a stretching surface in the presence of carbon nanotubes. Phys. B Condens. Matter 2015, 457, 40–47. [Google Scholar] [CrossRef]
  27. Ahmed, N.; Mohyud-Din, S.T.; Hassan, S.M. Flow and heat transfer of nanofluid in an asymmetric channel with expanding and contracting walls suspended by carbon nanotubes: A numerical investigation. Aerosp. Sci. Technol. 2016, 48, 53–60. [Google Scholar] [CrossRef]
  28. Khan, U.; Ahmed, N.; Mohyud-Din, S.T. Heat transfer effects on carbon nanotubes suspended nanofluid flow in a channel with non-parallel walls under the effect of velocity slip boundary condition: A numerical study. Neural Comput. Appl. 2017, 28, 37–46. [Google Scholar] [CrossRef]
  29. Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7, 26–28. [Google Scholar] [CrossRef]
  30. Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. 1970, 21, 645–647. [Google Scholar] [CrossRef]
  31. McLeod, J.B.; Rajagopal, K.R. On the uniqueness of flow of a Navier-Stokes fluid due to a stretching boundary. Arch. Ration. Mech. Anal. 1987, 98, 385–393. [Google Scholar] [CrossRef]
  32. Gupta, P.S.; Gupta, A.S. Heat and mass transfer on a stretching sheet with suction or blowing. Can. J. Chem. Eng. 1977, 55, 744–746. [Google Scholar] [CrossRef]
  33. Brady, J.F.; Acrivos, A. Steady flow in a channel or tube with an accelerating surface velocity. An exact solution to the Navier-Stokes equations with reverse flow. J. Fluid Mech. 1981, 112, 127–150. [Google Scholar] [CrossRef]
  34. Wang, C.Y. The three-dimensional flow due to a stretching flat surface. Phys. Fluids 1984, 27, 1915–1917. [Google Scholar] [CrossRef]
  35. Wang, C.Y. Fluid flow due to a stretching cylinder. Phys. Fluids 1988, 31, 466–468. [Google Scholar] [CrossRef]
  36. Chamkha, A.J.; Aly, A.M.; Mansour, M.A. Similarity solution for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects. Chem. Eng. Commun. 2010, 197, 846–858. [Google Scholar] [CrossRef]
  37. Bhattacharyya, K.; Uddin, M.S. Reactive Solute Diffusion in Boundary Layer Flow through a Porous Medium over a Permeable Flat Plate with Power-Law Variation in Surface Concentration. J. Eng. 2013, 2013, 1–7. [Google Scholar] [CrossRef]
  38. Rashidi, M.M.; Vishnu Ganesh, N.; Abdul Hakeem, A.K.; Ganga, B. Buoyancy effect on MHD flow of nanofluid over a stretching sheet in the presence of thermal radiation. J. Mol. Liq. 2014, 198, 234–238. [Google Scholar] [CrossRef]
  39. Ul Haq, R.; Nadeem, S.; Khan, Z.H.; Akbar, N.S. Thermal radiation and slip effects on MHD stagnation point flow of nanofluid over a stretching sheet. Phys. E Low Dimens. Syst. Nanostruct. 2015, 65, 17–23. [Google Scholar] [CrossRef]
  40. Malvandi, A.; Hedayati, F.; Nobari, M.R.H. An HAM Analysis of Stagnation-Point Flow of a Nanofluid over a Porous Stretching Sheet with Heat Generation. J. Appl. Fluid Mech. 2014, 7, 135–145. [Google Scholar]
  41. Tsai, R.; Huang, K.H.; Huang, J.S. Flow and heat transfer over an unsteady stretching surface with non-uniform heat source. Int. Commun. Heat Mass Transf. 2008, 35, 1340–1343. [Google Scholar] [CrossRef]
  42. Pal, D. Combined effects of non-uniform heat source/sink and thermal radiation on heat transfer over an unsteady stretching permeable surface. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1890–1904. [Google Scholar] [CrossRef]
  43. Sajid, M.; Ali, N.; Javed, T.; Abbas, Z. Stretching a Curved Surface in a Viscous Fluid. Chin. Phys. Lett. 2010, 27, 24703. [Google Scholar] [CrossRef]
  44. Abbas, Z.; Naveed, M.; Sajid, M. Heat transfer analysis for stretching flow over a curved surface with magnetic field. J. Eng. Thermophys. 2013, 22, 337–345. [Google Scholar] [CrossRef]
  45. Abbas, Z.; Naveed, M.; Sajid, M. Hydromagnetic slip flow of nanofluid over a curved stretching surface with heat generation and thermal radiation. J. Mol. Liq. 2016, 215, 756–762. [Google Scholar] [CrossRef]
  46. Rosca, N.C.; Pop, I. Unsteady boundary layer flow over a permeable curved stretching/shrinking surface. Eur. J. Mech.—B/Fluids 2015, 51, 61–67. [Google Scholar] [CrossRef]
  47. Naveed, M.; Abbas, Z.; Sajid, M. MHD flow of Micropolar fluid due to a curved stretching sheet with thermal radiation. J. Appl. Fluid Mech. 2016, 9, 131–138. [Google Scholar] [CrossRef]
  48. Hayat, T.; Rashid, M.; Imtiaz, M.; Alsaedi, A. MHD convective flow due to a curved surface with thermal radiation and chemical reaction. J. Mol. Liq. 2017, 225, 482–489. [Google Scholar] [CrossRef]
  49. Khan, U.; Ahmed, N.; Mohyud-Din, S.T.; Sikander, W. Flow of carbon nanotubes suspended nanofluid in stretchable non-parallel walls. Neural Comput. Appl. 2017, 1–13. [Google Scholar] [CrossRef]
  50. Rosseland, S. Astrophysik Und Atom-Theoretische Grundlagen; Springer: Berlin, Germany, 1931. [Google Scholar]
  51. Magyari, E.; Pantokratoras, A. Note on the effect of thermal radiation in the linearized Rosseland approximation on the heat transfer characteristics of various boundary layer flows. Int. Commun. Heat Mass Transf. 2011, 38, 554–556. [Google Scholar] [CrossRef]
  52. Ali, M.E. Heat transfer characteristics of a continuous stretching surface. Heat Mass Transf. 1994, 429, 227–234. [Google Scholar] [CrossRef]
  53. Grubka, L.J.; Bobba, K.M. Heat transfer characteristics of a continuous, stretching surface with variable temperature. ASME J. Heat Transf. 1985, 107, 248–250. [Google Scholar] [CrossRef]
  54. Ishak, A.; Nazar, R.; Pop, I. Boundary layer flow and heat transfer over an unsteady stretching vertical surface. Meccanica 2009, 44, 369–375. [Google Scholar] [CrossRef]
Figure 1. Physical illustration of the flow problem.
Figure 1. Physical illustration of the flow problem.
Applsci 08 00395 g001
Figure 2. (a) f ( ζ ) , f ( ζ ) and Θ ( ζ ) ; and (b) P ( ζ ) for some values of κ .
Figure 2. (a) f ( ζ ) , f ( ζ ) and Θ ( ζ ) ; and (b) P ( ζ ) for some values of κ .
Applsci 08 00395 g002
Figure 3. (a) f ( ζ ) and Θ ( ζ ) ; and (b) P ( ζ ) for some values of Φ .
Figure 3. (a) f ( ζ ) and Θ ( ζ ) ; and (b) P ( ζ ) for some values of Φ .
Applsci 08 00395 g003
Figure 4. Θ ( ζ ) for some values of R d .
Figure 4. Θ ( ζ ) for some values of R d .
Applsci 08 00395 g004
Figure 5. Θ ( ζ ) for some values of λ .
Figure 5. Θ ( ζ ) for some values of λ .
Applsci 08 00395 g005
Figure 6. Skin friction coefficient for some values of: (a) Φ ; and (b) κ .
Figure 6. Skin friction coefficient for some values of: (a) Φ ; and (b) κ .
Applsci 08 00395 g006
Figure 7. Nusselt number for some values of (a) λ ; (b) κ .
Figure 7. Nusselt number for some values of (a) λ ; (b) κ .
Applsci 08 00395 g007
Table 1. Physical and thermal properties of pure water and CNTs [11,27]. CNTs: carbon nanotubes; SWCNTs: Single-walled carbon nanotubes; MWCNTs: multi-walled carbon nanotubes.
Table 1. Physical and thermal properties of pure water and CNTs [11,27]. CNTs: carbon nanotubes; SWCNTs: Single-walled carbon nanotubes; MWCNTs: multi-walled carbon nanotubes.
Pure WaterSWCNTsMWCNTs
C p   ( J ( kgK ) 1 ) 4179 425 796
ρ `   ( kg ( m ) 3 ) 997.1 2600 1600
k   ( W ( mK ) 1 ) 0.613 6600 3000
Table 2. Variation in Thermo-physical properties of CNTs based nanofluid.
Table 2. Variation in Thermo-physical properties of CNTs based nanofluid.
Φ ρ ` ρ ` C p × 10 6 k
SWCNT0.0997.14.16690.613
0.051077.254.01381.1673
0.101157.393.86071.7831
0.151237.543.70762.4712
0.201317.683.55453.2451
MWCNT0.0997.14.16690.613
0.051027.254.02221.1164
0.101057.393.87761.6755
0.151087.543.73292.3002
0.201117.683.58823.0027
Table 3. Comparison of R e s * 1 / 2 C f between present results and previously existing results of Abbas et al. [44] for various values of Φ and κ .
Table 3. Comparison of R e s * 1 / 2 C f between present results and previously existing results of Abbas et al. [44] for various values of Φ and κ .
Φ κ Abbas et al. [44]SWCNTMWCNT
0.051.157631.157631.15763
101.073491.073491.07349
201.035611.035611.03561
501.014051.014051.01405
1001.007041.007041.00704
2001.003561.003561.00356
10001.000791.000791.00079
1.000001.000001.00000
0.15-1.437811.38677
10-1.325771.27251
20-1.275791.22190
0.25-1.808021.69810
10-1.651951.53544
20-1.583231.46471
Table 4. Comparison of R e s * 1 / 2 C f and R e s * 1 / 2 N u s * for different values of Φ and λ by setting R d = 0 and κ = 1000 fixed.
Table 4. Comparison of R e s * 1 / 2 C f and R e s * 1 / 2 N u s * for different values of Φ and λ by setting R d = 0 and κ = 1000 fixed.
λ = 0.3 λ = 0.0 λ = 0.3
Φ R e s * 1 / 2 C f R e s * 1 / 2 N u s * R e s * 1 / 2 C f R e s * 1 / 2 N u s * R e s * 1 / 2 C f R e s * 1 / 2 N u s *
SWCNT0.01.000013.193581.000012.876611.000012.51646
0.11.229065.063711.229064.439951.229063.66480
0.21.519456.423091.519455.472681.519454.01986
MWCNT0.01.000013.193581.000012.876611.000012.51646
0.11.174774.956141.174774.365491.174773.64727
0.21.399446.297341.399445.427431.399444.25578
Table 5. Comparison of Θ ( 0 ) between present results and previously existing results for different values of P r by setting R d = λ = Φ = 0 and κ .
Table 5. Comparison of Θ ( 0 ) between present results and previously existing results for different values of P r by setting R d = λ = Φ = 0 and κ .
Pr Ali [52]Grubka and Bobba [53]Ishak et al. [54]Present Results
0.720.80580.80860.80860.80884
11.00001.00001.00001.00001
31.92371.91441.92371.92368
7--3.07233.07226
103.70063.72073.72073.72068
100-12.294012.294112.29407

Share and Cite

MDPI and ACS Style

Saba, F.; Ahmed, N.; Hussain, S.; Khan, U.; Mohyud-Din, S.T.; Darus, M. Thermal Analysis of Nanofluid Flow over a Curved Stretching Surface Suspended by Carbon Nanotubes with Internal Heat Generation. Appl. Sci. 2018, 8, 395. https://doi.org/10.3390/app8030395

AMA Style

Saba F, Ahmed N, Hussain S, Khan U, Mohyud-Din ST, Darus M. Thermal Analysis of Nanofluid Flow over a Curved Stretching Surface Suspended by Carbon Nanotubes with Internal Heat Generation. Applied Sciences. 2018; 8(3):395. https://doi.org/10.3390/app8030395

Chicago/Turabian Style

Saba, Fitnat, Naveed Ahmed, Saqib Hussain, Umar Khan, Syed Tauseef Mohyud-Din, and Maslina Darus. 2018. "Thermal Analysis of Nanofluid Flow over a Curved Stretching Surface Suspended by Carbon Nanotubes with Internal Heat Generation" Applied Sciences 8, no. 3: 395. https://doi.org/10.3390/app8030395

APA Style

Saba, F., Ahmed, N., Hussain, S., Khan, U., Mohyud-Din, S. T., & Darus, M. (2018). Thermal Analysis of Nanofluid Flow over a Curved Stretching Surface Suspended by Carbon Nanotubes with Internal Heat Generation. Applied Sciences, 8(3), 395. https://doi.org/10.3390/app8030395

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