Free convection occurs due to the impact of gravity on the heated fluid of variable density and thus relative buoyancy of the fluid. It occurs at small velocities in the presence of large temperature difference. Example of the free convection is a pot of boiling water in which the hot and less-dense water on the bottom layer moves upwards in plumes while the cool and denser water near the top of a pot likewise sinks. The fully developed laminar free convective flow of a fluid heated asymmetrically between two vertical plates has been described very well by Aung  . Nelson and Wood  have obtained the solution in the case of the combined influence of heat and mass transfer on natural convection in vertical channels.
The heat transfer phenomena in non-Newtonian fluid flows through ducts are found in analysis of the polymer in chemical engineering and applications of fluid film lubrication etc. The behaviors and theories of many non-Newtonian fluids have been proposed by many researchers. The concepts of such fluids demand a new idea which differs from the real fluids. Here, we shall study a special family of non-Newtonian fluid such as polar fluid, whose substantial equations with couple stresses have been described by Aero et al.  , D’ep  and Cowin  . Some examples of a polar fluid are mud, crude oil, body fluids, lubricants with polymer additives etc. Comparative studies of the Poiseuille flow corresponding to a polar fluid under various boundary conditions with applications to blood flow have been studied by Chaturani and Biswas  . Singh et al.  have studied the natural convection under a radial magnetic field in vertical concentric annuli. Using the generalized Forchheimer’s model, Raptis and Takhar  have studied the flow of a polar fluid in a porous medium bounded by an infinite plate. Ogulu  has obtained the analytical solutions for the problem concerned with free convective flow of a polar fluid past a vertical porous plate with couple stresses and radiation, where the thermal energy of the plate is assumed to oscillate about a mean value. Later, the studies of flow formation in different geometries with and without heat source have been attempted by Cheng   . Patil  has obtained analytic solution for the problem of free convection on the oscillatory flow of a polar fluid within a porous medium with variable heat flux.
Fluid flow with internal heat sources/sinks finds physical as well as theoretical interest. The fluid motion develops slowly following the development of non-uniformity in the temperature field. However, when the thermal energy differences are admirably large, the volumetric heat generation or absorption term exerts strong impact on the heat transfer and as a consequence on the fluid flow as well. The natural convective flows with heat generation from an isothermal horizontal circular cylinder have been studied by Molla et al.  . The effects of chemical reaction on free convective flow of a polar fluid in a porous medium with internal heat generation have been described by Patil and Kulkarni  . The natural convective flows of a polar fluid between vertical concentric annuli have been examined by Ravi and Singh  . Mohamed et al.  have obtained the solutions for thermal radiation and magnetic effects on free convective flow of a polar fluid in a porous medium with internal heat generation and chemical reaction. Moreover, Ravi et al.  have presented the free convection of a micropolar fluid between two vertical walls by taking into account the temperature dependent source/sink. Perdikis and Rapti  have studied the flow of an incompressible and electrically conducting polar fluid in the presence of a magnetic field. Motivated by the above papers, a boundary value problem is considered here to investigate analytically the effect of heat source/sink on the velocity and microrotational velocity when outer cylinder is isothermal while thermal boundary condition at the inner one is mixed kind. Analytical solutions obtained here are compared to the solution of natural convection of a polar fluid in vertical concentric annuli by Ravi and Singh  to check the accuracy of the solution. Finally, the effects of heat source or sink, vertex viscosity and radii ratio parameters have been shown by graphs.
2. Mathematical Formulations
We have considered steady and laminar free convective flow of a polar fluid inside vertical concentric annulus of infinite length with temperature dependent heat source/sink. The -axis is taken along the axis of cylinder in the vertical upward direction and -axis is in the radial direction. Let a and b be the radii of inner and outer cylinders respectively. In this case the inner cylinder is cooled or heated either at a constant heat flux or isothermally so that its temperature (i.e. temperature of outer surface of the inner cylinder) is different from the ambient temperature . Since cylinders are of infinite length, the flow depends only on . Using the Boussinesq approximation, the governing equation for free convective flow of a polar fluid with temperature dependent heat source or sink can be obtained as follows:
The boundary conditions for the model are as follows:
For the non-dimensionalization of governing Equations (1) - (3) and boundary conditions given by (4), we use the non-dimensional variables and parameters for the considered flow, defined as follows:
In this equation, is or according as inner cylinder is
maintained at constant temperature or constant heat flux . Use of Equation (5) into Equations (1) - (3) has produced the following non-dimensional equations:
Further, the corresponding boundary conditions in dimensionless form are obtained as
The boundary conditions for the temperature field at the inner cylinder are of the Dirichlet and Newmann types. For a unified solution for both boundary conditions, these boundary conditions for temperature field at the inner cylinder can be combined into a single boundary condition known as boundary condition of third kind at the inner cylinder. By doing so, a combined condition is obtained as follows:
The corresponding cases can be obtained by assigning appropriate values to A, B and C in Equation (10). The solution of Equations (6) - (8) subject to their suitable boundary conditions are obtained as for the cases of source and sink respectively.
Case 1 When source is present
In this case, and as a result of which the solutions for T, u and N are obtained as follows:
Using Equation (10), the dimensionless skin-frictions at outer surface of inner cylinder and inner surface of outer cylinder in the presence of source are derived as
Case 2 When sink is present
In the case of sink, and to find analytical solution, we take . Corresponding to this case, the solutions for T, u and N in terms of Si are derived as follows:
In addition, the dimensionless skin-frictions at outer surface of inner cylinder and inner surface of outer cylinder with sink are given by
3. Results and Discussion
The impacts of physical parameters such as the source (S)/sink (Si), vertex viscosity (R) and radii ratio (λ) on the velocity (u), microrotational velocity (N) are analyzed by presenting the numerical calculation of the analytical solution for constant heat flux case and isothermal case on the graphs.
Graphical representation in Figures 1-3 shows the variation of the velocity profiles with source, vertex viscosity and gap between the cylinders for both cases such as constant heat flux and isothermal. It is clear from these figures that the velocity of the fluid increases when the source and gap between the cylinders increases for both cases. When , the velocity is greater for the constant heat flux case compared to the isothermal case while reverse phenomenon occurs when . In light of this effect, we obtain the critical value of at which the velocity profiles are the same for both cases. The nature of the velocity of the fluid is almost the same for constant heat flux and isothermal cases when the critical value of . The critical values 2.5 and 2.18 of corresponding to and 0.4, clearly show that the effect of source parameter is to decrease it. As in the case above, the velocity of the fluid decreases with increasing vertex viscosity parameters for both cases.
Figures 4-6 presents the variation of the microrotational velocity profiles with the source, vertex viscosity and gap between the cylinders for both the cases of
Figure 1. Velocity profiles for .
Figure 2. Velocity profiles for .
Figure 3. Velocity profiles for .
Figure 4. Influences of source and vertex viscosity parameters on microrotational velocity for .
Figure 5. Effect of source and vertex viscosity parameters on microrotational velocity for .
Figure 6. Effect of source and vertex viscosity parameter on microrotational velocity for λ < 2.5.
thermal heating. These figures demonstrate that the effect of the source parameters and gap between the cylinders on the microrotational velocity is to increases it for both cases. When , the microrotational velocity is more for the constant heat flux case compared to the isothermal case while reverse trend occurs when . Again, we reveal that the microrotational velocity profiles of a polar fluid for are approximately the same for both cases and the critical values of decrease with source ( ). The microrotational velocity of the fluid decreases with increasing vertex viscosity parameter for both types of thermal conditions. Figure 7 and Figure 8 describe that the velocity and microrotational velocity of a polar fluid with sink are the same for constant heat flux and isothermal cases when . In this case, we find that the critical value of for the velocity and microrotational velocity of a polar fluid enhances with heat sink for constant heat flux and isothermal cases. As in the case above, the velocity and microrotational velocity of the fluid decreases with increasing vertex viscosity parameter for both cases.
The numerical values of skin-friction and for isothermal and constant heat flux cases with temperature dependent heat source/sink are given in Table 1. This table clearly reveals that the skin-friction on inner cylinder increases when heat source increases and decreases with increasing heat sink for
Figure 7. Velocity profiles for .
Figure 8. Impact of sink and vertex viscosity parameters on microrotational velocity for .
Table 1. Numerical values of skin-friction for Isothermal and constant heat flux cases.
isothermal and constant heat flux cases while reverse tendency occurs on outer cylinder. The skin-friction on both the cylinders increases with increasing values of (gap between the inner and outer cylinder) for both cases. In Table 1, the comparative study of isothermal and constant heat flux cases shows that the skin-friction and are greater for isothermal heating of the inner cylinder than constant heat flux case when the gap between cylinders are less than and 3.26 times the radius of the inner cylinder with heat source/sink, while the effect is just reverse when the gap between the cylinders are greater than and 3.26 times the radius of inner cylinder respectively. The numerical values of skin-frictions on both cylinders are almost same at and 3.26 when and −0.1.
The analytic solutions of the governing differential equations in terms of the source/sink, vertex viscosity and radii ratio parameters have been obtained. We have examined the impact of source/sink and radii ratio parameters on the velocity and microrotational velocity of a polar fluid by taking mixed type of thermal condition on the surface of inner cylinder. It is established that the gap between cylinders and source or sink parameter plays an important role in controlling the properties of fluid flow. The impacts of the source parameter on the velocity and microrotational velocity have increasing tendency while decreasing tendency is due to the sink. The influence of thermal boundary conditions on the inner cylinder strongly depends on the gap between the cylinders and this result can be used for modeling in such types of system. In this case, we obtain that the skin-friction of a polar fluid on outer cylinder decreases with heat source and increases with sink for isothermal and constant heat flux cases while reverse phenomenon occurs on inner cylinder.
The author is thankful to UGC, New Delhi, India for the Economical endorsement in the form of a Junior Research Fellowship to finish this work.
Inner radius of the cylinder
Outer radius of the cylinder
Dimensionless source parameter
Dimensionless sink parameter
Dimensionless angular velocity
Stream wise velocity
Dimensionless stream wise velocity
#Math_91# Dimensionless transverse coordinate
Constant heat flux
Rate of heat flux
Temperature of the fluid
Temperature of the wall
Temperature of the free steam
Specific heat at a constant pressure
Acceleration due to gravity
Vertex viscosity parameter
Density of fluid
Coefficient of thermal expansion
Ratio of outer radius and inner radius, b/a
Dynamic viscosity of the fluid
Skin friction at the wall
Skin friction at the wall
Kinematic viscosity of the fluid
Submit or recommend next manuscript to SCIRP and we will provide best service for you:
Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.
A wide selection of journals (inclusive of 9 subjects, more than 200 journals)
Providing 24-hour high-quality service
User-friendly online submission system
Fair and swift peer-review system
Efficient typesetting and proofreading procedure
Display of the result of downloads and visits, as well as the number of cited articles
Maximum dissemination of your research work
Submit your manuscript at: http://papersubmission.scirp.org/
Or contact email@example.com
 Aung, W. (1972) Fully Developed Laminar Free Convection between Vertical Plates Heated Asymmetrically. International Journal of Heat and Mass Transfer, 15, 1577-1580.
 Nelson, D.J. and Wood, B.D. (1973) Combined Heat and Mass Transfer Natural Convection between Vertical Parallel Plates. International Journal of Heat and Mass Transfer, 32, 1779-1787. https://doi.org/10.1016/0017-9310(89)90059-8
 Chaturani, P. and Biswas, D. (1984) A Comparative Study of Poiseuille Flow of a Polar Fluid Under Various Boundary Conditions with Applications to Blood Flow. Rheologica Acta, 23, 435-445.
 Singh, S.K., Jha, B.K. and Singh, A.K. (1997) Natural Convection in Vertical Concentric Annuli under a Radial Magnetic Field. Heat and Mass Transfer, 32, 399-401.
 Ogulu, A. (2005) On the Oscillating Plate Temperature Flow of a Polar Fluid past a Vertical Porous Plate in the Presence of Couple Stresses and Radiation. International Communications in Heat Mass Transfer, 32, 1231-1243.
 Cheng, C.Y. (2006) The Effect of Temperature-Dependent Viscosity on the Natural Convection Heat Transfer from a Horizontal Isothermal Cylinder of Elliptic Cross Section. International Communications in Heat and Mass Transfer, 33, 1021-1028.
 Cheng, C.Y. (2009) Natural Convection Heat Transfer from a Horizontal Isothermal Elliptical Cylinder with Internal Heat Generation. International Communications in Heat and Mass Transfer, 36, 346-350.
 Patil, P.M. (2008) Effects of Free Convection on the Oscillatory Flow of a Polar Fluid through a Porous Medium in the Presence of Variable Heat Flux. Journal of Engineering Physics and Thermo Physics, 81, 905-922.
 Molla, M.M., Hossain, M.A. and Paul, M.C. (2006) Natural Convection Flow from an Isothermal Horizontal Circular Cylinder in Presence of Heat Generation. International Journal of Engineering Science, 44, 949-958.
 Patil, P.M. and Kulkarni, P.S. (2008) Effects of Chemical Reaction on Free Convective Flow of a Polar Fluid through a Porous Medium in the Presence of Internal Heat Generation. International Journal of Thermal Sciences, 47, 1043-1054.
 Mohamed, R.A., Abo-Dahab, S.M. and Nofal, T.A. (2010) Thermal Radiation and MHD Effects on Free Convective Flow of a Polar Fluid through a Porous Medium in the Presence of Internal Heat Generation and Chemical Reaction. Mathematical Problems in Engineering, 2010, 27. https://doi.org/10.1155/2010/804719
 Ravi, S.K., Singh, A.K. and Alawadhi, K.A. (2011) Effect of Temperature Dependent Heat Source/Sink on Free Convective Flow of a Micropolar Fluid between Two Vertical Walls. International Journal of Energy and Technology, 3, 1-8.