Next Article in Journal
Monitoring Techniques of Cerium Stabilized Zirconia for Medical Prosthesis
Previous Article in Journal
Lean-Agile Adaptations in Clinical Laboratory Accredited ISO 15189
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels

1
Department of Mathematics, Faculty of Sciences, HITEC University, 47080 Taxila, Pakistan
2
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2015, 5(4), 1639-1664; https://doi.org/10.3390/app5041639
Submission received: 21 September 2015 / Revised: 18 November 2015 / Accepted: 19 November 2015 / Published: 4 December 2015
(This article belongs to the Section Nanotechnology and Applied Nanosciences)

Abstract

:
This article is dedicated to analyzing the heat transfer in the flow of water-based nanofluids in a channel with non-parallel stretchable walls. The magnetohydrodynamic (MHD) nature of the flow is considered. Equations governing the flow are transformed into a system of nonlinear ordinary differential equations. The said system is solved by employing two different techniques, the variational iteration method (VIM) and the Runge-Kutta-Fehlberg method (RKF). The influence of the emerging parameters on the velocity and temperature profiles is highlighted with the help of graphs coupled with comprehensive discussions. A comparison with the already existing solutions is also made, which are the special cases of the current problem. It is observed that the temperature profile decreases with an increase in the nanoparticle volume fraction. Furthermore, a magnetic field can be used to control the possible separation caused by the backflows in the case of diverging channels. The effects of parameters on the skin friction coefficient and Nusselt number are also presented using graphical aid. The nanoparticle volume fraction helps to reduce the temperature of the channel and to enhance the rate of heat transfer at the wall.

1. Introduction

Flow through non-parallel walls is an important area of research due to its many practical, industrial, physical and biological applications, such as flows through rivers and canals. Blood also exhibits such a flow, where the arteries and capillaries are linked with each other. Pioneering work in this area was done by Jeffery and Hamel [1,2]. Plenty of studies are available in this area inspecting the various properties of flow by considering various external forces, such as external magnetic fields and the other forces. One can easily find enough literature on this topic in [3,4] and the references therein.
The study of stretching or shrinking sheets has been an area of modern research due to the practical applications in the real world. These applications include glass fiber production, engineering melt spinning, the manufacturing of rubber sheets, etc. The work presented by Crane [5] was the first study on stretching flows. Since then, several researchers from all over the world have used the said problem to formulate and present the various properties of flows over the stretching/shrinking surfaces that can be seen in [6,7,8] and some of the references therein. They used heat and mass transfer effects and analyzed the various properties of flows over stretching/shrinking surfaces.
Traditional fluids, like water, kerosene, ethylene glycol, etc., in general, are not good conductors. To cope with this problem and to enhance the thermal properties of various fluids, Choi [9,10] presented a model that used the term nanofluids. The idea is to add an appropriate amount of nanoparticles to the traditional fluids. The experimental results have proven that the thermal properties of traditional fluids can be enhanced appreciably by using the said technique. After this groundbreaking idea, many models were presented by various researchers that had also incorporated the Brownian motion and thermophoresis effects. Buongiorno [11], Xue [12], Hamilton and Crosser [13] and Maxwell [14] are the major contributors who have presented some useful and efficient models for nanofluids. In recent times, many authors have been inclined to work in the field of nanofluids. Using Buongiorno’s model, Khan et al. [15] extended the work on nanofluids to a stretching surface. Since then, many studies have been presented using the models discussed earlier for different geometries and various types of nanoparticles. One can find enough literature on these topics [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35].
The complex nature of the equations in daily life problems makes exact solutions least likely. Therefore, many analytical, as well as numerical techniques are developed to solve the problems representing the physical nature. Some of the analytical techniques are Adomian’s decomposition method (ADM), the homotopy perturbation method (HPM), the variational iteration method (VIM), the variation of parameters method (VPM), the homotopy analysis method (HAM), etc. In this article, VIM is successfully applied to the problem, and the efficiency and accuracy of VIM can be seen in [36,37] and the references therein.
In a recent study, Turkyilmazoglu [38] extended the traditional Jeffery-Hamel flow problem to a case incorporating stretching/shrinking walls. He also studied the heat transfer effects on the flow. In the study under consideration, we examined a further extended form of the same problem dealing with the nanofluid flow under the influence of magnetohydrodynamic forces. Due to the very abstract nature of the governing equations and the expected non-availability of the exact solutions, an analytical approximation technique called the variational iteration method (VIM) has been used to approximate the solution. A comparison of analytical results is made with the numerical ones obtained by employing theRunge-Kutta-Fehlberg technique and is presented in the form of tables. Existing solutions (of the limited cases) and the ones obtained in this article are also compared to verify the solutions obtained. A comprehensive graphical description of the effects of the various parameters on the velocity and temperature profiles is presented coupled with some detailed discussions.

2. Governing Equations

Consider the flow from the source or sink at the intersection of two plane walls. 2α is taken as the angle between the walls, as depicted in Figure 1. The walls are assumed to be stretching/shrinking at a rate s   in such a way that the velocity at the wall can be written as u r = U w = s r . Furthermore, α > 0 represents the case of a divergent channel, while for α < 0 , the same channel is convergent. The channel is assumed to be saturated by the nanofluid consisting of water (the base fluid) and the copper (Cu) or silver (Ag) nanoparticles. A thermal equilibrium between the base fluid and the nanoparticles is assumed. A uniform transverse magnetic field of strength B 0 is applied across the channel to analyze the flow under the influence of an external body force.
Figure 1. Schematic diagram of the flow problem.
Figure 1. Schematic diagram of the flow problem.
Applsci 05 01639 g001
For the said flow, a polar coordinates system ( r , θ ) is used to formulate the problem. Under the said conventions, the velocity field takes the form V = ( u r , 0, 0), where u r is a function of both r and θ .
The equations of continuity, momentum and energy in polar coordinates under the imposed assumptions become [38]:
1 r r ( r u r ) = 0
ρ n f ( u r u r r ) = p r + μ n f [ 2 u r r 2 + 1 r u r r + 1 r 2 2 u r θ 2 u r r 2 ] σ n f B 0 2 ρ n f r 2 u r
1 ρ n f r p r + 2 r 2 μ n f ρ n f u r θ = 0
u r T r = k n f ( ρ C p ) n f [ 2 T r 2 + 1 r T r + 1 r 2 2 T θ 2 ] + μ n f ( ρ C p ) n f [ 4 ( u r r ) 2 + 1 r 2 ( u r θ ) 2 ]
Boundary conditions for the problem are:
u r = U c r , u r θ = 0 , T θ = 0   at θ = 0 u r = U w = s r , T = T w r 2   at θ = ± α
where U c is the velocity at the centerline of the channel and U w , T w , respectively, represent the velocity and temperature at the wall of the channel.
Further, ρ n f , μ n f , σ n f , ( ρ C p ) n f and k n f denote the density, dynamic viscosity, electric conductivity, heat capacity and the thermal conductivity of the fluid, respectively [27,28,29]; where,
υ n f = μ n f ρ n f , μ n f = μ f ( 1 φ ) 2.5 , α n f = k n f ( ρ C p ) n f , ρ n f = ( 1 φ ) + φ ρ s ρ f , ( ρ C p ) n f = ( 1 φ ) + φ ( ρ C p ) s ( ρ C p ) f
The Maxwell-Garnett model for the effective thermal conductivity of nanofluids is used in this problem as [14]:
k n f k f = k s + 2 k f 2 φ ( k f k s ) k s + 2 k f + φ ( k f k s )
As proposed by Maxwell [18], the effective electrical conductivity of the nanofluid can be written as:
σ n f σ f = 1 + 3 ( σ s σ f 1 ) φ ( σ s σ f + 2 ) ( σ s σ f 1 ) φ
Here, μ f is the viscosity of the base fluid and φ is the volume fraction for the nanoparticles. k f , k s are the thermal conductivities, and ρ f , ρ s are the densities of the base fluid and nanoparticles, respectively. Furthermore, σ f and σ s , in order, represent the electrical conductivity of the base fluid and nanoparticles.
The continuity equation (Equation (1)) gives a clear description of the radial velocity of the form,
f ( θ ) = r u r ( r , θ )
Using the following non-dimensional variables,
F ( η ) = f ( θ ) U c ,   η = θ α , Θ = r 2 T T w
Elimination of the pressure terms from Equations (2) and (3) and the implementation of similarity transform Equation (9) gives us the coupled equations for the velocity and temperature profiles as:
F ( η ) + 2 α Re ( 1 φ ) 2.5 K 1 F ( η ) F ( η ) + ( 4 H a ( 1 φ ) 2.5 A 1 ) α 2 F ( η ) = 0
Θ ( η ) + 4 α 2 Θ ( η ) + 2 K 2 P r K 3 α 2 Θ ( η ) F ( η ) + K 2 E c P r K 3 R e ( 1 ϕ ) 2.5 [ 4 α 2 F 2 ( η ) + ( F ( η ) ) 2 ] = 0
Use of Equation (9) reduces the boundary conditions for the problem to a more compact form as:
F ( 0 ) = 1 , F ( 0 ) = 0 , F ( 1 ) = S    Θ ( 0 ) = 0 , Θ ( 1 ) = 1
In the above equations, S denotes the stretching/shrinking parameter and S > 0 stands for stretchable walls, while S < 0 represents the case of shrinking walls. Re here is the Reynolds number given by:
Re = U c α υ f ( D i v e r g e n t C h a n n e l :   α > 0 , U > 0 C o n v e r g e n t C h a n n e l : α < 0 , U < 0 )
Further, H a = σ f B 0 2 μ f , E c = U c 2 α ( C p ) f T w , Pr = ( ρ C p ) f U c k η f represent the Hartmann number, the Eckert number and the Prandtl number, respectively; also:
K 1 = ( ( 1 φ ) + φ ρ s ρ f ) , K 2 = ( 1 φ ) + φ ( ρ C p ) s ( ρ C p ) f , K 3 = k n f k f = k s + 2 k f 2 φ ( k f k s ) k s + 2 k f + φ ( k f k s ) , A 1 = σ n f σ f = 1 + 3 ( σ s σ f 1 ) φ ( σ s σ f + 2 ) ( σ s σ f 1 ) φ
Quantities of physical interest are the skin friction coefficient and the Nusselt number, defined as:
C F = μ n f ( τ r θ ) η = 1 ρ n f U 2 , N u = l k f ( q w ) η = 1 k T w
In terms of Equations (8)–(11), we have:
Re r C f = 1 K 1 ( 1 φ ) 2.5 F ( 1 )
α N u = K 3 Θ ( 1 )

3. Solution Procedure

The solution of the system of equations governing the flow is obtained by using VIM. Implementing the traditional form of VIM, the velocity and temperature profiles in recursive form can be written in the form,
F n + 1 ( η ) = F n ( η ) + 0 η λ F ( s ) ( F ( s ) + ( 4 ( 1 φ ) 2.5 A 1 H a ) α 2 F ˜ ( s ) + ( 1 φ ) 2.5 A 1 2 α Re F ˜ ( s ) F ˜ ( s ) ) d s
Θ n + 1 ( η ) = Θ n ( η ) 0 η λ Θ [ [ Θ ( s ) + 4 α 2 Θ ( s ) + 2 K 2 P r K 3 α 2 Θ ( s ) F ( s ) K 2 E c P r K 3 R e ( 1 ϕ ) 2.5 [ 4 α 2 F 2 ( s ) + ( F ( s ) ) 2 ] ] d s
where the Lagrange multipliers for velocity and temperature profile are λ F ( s ) and λ Θ ( s ) , respectively. Approximate Lagrange multipliers for the above problem can be obtained using the variational approach as λ F ( s ) = ( s η ) 2 2 ! and λ Θ ( s ) = ( s η ) , so that the above iterative scheme can be reformulated as:
F n + 1 ( η ) = F n ( η ) 0 η ( s η ) 2 2 ! ( F ( s ) + ( 4 ( 1 φ ) 2.5 A 1 H a ) α 2 F ( s ) + ( 1 φ ) 2.5 K 1 2 α Re F ( s ) F ( s ) ) d s
Θ n + 1 ( η ) = Θ n ( η ) + 0 η ( s η ) [ [ Θ ( η ) + 4 α 2 Θ ( η ) + 2 K 2 P r K 3 α 2 Θ ( η ) F ( η ) K 2 E c P r K 3 R e ( 1 ϕ ) 2.5 [ 4 α 2 F 2 ( η ) + ( F ( η ) ) 2 ] ] d s + 0 η ( s η ) [ [ Θ ( η ) + 4 α 2 Θ ( η ) + 2 K 2 P r K 3 α 2 Θ ( η ) F ( η ) K 2 E c P r K 3 R e ( 1 ϕ ) 2.5 [ 4 α 2 F 2 ( η ) + ( F ( η ) ) 2 ] ] d s
For the initial guess, using boundary conditions for velocity and temperature profiles, we get:
F 0 ( η ) = 1 + β η 2 2 ,
Θ 0 ( η ) = γ
where β and γ are the constants that will be calculated by using the boundary conditions F ( 1 ) = S and Θ ( 1 ) = 1 .
Next, a few terms of the solution are given as:
F 1 ( η ) = 1 + β η 2 2 0.008333333333333318 R e α β 2 η 6 ( 1. 1. ϕ ) 2.5 K 1 + η 4 ( 0.16666666666666674 α 2 + 0.041666666666666685 Ha α 2 β ( 1. 1. ϕ ) 2.5 A 1 0.08333333333333337 Re α β ( 1. 1. ϕ ) 2.5 K 1 )
Θ 0 ( η ) = γ + η 2 ( 2. α 2 γ 1. Pr α 2 γ K 2 K 3 2. EcPr α 2 K 2 ( 1. 1. ϕ ) 2.5 K 3 ) + η 4 ( 0.08333333333333331 Pr α 2 β γ K 2 K 3 0.33333333333333326 EcPr α 2 β K 2 ( 1. 1. ϕ ) 2.5 K 3 0.08333333333333331 EcPr β 2 K 2 ( 1. 1. ϕ ) 2.5 K 3 ) 0.033333333333333354 EcPr α 2 β 2 η 6 K 2 ( 1. 1. ϕ ) 2.5 K 3
For the sake of simplicity, other iterations of the solutions are omitted here. One can calculate the next iterations in a similar manner.

4. Results and Discussion

To study the flow behavior and the variations in velocity and temperature due to the varying values of the involved parameters, this section is divided into two subsections. One is for the velocity profile, and the other is for the temperature distribution. Both the stretching and shrinking cases combined with the convergence or the divergence of the channel are discussed in detail using the graphical aid. As water is taken to be the base fluid, the numerical values of the Prandtl number Pr shall be fixed at 6.2 throughout the study, unless stated otherwise. Copper (Cu) and silver (Ag) are the nanoparticles incorporated in the base fluid. The thermophysical properties of the base fluid and the nanoparticles are highlighted in Table 1. Furthermore, in all of the graphical descriptions, a solid line represents the graph for Cu, while the dashed line represents the same for Ag.
Table 1. Thermophysical properties of water, copper and silver nanoparticles [22].
Table 1. Thermophysical properties of water, copper and silver nanoparticles [22].
ρ ( k g / m 3   ) C p ( J / K g K ) k ( W / m K ) Electrical Conductivity σ ( Ω m ) 1
Pure Water997.141790.6130.05
Copper (Cu)89333854015.96 × 107
Silver10,5002354296.3 × 107

4.1. Velocity Profile

This section is dedicated to exploring the behavior of the velocity profile for both convergent and divergent channels under the variations in the parameters, like the channel opening   α , the Reynolds number   R e , the volume fraction of nanoparticles ϕ, the Hartmann number   H a and the stretching/shrinking parameter S. For the said purpose, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 are plotted. Figure 2 gives a graphical description of the variations in the velocity profile for a divergent channel for different values of channel openin   α . Both the stretching and shrinking cases are discussed in the same figure. A clear drop in the velocity profile is evident for increasing values of α. Here, clearly, for the stretching walls, the change in velocity is seen to be quite rapid at the central portion; there comes a point at which the backflow phenomena are observed particularly for α = 5 ° . This backflow behavior may result in separation at some point. In Figure 3, the same effects of α for the convergent channel are portrayed. A quite opposite behavior is observed for this case, and increasing values of velocity are observed. For the stretching walls, the change in velocity is on the lower side as compared to the shrinking walls. This means that for a shrinking, convergent channel, more fluid is moved towards the walls as compared to the same for a diverging channel. Copper nanoparticles are seen to have a slightly higher velocity for divergent channels as compared to the silver nanoparticles. The behavior is reversed for a convergent channel.
Figure 2. Change in velocity for varying α (divergent channel).
Figure 2. Change in velocity for varying α (divergent channel).
Applsci 05 01639 g002
Figure 3. Change in velocity for varying α (convergent channel).
Figure 3. Change in velocity for varying α (convergent channel).
Applsci 05 01639 g003
Figure 4. Change in velocity for varying Re (divergent channel).
Figure 4. Change in velocity for varying Re (divergent channel).
Applsci 05 01639 g004
Figure 5. Change in velocity for varying Re (convergent channel).
Figure 5. Change in velocity for varying Re (convergent channel).
Applsci 05 01639 g005
Figure 6. Change in velocity for varying H a (divergent channel).
Figure 6. Change in velocity for varying H a (divergent channel).
Applsci 05 01639 g006
Figure 7. Change in velocity for varying H a (convergent channel).
Figure 7. Change in velocity for varying H a (convergent channel).
Applsci 05 01639 g007
Figure 8. Change in velocity for varying φ (divergent channel).
Figure 8. Change in velocity for varying φ (divergent channel).
Applsci 05 01639 g008
Figure 9. Change in velocity for varying φ (convergent channel).
Figure 9. Change in velocity for varying φ (convergent channel).
Applsci 05 01639 g009
Figure 10. Change in velocity for varying S (divergent channel).
Figure 10. Change in velocity for varying S (divergent channel).
Applsci 05 01639 g010
Figure 11. Change in velocity for varying S (convergent channel).
Figure 11. Change in velocity for varying S (convergent channel).
Applsci 05 01639 g011
The velocity profile for the increasing values of the Reynolds number Re for the case of a divergent channel is plotted in Figure 4. The higher the values of Re, the lower the velocity becomes in the diverging channel. This results in a backflow near the walls of the channel for stretching walls. A steep decrease in velocity is seen for the stretching/diverging channel as compared to the shrinking one. A quite opposite scenario is seen for the convergent channel, and velocity is an increasing function of Re. This entire discussion shows that the upsurge in viscous forces results in lower values of the Reynolds number, and an increase in the mainstream velocity of the channel increases the values of Re that result in the increasing of the velocity for convergent channels and a decrease in velocity for the divergent channel. Again, for increasing Re, the copper nanoparticles have slightly higher values for the velocity of the divergent channel. This is due to the lower density of copper as compared to silver.
The Hartmann number that represents the strength of the magnetic field results in a growth in the velocity profile when the channel is divergent, as is clear from Figure 6. The behavior is similar for both the stretching and shrinking of walls. Furthermore, whenever there is a backflow in the divergent channel, the Hartmann number can be used to reduce the backflows. It controls the separation phenomena, resulting in a much smoother flow. The velocity behaves in a similar manner for a convergent channel and the increasing values of Ha. However, in this case, the change is slower as compared to the divergent channel. Furthermore, a stretching/converging channel has lower variations in velocity as compared to the shrinking channel (see Figure 7).
To discuss the strength of nanoparticles on the velocity, Figure 8 and Figure 9 are portrayed. A rise in the nanoparticle volume fraction reduces the velocity in a divergent channel. For the stretching of a divergent channel, again, a backflow is observed. However, when the channel is shrinking, there is no backflow. Again, a slightly higher value of the velocity can be observed for the copper nanoparticles than the silver nanoparticles. From all of this, one can easily say for a divergent channel that the higher volume fraction of the nanoparticles provides a low speed flow. On the other hand, a slow rise in the velocity for a convergent channel is evident: it keeps growing with the rising nanoparticle volume fraction.
The influence of stretching and shrinking parameter S on the velocity profile is depicted in Figure 10 and Figure 11. For a divergent channel, an upsurge in the velocity is witnessed. The stretching of a divergent channel marks a higher velocity near the walls as compared to the center of the channel. On the contrary, when a divergent channel is shrinking, the velocity decreases quite sharply near the walls, as shown in Figure 10. In Figure 11, the effects of the same parameter on a convergent channel are portrayed. It gives a clear picture of an accelerated flow near the walls for a stretching channel. On the contrary, a significant decrease in velocity is seen in a convergent channel with shrinking walls. For both of these cases, the changes are quite prominent near the walls of the channel. As we move towards the center, the effects get weaker and less influential.

4.2. Temperature Profile

Variations in the temperature profile under the influence of various parameters, like the angle of opening   α , the volume fraction of nanoparticle   ϕ , the stretching/shrinking parameter   S and the Eckert numbe   E c , are plotted in Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23, Figure 24, Figure 25, Figure 26 and Figure 27. Here, the value of the Prandtl number   P r is fixed at 6.2 because we are taking water as a base fluid. The temperature profile under the effect of angle opening α is plotted in Figure 12, Figure 13, Figure 14 and Figure 15. Almost an identical behavior to the temperature profile is observed for both the divergent and convergent channels due to the increasing values of   α . When the walls are stretching, the temperature is seen to be on a slightly lower side for rising   α . On the other hand, for shrinking walls, some higher values of temperature are seen when there is no opening of the walls   ( α = 0 ) . This means whether the walls are stretching or shrinking, with the change in angle, the temperature of the system will rise. There is not much of a difference between the temperature of Cu and Ag.
Figure 12. Change in temperature for varying α (divergent channel/stretching case).
Figure 12. Change in temperature for varying α (divergent channel/stretching case).
Applsci 05 01639 g012
Figure 13. Change in temperature for varying α (divergent channel/shrining case).
Figure 13. Change in temperature for varying α (divergent channel/shrining case).
Applsci 05 01639 g013
Figure 14. Change in temperature for varying α (convergent channel/stretching case).
Figure 14. Change in temperature for varying α (convergent channel/stretching case).
Applsci 05 01639 g014
Figure 15. Change in temperature for varying α (convergent channel/shrinking case).
Figure 15. Change in temperature for varying α (convergent channel/shrinking case).
Applsci 05 01639 g015
Figure 16. Change in temperature for varying φ (divergent channel/stretching case).
Figure 16. Change in temperature for varying φ (divergent channel/stretching case).
Applsci 05 01639 g016
Figure 17. Change in temperature for varying φ (divergent channel/shrinking case).
Figure 17. Change in temperature for varying φ (divergent channel/shrinking case).
Applsci 05 01639 g017
Figure 18. Change in temperature for varying φ (convergent channel/stretching case).
Figure 18. Change in temperature for varying φ (convergent channel/stretching case).
Applsci 05 01639 g018
Figure 19. Change in temperature for varying φ (convergent channel/shrinking case).
Figure 19. Change in temperature for varying φ (convergent channel/shrinking case).
Applsci 05 01639 g019
Figure 20. Change in temperature for varying S (divergent channel/stretching case).
Figure 20. Change in temperature for varying S (divergent channel/stretching case).
Applsci 05 01639 g020
Figure 21. Change in temperature for varying S (divergent channel/shrinking case).
Figure 21. Change in temperature for varying S (divergent channel/shrinking case).
Applsci 05 01639 g021
Figure 22. Change in temperature for varying S (convergent channel/stretching case).
Figure 22. Change in temperature for varying S (convergent channel/stretching case).
Applsci 05 01639 g022
Figure 23. Change in temperature for varying S (convergent channel/shrinking case).
Figure 23. Change in temperature for varying S (convergent channel/shrinking case).
Applsci 05 01639 g023
Figure 24. Change in temperature for varying E c (divergent channel/stretching case).
Figure 24. Change in temperature for varying E c (divergent channel/stretching case).
Applsci 05 01639 g024
Figure 25. Change in temperature for varying E c (divergent channel/shrinking case).
Figure 25. Change in temperature for varying E c (divergent channel/shrinking case).
Applsci 05 01639 g025
Figure 26. Change in temperature for varying E c (convergent channel/stretching case).
Figure 26. Change in temperature for varying E c (convergent channel/stretching case).
Applsci 05 01639 g026
Figure 27. Change in temperature for varying E c (convergent channel/shrinking case).
Figure 27. Change in temperature for varying E c (convergent channel/shrinking case).
Applsci 05 01639 g027
The influence of the nanoparticle volume fraction on the temperature profile is plotted in Figure 16, Figure 17, Figure 18 and Figure 19. An addition of the nanoparticles in the base fluid decreases the temperature of the system quite considerably. This change is prominent at the centerline of the channel. For a divergent channel, the temperature profile is somewhat on the higher side when the walls are stretching as compared to the shrinking of the walls. For a convergent channel, this behavior is reversed. For all of these cases, silver (Ag) nanoparticles tend to have slightly lower temperature values than copper (Cu) nanoparticles. The differences in the specific heat and the thermal conductivity of Cu and Ag are the major factors for the change in temperature.
The next set of figures describes the behavior of temperature with the stretching/shrinking of the walls. Both the convergent and divergent channel cases are discussed. Figure 20 and Figure 21 depict that the stretching and shrinking of the walls have a quite opposite effect on the temperature profile for a divergent channel. When the walls are stretching, a steep drop in the temperature is observed at the center of the channel. The change becomes smaller as we move near the walls of the channel. On the other hand, the shrinking walls tend to increase the temperature quite substantially at the central portion of the channel. For a convergent channel, the temperature behaves almost in an identical manner as the one seen in a divergent channel; Figure 22 and Figure 23 are here to back our remarks. The only difference is the slightly higher values of the temperature for the stretching walls in a divergent channel. Interestingly, Cu nanoparticles have marginally lower values for the stretching of a divergent channel as compared to the stretching of a convergent channel. For a convergent channel, the temperature of Cu seems to be on the higher side. However, when the walls are shrinking, Cu nanoparticles show a marginal advantage in temperature.
Figure 24, Figure 25, Figure 26 and Figure 27 are plotted to study the changes occurring in the temperature profile due to the variations in   E c . Here, Ec represents the influence of dissipation on the temperature profile. As expected, the Eckert number increases the temperature quite expressively. This change in the temperature is slightly lower for the case when the walls are stretching. For the shrinking walls, the temperature rises quite efficiently with increasing values of the Eckert number. Furthermore, for a convergent channel and the stretching walls, the temperature curves deviate relatively less.

4.3. Skin Friction Coefficient and the Nusselt Number

The effects of various parameters on the skin friction coefficient and the Nusselt number are plotted in Figure 28, Figure 29, Figure 30, Figure 31, Figure 32, Figure 33, Figure 34 and Figure 35. Both stretching and shrinking cases are taken into account. For the case of a divergent channel, Figure 28 gives a description of the varying angle and the nanoparticle volume fraction on the skin friction coefficient. Clearly, with an increase in the nanoparticle volume fraction, the skin friction is seen to be declining for both stretching and shrinking divergent channels. Furthermore, the greater the channel opening for the divergent channel, the lower are the skin friction values. It can be observed that the stretching channel has lower values of the skin friction as compared to the shrinking case. Due to the higher density values, silver nanoparticles have lower values as compared to copper nanoparticles. For the case of the convergent channel, the behavior of the skin friction coefficient can be observed from Figure 29. The skin friction is seen to be increasing with the increase in the angle for the convergent channel. Here, the stretching convergent channel has lower values of the skin friction as compared to the shrinking convergent channel. The manners in which skin friction behaves for increasing values of Re are plotted in Figure 30 and Figure 31. For the divergent channel, the increase in Re results in declining values of the skin friction coefficient, as well as for the nanoparticle volume fraction. Stretching divergent channels have lower skin friction values than shrinking divergent channels. Almost the reverse behavior is observed for the convergent channel, i.e., higher Re values result in higher skin friction values for both stretching and shrinking convergent channels.
Figure 28. Variation in the skin friction coefficient with the angle and the nanoparticle volume fraction (divergent channel).
Figure 28. Variation in the skin friction coefficient with the angle and the nanoparticle volume fraction (divergent channel).
Applsci 05 01639 g028
Figure 29. Variation in the skin friction coefficient with the angle and the nanoparticle volume fraction (convergent channel).
Figure 29. Variation in the skin friction coefficient with the angle and the nanoparticle volume fraction (convergent channel).
Applsci 05 01639 g029
Figure 30. Variation in the skin friction coefficient with the Reynold number and the nanoparticle volume fraction (divergent channel).
Figure 30. Variation in the skin friction coefficient with the Reynold number and the nanoparticle volume fraction (divergent channel).
Applsci 05 01639 g030
Figure 31. Variation in the skin friction coefficient with the Reynold number and the nanoparticle volume fraction (convergent channel).
Figure 31. Variation in the skin friction coefficient with the Reynold number and the nanoparticle volume fraction (convergent channel).
Applsci 05 01639 g031
Figure 32. Variation in the Nusselt number with the angle and the nanoparticle volume fraction (divergent channel).
Figure 32. Variation in the Nusselt number with the angle and the nanoparticle volume fraction (divergent channel).
Applsci 05 01639 g032
Figure 33. Variation in the Nusselt number with the angle and the nanoparticle volume fraction (convergent channel).
Figure 33. Variation in the Nusselt number with the angle and the nanoparticle volume fraction (convergent channel).
Applsci 05 01639 g033
Figure 34. Variation in the Nusselt number with the Eckert number and the nanoparticle volume fraction (divergent channel).
Figure 34. Variation in the Nusselt number with the Eckert number and the nanoparticle volume fraction (divergent channel).
Applsci 05 01639 g034
Figure 35. Variation in the Nusselt number with the Eckert number and the nanoparticle volume fraction (convergent channel).
Figure 35. Variation in the Nusselt number with the Eckert number and the nanoparticle volume fraction (convergent channel).
Applsci 05 01639 g035
Changes in the Nusselt number with varying parameters are plotted in Figure 32, Figure 33, Figure 34 and Figure 35 both for divergent and convergent channels. Figure 32 gives a description of the change in the Nusselt number with the varying nanoparticle volume fraction and the angle for the divergent channel case. The heat transfer rate at the wall is seen to be an increasing function of both the angle and the nanoparticle volume fraction. This behavior is quite fast for a shrinking divergent channel. A similar pattern is seen for the convergent channel, the difference being only that, for the convergent channel case, the values are slightly on the higher side as compared to the same values for the divergent channel, as depicted in Figure 33. Apparently, there is no difference in the values for copper and silver nanoparticles. The increase in the Eckert number results in a higher heat transfer rate at the wall for both divergent and convergent channels. The stretching of walls has a dominant effect as compared to the shrinking of the walls. The nanoparticle volume fraction increases the heat transfer rate and, thus, can be useful in many industrial processes. The problem under consideration is solved using two methods: the variational iteration method and the Runge-Kutta-Fehlberg method. The results obtained are compared to existing solutions when there is no magnetic field and the fluid is taken as water only. An excellent agreement between the solutions is seen, as given in Table 2 and Table 3.
Table 2. Comparison of the numerical values for the skin friction coefficient with the already existing solutions in the literature for Re = 50 ,   H a = 0 ,   φ = 0. VIM, variational iteration method; RKF, Runge–Kutta–Fehlberg method.
Table 2. Comparison of the numerical values for the skin friction coefficient with the already existing solutions in the literature for Re = 50 ,   H a = 0 ,   φ = 0. VIM, variational iteration method; RKF, Runge–Kutta–Fehlberg method.
S F ( 1 ) ( α = 5 ° ) S F ( 1 ) ( α = 5 ° )
VIMRKFTurkyilmazoglu [38]VIMRKFTurkyilmazoglu [38]
−1−3.508103−3.508103−3.508103−2−5.130922−5.130922−5.130922
−1/2−2.173044−2.173044−2.173044−1−4.652159−4.652159−4.652159
00.0000000.0000000.0000000−2.833951−2.833951−2.833951
1/2−0.361846−0.361846−0.36184610.0000000.0000000.000000
10.0000000.0000000.00000023.6697113.6697113.669711
Table 3. Comparison of the numerical values for the Nusselt number with the already existing solutions in the literature for Re = 50 , H a = 0 , φ = 0.
Table 3. Comparison of the numerical values for the Nusselt number with the already existing solutions in the literature for Re = 50 , H a = 0 , φ = 0.
S Θ ( 1 ) ( α = 5 ° ) S Θ ( 1 ) ( α = 5 ° )
VIMRKFTurkyilmazoglu [38]VIMRKFTurkyilmazoglu [38]
−10.0347750.0347750.034775−20.0315760.0315760.031576
−1/20.03726850.03726850.0372685−10.0373220.0373220.037322
00.0399820.0399820.03998200.0421510.0421510.042151
1/20.0429860.0429860.04298610.0464010.0464010.046401
10.0464010.0464010.04640120.0502420.0502420.050242

5. Conclusions

This article is presented to study the flow and heat transfer of water-based nanofluids in stretchable convergent/divergent channels. The magnetohydrodynamic nature of the flow is incorporated. Copper and silver nanoparticles are added to water to analyze the behavior of the velocity and temperature profiles with varying the parameters involved. Equations governing the flow are solved using two methods, namely VIM and RK-F. A comparison with existing solutions is also given to check the accuracy of our obtained results. The effects of the parameters on the skin friction coefficient and heat transfer at the wall are also plotted. The major outcomes of this study are as follows:
  • Angle opening and the Reynolds number have opposite effects on the velocity profile for convergent and divergent channels.
  • Stretching walls may result in backflow regimes for divergent channels, because this moves the particles near the wall away. An increase in the angle opening and Reynolds number may result in backflow and, thus, separation. This result might cause instabilities in the flow.
  • The Hartmann number gives a solution to the backflow regions. An increase in the Hartmann number removes the backflow for stretching divergent channels. This effect can be useful in several physical phenomena.
  • The nanoparticle volume fraction reduces the velocity of the fluid for both copper and silver nanoparticles in the case of a divergent channel. For a convergent channel, the increase in the volume fraction increases the velocity.
  • Stretching of the divergent channel increases the flow near the walls of the channel. Additionally, shrinking reduces the velocity of the fluid near the walls of the channel. Identical behavior is seen for the case of a convergent channel.
  • The almost identical behavior of the temperature profile for increasing channel opening is observed for both convergent and divergent channels, both for stretching and shrinking channels.
  • Temperature is seen to be decreasing for increasing values of the nanoparticle volume fraction. Both stretching and shrinking channels have an almost identical behavior.
  • Stretching of the walls results in lower temperature values for both convergent and divergent channels. Increasing values of temperature near the walls of the channel are observed for the shrinking case.
  • The Eckert number increases the temperature of the fluid for all of the cases.
  • The nanoparticle volume fraction reduces the skin friction coefficient for all of the cases. Stretching channels have lower skin friction values than shrinking channels.
  • The opposite behavior of the skin friction coefficient for increasing Re is observed for convergent and divergent channels.
  • The rate of heat transfer at the wall increases with a higher nanoparticle volume fraction and Eckert number. These two parameters have the same behavior for all of the cases.

Acknowledgments

The authors are highly grateful to the unknown referees for their valuable suggestions, which really helped to improve the standards of the manuscript. This project was supported by King Saud University, Deanship of Scientific Research, College of Sciences Research Center.

Author Contributions

Author S.T.M.-D. developed the problem and its MATHEMATICA code. U.K. in collaboration with first author, did the literature review, developed and implemented the computer code, and interpreted the subsequently obtained results. N.A. did the editing and removed the grammatical mistakes. S.M.H. 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 that they have no conflict of interest.

References

  1. Jeffery, G.B. The two-dimensional steady motion of a viscous fluid. Philos. Mag. 1915, 6, 455–465. [Google Scholar] [CrossRef]
  2. Hamel, G. Spiralförmige Bewgungen ZäherFlüssigkeiten. Jahresber. Deutsch. Math. Verein. 1916, 25, 34–60. [Google Scholar]
  3. Asadullah, M.; Khan, U.; Ahmed, N.; Manzoor, R.; Mohyud-Din, S.T. MHD Flow of a Jeffery Fluid in Converging and Diverging Channels. Int. J. Mod. Math. Sci. 2013, 6, 92–106. [Google Scholar]
  4. Motsa, S.S.; Sibanda, P.; Marewo, G.T. On a new analytical method for flow between two inclined walls. Numer. Algorithms 2012, 61, 499–514. [Google Scholar] [CrossRef]
  5. Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. 1970, 21, 645–647. [Google Scholar] [CrossRef]
  6. Akbar, N.S.; Nadeem, S.; Haq, R.U.; Khan, Z.H. Radiation effects on MHD stagnation point flow of nano fluid towards a stretching surface with convective boundary condition. Chin. J. Aeronaut. 2013, 26, 1389–1397. [Google Scholar] [CrossRef]
  7. Nadeem, S.; Haq, R.U. Effect of Thermal Radiation for Megnetohydrodynamic Boundary Layer Flow of a Nanofluid Past a Stretching Sheet with Convective Boundary Conditions. J. Comput. Theor. Nanosci. 2013, 11, 32–40. [Google Scholar] [CrossRef]
  8. Magyari, E.; Keller, B. Exact solutions for self-similarboundary-layerflows induced by permeable stretching walls. Eur. J. Mech. B Fluids 2000, 19, 109–122. [Google Scholar] [CrossRef]
  9. Choi, S.U.S.; Eastman, J.A. Enhancing thermal conductivity of fluids with nanoparticle. In Proceedings of ASME International Mechanical Engineering Congress & Exposition, San Francisco, CA, USA, 12–17 November 1995.
  10. Choi, S.U.S.; Zhang, Z.G.; Yu, W.; Lockwood, F.E.; Grulke, E.A. Anomalously thermal conductivity enhancement in nanotube suspensions. Appl. Phys. Lett. 2001, 79, 2252–2254. [Google Scholar] [CrossRef]
  11. Buongiorno, J. Convective transport in nanofluids. ASME J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  12. Xue, Q. Model for thermal conductivity of carbon nanotube-based composites. Phys. B Condens. Matter 2005, 368, 302–307. [Google Scholar] [CrossRef]
  13. Hamilton, R.L.; Crosser, O.K. Thermal conductivity of heterogeneous two component systems. Ind. Eng. Chem. Fundam. 1962, 1, 187–191. [Google Scholar] [CrossRef]
  14. Maxwell, J.C. A Treatise on Electricity and Magnetism, 3rd ed.; Clarendon: Oxford, UK, 1904. [Google Scholar]
  15. 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]
  16. Nadeem, S.; Haq, R.U.; Akbar, N.S. MHD three-dimensional boundary layer flow of casson nanofluid past a linearly stretching sheet with convective boundary condition. Nanotechnol. IEEE Trans. 2014, 13, 109–115. [Google Scholar] [CrossRef]
  17. Noor, N.F.M.; Haq, R.U.; Nadeem, S. Mixed convection stagnation flow of a micropolarnanofluid along a vertically stretching surface with slip effects. Meccanica 2015. [Google Scholar] [CrossRef]
  18. Khan, U.; Ahmed, N.; Khan, S.I.U.; Mohyud-Din, S.T. Thermo-diffusion and MHD effects on stagnation point flow towards a stretching sheet in a nanofluid. Propuls. Power Res. 2014, 3, 151–158. [Google Scholar] [CrossRef]
  19. 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. 2015. [Google Scholar] [CrossRef]
  20. 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]
  21. Khalid, A.; Khan, I.; Shafie, S. Exact solutions for free convection flow of nanofluids with ramped wall temperature. Eur. Phys. J. Plus 2015, 130. [Google Scholar] [CrossRef]
  22. Sheikholeslami, M.; Hatami, M.; Ganji, D.D. Nanofluid flow and heat transfer in a rotating system in the presence of a magnetic field. J. Mol. Liq. 2014, 190, 112–120. [Google Scholar] [CrossRef]
  23. Sheikholislami, M.; Ellahi, R. Electrohydrodynamic nanofluid hydrothermal treatment in an enclosure with sinusoidal upper wall. Appl. Sci. 2015, 5, 294–306. [Google Scholar] [CrossRef]
  24. 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]
  25. Ellahi, R.; Hassan, M.; Zeeshan, A. Study on magnetohydrodynamic nanofluid by means of single and multi-walled carbon nanotubes suspended in a salt water solution. IEEE Trans. Nanotechnol. 2015, 14, 726–734. [Google Scholar] [CrossRef]
  26. Sheikholeslami, M.; Ellahi, R. Simulation of ferrofluid flow for magnetic drug targeting using Lattice Boltzmann method. J. Z. Naturforsch. A 2015, 70, 115–124. [Google Scholar]
  27. 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]
  28. Ellahi, R.; Hassan, M.; Zeeshan, A. Shape effects of nanosize particles in Cu-H2O nanofluid on entropy generation. Int. J. Heat Mass Transf. 2015, 81, 449–456. [Google Scholar] [CrossRef]
  29. Sheikholeslami, M.; Ganji, D.D.; Javed, M.Y.; Ellahi, R. Effect of thermal radiation on nanofluid flow and heat transfer using two phase model. J. Magn. Magn. Mater. 2015, 374, 36–43. [Google Scholar] [CrossRef]
  30. Rashidi, S.; Dehghan, M.; Ellahi, R.; Riaz, M.; Jamal-Abad, M.T. Study of stream wise transverse magnetic fluid flow with heat transfer around a porous obstacle. J. Magn. Magn. Mater. 2015, 378, 128–137. [Google Scholar] [CrossRef]
  31. 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]
  32. Sheikholeslami, M.; Hatami, M.; Ganji, D.D. Nanofluid flow and heat transfer in a rotating system in the presence of a magnetic field. J. Mol. Liq. 2014, 190, 112–120. [Google Scholar] [CrossRef]
  33. Ellahi, R.; Bhatti, M.M.; Khan, A.A. Influence of wall compliances on unsteady peristaltic transport of a nanofluid in a vertical rectangular duct. Wulfenia 2015, 22, 248–265. [Google Scholar]
  34. Nawaz, M.; Zeeshan, A.; Ellahi, R.; Abbasbandy, S.; Rashidi, S. Joules heating effects on stagnation point flow over a stretching cylinder by means of genetic algorithm and Nelder-Mead method. Int. J. Numer. Methods Heat Fluid Flow 2015, 25, 665–684. [Google Scholar] [CrossRef]
  35. Khan, W.A.; Khan, Z.H.; Rahi, M. Fluid flow and heat transfer of carbon nanotubes along a flat plate with Navier slip boundary. Appl. Nanosci. 2014. [Google Scholar] [CrossRef]
  36. Noor, M.A.; Mohyud-Din, S.T. Variational iteration technique for solving higher order boundary value problems. Appl. Math. Comput. 2007, 189, 1929–1942. [Google Scholar] [CrossRef]
  37. Abdou, M.A.; Soliman, A.A. New applications of Variational iteration method. Phys. D Nonlinear Phenom. 2005, 211, 1–8. [Google Scholar] [CrossRef]
  38. Turkyilmazoglu, M. Extending the traditional Jeffery Hamel flow tostretchable convergent/divergent channels. Comput. Fluids 2014, 100, 196–203. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Mohyud-Din, S.T.; Khan, U.; Ahmed, N.; Hassan, S.M. Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels. Appl. Sci. 2015, 5, 1639-1664. https://doi.org/10.3390/app5041639

AMA Style

Mohyud-Din ST, Khan U, Ahmed N, Hassan SM. Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels. Applied Sciences. 2015; 5(4):1639-1664. https://doi.org/10.3390/app5041639

Chicago/Turabian Style

Mohyud-Din, Syed Tauseef, Umar Khan, Naveed Ahmed, and Saleh M. Hassan. 2015. "Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels" Applied Sciences 5, no. 4: 1639-1664. https://doi.org/10.3390/app5041639

APA Style

Mohyud-Din, S. T., Khan, U., Ahmed, N., & Hassan, S. M. (2015). Magnetohydrodynamic Flow and Heat Transfer of Nanofluids in Stretchable Convergent/Divergent Channels. Applied Sciences, 5(4), 1639-1664. https://doi.org/10.3390/app5041639

Article Metrics

Back to TopTop