International Journal of Engineering Mathematics
Volume 2016 (2016), Article ID 8915745, 12 pages
http://dx.doi.org/10.1155/2016/8915745

Research Article

Perturbation Solutions for Hagen-Poiseuille Flow and Heat Transfer of Third-Grade Fluid with Temperature-Dependent Viscosities and Internal Heat Generation

B. Y. Ogunmola, A. T. Akinshilo, and M. G. Sobamowo

Department of Mechanical Engineering, University of Lagos, Akoka-Yaba, Lagos, Nigeria

Received 30 June 2016; Accepted 6 September 2016

Academic Editor: Giuseppe Carbone

Copyright © 2016 B. Y. Ogunmola et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Regular perturbation technique is applied to analyze the fluid flow and heat transfer in a pipe containing third-grade fluid with temperature-dependent viscosities and heat generation under slip and no slip conditions. The obtained approximate solutions were used to investigate the effects of slip on the heat transfer characteristics of the laminar flow in a pipe under Reynolds’s and Vogel’s temperature-dependent viscosities. Also, the effects of parameters such as variable viscosity, non-Newtonian parameter, viscous dissipation, and pressure gradient at various values were established. The results of this work were compared with the numerical results found in literature and good agreements were established. The results can be used to advance the analysis and study of the behavior of third-grade fluid flow and steady state heat transfer processes such as those found in coal slurries, polymer solutions, textiles, ceramics, catalytic reactors, and oil recovery applications.

1. Introduction

The third-grade fluid being a subclass of non-Newtonian fluid whose viscosity varies based on the applied force is a favored fluid due to the exciting phenomena it captures such as the shear thinning and thickening effects. The mathematical model of a third-grade fluid represents a more realistic description of the behavior of non-Newtonian fluids. The model also represents a further attempt towards the study of the flow structure of non-Newtonian fluids. Consequently, considerable interest has been shown in the third-grade fluid over the past few decades by various researchers due to its potential application in industry and technology. It should be stated that the governing equations for the third-grade fluid model are nonlinear and much more complicated than those of Newtonian fluids. They require additional boundary conditions to obtain a physically meaningful solution. This issue has been discussed in detail by researchers [1–23] who made a complete thermodynamical analysis of a third-grade fluid and derived the restriction on the stress constitutive equation. They investigated some stability characteristics of third-grade fluids and showed that they exhibit features different from those of Newtonian and second-grade fluids. In an attempt to improve the characteristics of the viscoelastic properties of the fluid, the stability of the third-grade fluid model was studied by Fosdick and Rajagopal [1] while Majhi and Nair [2] investigated the effects of stenotic geometry and the non-Newtonian parameter of the third-grade fluid on the resistive impedance and wall shearing stress. Their results were compared with similar study submitted by Massoudi and Christie [3] who presented the numerical solutions on the effects of variable viscosity and viscous dissipation on the flow of third-grade fluid in a pipe using the finite difference method. Yürüsoy and Pakdemirli [4] developed approximate analytical solution for the flow of the third-grade fluid in a pipe using constant viscosity model and variable viscosity models under no slip condition while Vajravelu et al. [5] presented a numerical solution for the third-grade fluid flows between rotating cylinder using Schauder theory and perturbation technique. The fluctuating flow of a magnetohydrodynamic rotational flow of the third-grade fluid on a porous plate was studied by Hayat et al. [6]. The similarity solutions to boundary layer equation for the third-grade fluid were developed by Muhammet [7] using the special coordinate system generated by the potential flow. The steady flow analysis of the third-grade fluid between circular concentric cylinders with heat transfer was presented by Yürüsoy [8]. In the study, the pipe temperature is assumed to be higher than fluid temperature. Also, Pakdemirli and Yilbas [9] developed approximate analytical solution of non-Newtonian fluid using Vogel’s viscosity model and entropy generation in a pipe. The steady flow of a third-grade fluid past a porous horizontal plate with partial slip was investigated by Sajid et al. [10]. Ellahi et al. [11] used implicit finite difference method to analyze the unsteady free convective flow of a third-grade fluid past an infinite vertical plate when uniform suction is applied at the plate. Jayeoba and Okoya [12] presented the analytical approximate solution to determine the temperature fields for steady flow of a third-grade fluid in a pipe with models of viscosities including a heat generation term for the no slip boundary condition. Most of the studies previously carried out on the third-grade fluid are limited to no slip flow condition which is a simplified method of predicting the actual behavior of the fluid in various applications. In practice, most problems of fluid flow exist as either partial slip or slip condition, Hence, in this work, regular perturbation method is used to develop approximate analytical solutions that are used to investigate the effects of slip on the heat transfer analysis in pipe flow of third-grade fluid with temperature-dependent viscosities and internal heat generation.

2. Model Development and Analytical Solution

Consider the flow of a third-grade fluid in a pipe as shown in Figure 1. The momentum and temperature equations with incorporated source term are given by [3, 4, 13]The slip condition at the pipe wall can be introduced in terms of shear stress which makes the boundary condition nonlinear. Linearizing the boundary condition at the wall gives a corresponding linear boundary condition at the wall under slip condition as At the center of the pipe, the boundary condition is given asThe boundary conditions under slip condition are given as With the dimensionless parameters as stated in the nomenclature, the dimensionless equations (leaving out the bars on the equations) for (3) and (4) are given as And the dimensionless boundary conditions under slip condition are given as where and

Figure 1: Physical model and coordinate system [12].

In this work, two equations of temperature-dependent viscosity are used, Reynolds and Vogel’s models. Reynolds model is a model of viscosity varying constantly with temperature while Vogel’s model of viscosity is varying exponentially with temperature. The nondimensionalized form of both models yields

2.1. Reynolds Model of Viscosity

The regular perturbation technique is used to determine the approximate analytical solution for the Reynolds viscosity model.

Taking Maclaurin’s series, (13a) takes the formSeries solutions of the equation of the velocity and temperature may be obtained by using perturbation method taking as the perturbation parameter. Solutions are obtained in the following form:Substituting (15) into (9) and changing parameters into terms yieldAlso, substituting (14) and (15) into (10) and changing parameters into terms yieldThe boundary conditions for the leading order equation arewhere , , .

With the boundary conditions, it could be easily shown that (16) and (18) are given asAlso, for the first-order equation, the boundary conditions are given as where and

Substituting the corresponding terms from the solutions in (21) into (17) and (19) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (22) yieldsSubstituting (21) and (23) into (15) the series solution at the order of expansion and then changing terms back to original parameters finally give

2.2. Vogel’s Model of Viscosity

With the aid of Maclaurin series, Vogel’s viscosity model can be written asTaking the series solution of the velocity and temperature fields yields the expansion Therefore substituting (26) and (27) into (9) and changing parameters into term yieldwhere ,where .

Substituting (26), (27), and (28) into (10) and changing parameters into term yieldwhere ,  ,  ,  , The boundary conditions for the leading order are as follows:Using the boundary condition, it could be easily shown that (29) and (30) are given asAlso, the first-order boundary condition is as follows:where and .

Substituting the corresponding terms from the solutions in (33) and (34) into (30) and (32) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (36) yieldsSubstituting (34) and (37), (35), and (38) into (27) and (28), the series solution at the order of expansion, and then changing terms back to original parameters finally givewhere

3. Results and Discussion

Plotting the results for the temperature distribution for both Reynolds and Vogel viscosity models of Jayeoba and Okoya [12] who studied the flow and heat transfer under no slip and no temperature jump conditions and that of the temperature distribution obtained in this work under slip and temperature jump conditions are reported graphically below at varying parameters.

Figure 2 shows the effect of pressure drop parameter for the Reynolds viscosity model on the temperature distribution. From the result, it shows that as becomes more negative for no slip condition, the maximum temperature of the fluid which is at the center of the pipe increases, while for the slip condition the maximum temperature increases but the fully developed temperature profile does not begin at the origin which was shown in Figure 3.

Figure 2: Effect of pressure gradient parameter on the temperature distribution of Reynolds viscosity model under no slip condition when .

Figure 3: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under slip condition when .

When heat generation term is excluded, the effect is seen as . For the Reynolds viscosity model as seen in Figures 4 and 5, as increases the temperature distribution of the fluid increases and the effect of the heat generation on the fluid is maximum at the center of the pipe for no slip condition, while for slip condition the effect of increasing causes an increase in temperature distribution but the temperature profile does not begin at the origin. Increasing values of viscous dissipation gives a corresponding increase in temperature distribution for the Reynolds viscosity model because mechanical energy is converted into thermal energy. When the condition is no slip, is maximum at the pipe center as depicted in Figure 6. However, when the flow is under the slip, temperature profile is not from the origin as shown in Figure 7. The third-grade fluid can be seen to exhibit Newtonian character for the Reynolds viscosity model when and also the effect when the fluid behaves in non-Newtonian manner at increasing values of . It can be seen from Figure 8 that, at increasing values of , the temperature distribution of the fluid decreases under the no slip condition and the maximum effect of occurs at the center of the pipe, while for slip condition at increasing values of the temperature distribution decreases and the fully developed temperature profile does not begin at the origin which can be seen in Figure 9. Viscosity variation parameter for the Reynolds viscosity model can be viewed from Figures 10 and 11. As viscosity variation parameter increases, the temperature distribution increases for no slip condition and the effect of is maximum at the center of the pipe, in Figure 10. As increases for the slip condition, the temperature distribution increases and the fully developed profile does not begin at the origin as seen from Figure 11.

Figure 4: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under no slip condition when .

Figure 5: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under slip condition when .

Figure 6: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under no slip condition when

Figure 7: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under slip condition when .

Figure 8: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under no slip condition when .

Figure 9: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under slip condition when .

Figure 10: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under no slip condition for

Figure 11: Effect of varying parameter for the temperature distribution of Reynolds viscosity model under slip condition when .

Increasing values of parameter in Vogel’s viscosity model gives decreasing temperature distribution and the effect of Vogel’s parameter is maximum at the center of the pipe for the no slip condition shown in Figure 12. But the temperature distribution profile does not begin at the origin for slip condition which can be seen in Figure 13. Also, for the Vogel viscosity parameter the influence of increasing value of shows a corresponding increase in temperature distribution, where the maximum temperature occurs at the center of the pipe. Fully developed temperature profile begins from the origin which can be glanced in Figure 14. And for slip condition, the temperature distribution profile is not from the origin as in Figure 15. As pressure drop parameter () becomes more negative for Vogel’s viscosity model, the temperature distribution increases with the maximum temperature at the center of the pipe and temperature profile begins at the origin for no slip condition as shown in Figure 16, while effects of pressure drop parameter on the the temperature distribution for slip condition are shown in Figure 17. The fully developed temperature profile is not from the origin. When the heat generation term is excluded, the effect is seen as for Vogel’s viscosity model. When the condition is no slip, temperature profile begins at the origin and effect of heat generation term is maximum at the center of the pipe as shown in Figure 18. And when the condition is slip, the fully developed profile does not begin at the origin, the effect of heat generation is also maximum at the center, and the temperature distribution increases as shown in Figure 19.

Figure 12: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when .

Figure 13: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when .

Figure 14: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when

Figure 15: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when .

Figure 16: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when .

Figure 17: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when

Figure 18: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when .

Figure 19: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when .

The viscous dissipation parameter () for Vogel’s viscosity model shows that, for increasing values of , the temperature distribution increases. The effect of is maximum at the pipe center and the temperature profile begins at the origin for no slip condition which can be seen from Figure 20, while for slip condition the effect of is maximum at the center of the pipe and the fully developed temperature profile does not begin from the origin as shown in Figure 21. Increasing values of non-Newtonian parameter () for Vogel’s model gives decreasing values of temperature distribution where maximum effect occurs at the center of the pipe. The fluid behaves in a Newtonian manner when , as can be seen from Figure 22. When condition is slip, the effect of increasing values of gives a corresponding decrease in temperature distribution and the fully developed temperature profile does not begin at the origin which can be seen in Figure 23. When the initial temperature term for Vogel’s viscosity model is increased, there is a corresponding increase in temperature distribution. The effect of is maximum at the center of the pipe. When condition is no slip, the fully developed temperature profile begins from the origin as shown in Figure 24 and for slip condition increase in gives a corresponding increase in temperature distribution with the fully developed temperature profile not beginning from the origin, as depicted in Figure 25. From the velocity distribution in Figure 26, as the Knudsen number increases from zero, it can be seen that the velocity profile reduces, with the arrow indicating the direction of increasing Knudsen number. Thus, as the Knudsen number increases, the profile of the velocity decreases.

Figure 20: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when

Figure 21: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when for .

Figure 22: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when

Figure 23: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when for .

Figure 24: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under no slip condition when .

Figure 25: Effect of varying parameter for the temperature distribution of Vogel’s viscosity model under slip condition when .

Figure 26: Velocity distribution profile as the Knudsen number increases from 0 to 0.08.

4. Conclusion

In this work, approximate analytical solutions for the nonlinear analysis of heat transfer in a pipe flow of a third-grade fluid with temperature-dependent viscosities and heat generation under slip and no slip conditions have been developed using regular perturbation technique. The obtained analytical solutions were used to investigate the effects of slip boundary conditions on the heat transfer in the pipe flow. Also, the effects of parameters such as variable viscosity, non-Newtonian parameter, viscous dissipation, and pressure gradient at various values (which indicate different fluid flow situations) are established. The results can be used to advance the analysis and study of the behavior of third-grade fluid flow and steady state heat transfer processes such as those found in coal slurries, polymer solutions, textiles, ceramics, catalytic reactors, and oil recovery applications.

Nomenclature

and : | Constants from dimensional Vogel’s viscosity model
: | Constant from dimensional Reynolds viscosity model
: | Dimensionless constant from Vogel’s viscosity model
: | Dimensionless constant from Vogel’s viscosity model
: | Initial concentration of reacting species
: | Pressure gradient parameter
: | Constant thermal conductivity
: | Pressure gradient along the normal to the pipe axis
: | Pressure gradient in the axial direction
: | Pressure gradient in rotational direction
: | Heat generation constant
: | Dimensional perpendicular distance in the pipe axis
: | Dimensionless perpendicular distance in the pipe axis
: | Radius of the pipe
: | Initial temperature
: | Dimensional velocity component in the -axis
: | Dimensionless velocity component in the -axis
: | Dimensional reference velocity
: | Axis of the cylinder
: | Material constants
: | Activation energy parameter
: | Reynolds viscosity variational parameter
: | Dynamic shear viscosity
: | Dimensionless viscosity
: | Viscous heating parameter
: | Dimensionless temperature excess
: | Non-Newtonian material parameter of the fluid
: | Heat generation parameter.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

R. L. Fosdick and K. R. Rajagopal, “Thermodynamics and stability of fluids of third grade,” Procter Society London, vol. 339, pp. 351–377, 1980.

S. N. Majhi and V. R. Nair, “Flow of a third grade fluid over a sternosed tubes,” Indian National Science Academy, vol. 60, no. 3, p. 535, 1994.

M. Massoudi and I. Christie, “Effects of variable viscosity and viscous dissipation on the flow of a third grade fluid in a pipe,” International Journal of Non-Linear Mechanics, vol. 30, no. 5, pp. 687–699, 1995.

M. Yürüsoy and M. Pakdemirli, “Approximate analytical solutions for the flow of a third-grade fluid in a pipe,” International Journal of Non-Linear Mechanics, vol. 37, no. 2, pp. 187–195, 2002.

K. Vajravelu, J. R. Cannon, D. Rollins, and J. Leto, “On solutions of some non-linear differential equations arising in third grade fluid flows,” International Journal of Engineering Science, vol. 40, no. 16, pp. 1791–1805, 2002.

T. Hayat, S. Nadeem, S. Asghar, and A. M. Siddiqui, “Fluctuating flow of a third-grade fluid on a porous plate in a rotating medium,” International Journal of Non-Linear Mechanics, vol. 36, no. 6, pp. 901–916, 2001.

Y. Muhammet, “Similarity solutions to boundary layer equations for third grade non-Newtonian fluid in special channel coordinate system,” Journal of Theoretical and Applied Mechanics, vol. 41, no. 4, pp. 775–787, 2003.

M. Yürüsoy, “Flow of a third grade fluid between concentric circular cylinders,” Mathematical and Computational Applications, vol. 9, no. 1, pp. 11–17, 2004.

M. Pakdemirli and B. S. Yilbas, “Entropy generation for pipe flow of a third grade fluid with Vogel model viscosity,” International Journal of Non-Linear Mechanics, vol. 41, no. 3, pp. 432–437, 2006.

M. Sajid, R. Mahmood, and T. Hayat, “Finite element solution for flow of a third grade fluid past a horizontal porous plate with partial slip,” Computers & Mathematics with Applications, vol. 56, no. 5, pp. 1236–1244, 2008.

R. Ellahi, T. Hayat, F. M. Mahomed, and S. Asghar, “Effects of slip on the non-linear flows of a third grade fluid,” Nonlinear Analysis: Real World Applications, vol. 11, no. 1, pp. 139–146, 2010.

O. J. Jayeoba and S. S. Okoya, “Approximate analytical solutions for pipe flow of a third grade fluid with variable model of viscosities and heat generation/absorption,” Journal of the Nigerian Mathematical Society, vol. 31, pp. 207–227, 2012.

S. S. Abbasbandy, T. Hayat, R. Ellahi, and S. Asghar, “Numerical results of a flow in a third grade fluid between two porous walls,” Zeitschrift für Naturforschung A, vol. 64, no. 1-2, pp. 59–64, 2009.

K. S. Adegbie and F. I. Alao, “Flow of a newtonian fluid in a symmetrically heated channel: effect of viscosity and viscous dissipation,” Mathematical Problems in Engineering, vol. 2006, Article ID 29314, 7 pages, 2006.

Y. M. Aiyesimi, G. T. Okedayo, and O. W. Lawal, “Unsteady magnetohydrodynamic (MHD) thin film flow of a third grade fluid with heat transfer and no slip boundary condition down an inclined plane,” International Journal of Physical Sciences, vol. 8, no. 19, pp. 946–955, 2013.

R. P. Chhabra, Non Newtonian Fluid: An Introduction, Indian Institute of Technology, New Delhi, India, 2010.

R. P. Chhbra and J. F. Richardson, Non Newtonian Flow and Rheology, Butterworth-Heinemann, Oxford, UK, 2nd edition, 1999.

R. Ellahi, M. Raza, and K. Vafai, “Series solutions of non-Newtonian nanofluids with Reynolds' model and Vogel's model by means of the homotopy analysis method,” Mathematical and Computer Modelling, vol. 55, no. 7-8, pp. 1876–1891, 2012.

T. Hayat, Z. Iqbal, M. Mustafa, and A. A. Hendi, “Melting heat transfer in the stagnation-point flow of third grade fluid past a stretching sheet with viscous dissipation,” Thermal Science, vol. 17, no. 3, pp. 865–875, 2013.

T. Hayat, A. H. Kara, and E. Momoniat, “The unsteady flow of a fourth-grade fluid past a porous plate,” Mathematical and Computer Modelling, vol. 41, no. 11-12, pp. 1347–1353, 2005.

M. Kothandapani, “The influence of exponential viscosity on a MHD viscoelastic fluid over a stretching sheet,” International Journal of Applied Mathematics and Mechanics, vol. 9, no. 5, p. 81, 2012.

S. Nadeem and M. Ali, “Analytical solutions for pipe flow of a fourth grade fluid with Reynolds and Vogel’s model of viscosities,” Communications in Nonlinear Mechanics, vol. 43, pp. 558–599, 2009.

S. Nadeem and H. Sadia, “Influence of temperature dependent viscosity and entropy generation on the flow of a Johnson-Segalman fluid,” Walailak Journal, vol. 10, no. 5, pp. 553–557, 2013.