Free convective heat transfers inside concentric and eccentric annulus has received an intensive attention of researchers owing to its scientific, technological and engineering applications such as in the construction of electrical motors and generators, cooling of electronic components, nuclear reactors, thermal storage systems, heating and cooling of underground electric cables, completion of oil source and aircraft fuselage insulation.  considered transient natural convection heat transfer problem between two horizontal isothermal cylinders formed within the Boussinesq approximation and solved numerically through the vorticity-stream function approach. Numerical prediction of natural convection heat transfer in horizontal annulus where the inner cylinder is hotter than outer was reported by  .  investigated natural convection of gasses in a horizontal annulus, where the inner cylinder is heated by the application of a constant heat flux and the outer cylinder is isothermally cooled.  conducted an investigation on natural convection in concentric and eccentric horizontal cylindrical annuli with mixed boundary conditions. In  , natural convection in a narrow horizontal cylindrical annulus for fluid of been investigated. Also  discussed the existence of dual solution and role of Prandtl number on bifurcation.  performed three-dimensional linear stability analysis of air in horizontal concentric annuli between concentric cylinders by using a spectral method. In  , finite difference method in conjunction with least-squares scheme and experimental temperature data is used to predict the average heat transfer coefficient and fin efficiency on the fin of annular-fined tube.  carried out experimental investigation on heat transfer characteristics of Taylor-Couette- Poiseuille flow in an annular channel by mounting longitudinal ribs on the rotating inner cylinder.  studied natural convection flow in a horizontal annulus enclosure with a transversely oscillating inner cylinder.  investigated transition of natural convection in an annulus between horizontal concentric cylinders theoretically by assuming two-dimensional and incompressible flow fields.  examined a numerical investigation on horizontal concentric annulus with open ends and conditions of either adiabatic or isothermal outer cylinder surface.  showed that dual steady state solutions exist above a critical Rayleigh number for free convective flows in a horizontal annulus with constant heat-flux wall.  discussed the effect of nanofluid on the natural convection heat transfer and fluid flow through an annular tube with an inner heat generating solid circular rod.  compiled a comprehensive theoretical study on natural convection heat transfer in nanofluid contained inside the horizontal annular space existing between two long concentric cylinder whose surface are maintained at different uniform temperatures, with primary scope to determine the main heat transfer features for various operating conditions, nanoparticle diameters, and solid-liquid combination.  conducted a numerical investigation of natural convection heat transfer in a semi-annulus enclosure filled with nanofluid using the control volume based finite element method.  developed and tested a discrete phase model for forced convection of nanofluids in a circular tube subjected to a uniform heat flux.  and  examined the role of magnetic field on natural convection flow of Nano-fluid using Lattice Boltzmann method.  studied numerically natural convection flow and heat transfer of Copper water nanofluid inside eccentric horizontal annulus while the inner and outer cylinder are kept at constant temperatures.  conducted an investigation on effect of static radial magnetic field on natural convection heat transfer in a horizontal cylindrical annulus enclosure filled with nanofluid using Boltzmann method (LBM).
2. Mathematical Formulation
The physical problem under consideration consists of a transient natural convection flow in an infinite vertical annulus formed by two infinite concentric vertical cylinders in presence of thermal radiation. The transient flow formation is due to sudden heating of outer surface of inner cylinder in presence of thermal radiation. The physical properties of the working fluid are assumed to be constant. A schematic diagram of the present problem is shown in Figure 1. At time , both the fluid cylinder are assumed to be at thesame temperature with absence of fluid motion. At time , the temperature of the outer surface of inner cylinder at is suddenly raised to causing transient free convection current. Since the flow is fully developed and the cylinders are of infinite length, the flow depends on radial coordinate ( ) and time ( ). Using of Boussinesqs approximation, the governing equations in dimensional form are:
The radiation heat flux term in the problem is simplified by using the Rosseland approximation:
Figure 1. Schematic diagram for problem under consideration.
The initial and boundary conditions for the present problem assumed the form:
To obtain the solutions of Equations ((1) and (2)) subject to the conditions (4) in dimensionless form, the following appropriate dimensionless quantities are introduce to the problem:
Using the dimensionless quantities introduced in Equation (5), the dimensionless form of equation in (1) and (2) are:
while the dimensionless initial and boundary conditions are:
3. Analytical Solution
Analytical solution is often an opportunity to validate computer routines of complicated problems and comparison with data as well as inspecting the internal consistency of mathematical models. It is of interest to reduce the nonlinear governing equations presented in the previous section to a form that can be solved analytically. A special case of the present problem that can exhibits analytical solution is the problem of steady state free-convection flow in vertical annulus in presence of thermal radiation. The resulting steady state equations and the boundary conditions for the special case can be written as:
The relevant boundary conditions to be satisfied are:
We approximate solution to Equations ((9) and (10)) subject to (11) using regular perturbation method by taking power series expansion in the radiation parameter R.
Substituting Equation (12) into Equations ((9) and (10)) and equating the like power of R, the required analytical expressions for the steady state velocity and temperature fields subject to boundary conditions (11) are as follows:
From (13), the steady-state skin frictions at the outer surface of inner cylinder and inner surface of the outer cylinder of the annulus are:
Equally from (14), the steady state of heat transfers at the outer surface of inner cylinder and inner surface of outer cylinder of the annulus are:
4. Numerical Procedure
The momentum and energy equations given in Equations ((6) and (7)) are solved numerically using implicit finite difference method. The time derivatives in both equations are approximated using backward difference formula as:
while the first and second order space derivatives are approximated by the central difference formula.
Replacing by in (19), (20), (22) and (23) gives an iterative system, which does not restrict the time step. Thus the transport Equations ((6) and (7)) at the grid point are linearized. The momentum Equation reads:
5. Results and Discussion
The dimensionless mathematical model representing the transient natural convection flow in an annulus with thermal radiation is solved numerically using finite difference method. The problem examines two fluids air and water taking into account the influence of thermal radiation parameter R and temperature difference parameter . The numerical scheme is validated using the steady state analytical solution derived by perturbation method.Choosing small perturbation parameter , the results are found in good agreement between numerical and analytical solution at large value of time as depicted in Figure 2. The present parametric study has been performed over reasonable ranges of and with so as to capture the transient behavior of both velocity and temperature. Besides, all other parameters are taken arbitrarily. The result that clearly reports the influence of the flow governing parameters on velocity, temperature, skin friction and Nusselt number has been shown graphically in Figures 3-10.
Figure 3 represents transient and steady state velocity for different values of temperature difference parameter ( ) and time (t). It is observed that as the values of increases the velocity increases. Also as time increases the velocity increases and finally attains its steady state value. A comparison of Figure 3(a) and Figure 3(b) reveals that velocity is inversely proportional to Prandtl number of the fluid. This is physically true since penetration of heat is low in case of
Figure 2. Temperature and Velocity Profiles
Figure 3. Velocity profile
Figure 4. Velocity profile
higher value of Pr which create weak convection current inside the annulus.
Figure 4 depicts the variations of transient and steady state velocity for different values of the radiation parameter (R) and time (t). It is recorded from these Figure 4(a) and Figure 4(b) that, as R increases, velocity increases. The physical fact is that, an increase in radiation adds more heat to the fluid leading to an increase in convection current, which causes velocity increase. These figures also report that as time increases the velocity for both air and water increases and attains steady state. During the course of numerical computation, it was observed that the time required to reach steady state is directly proportional to Prandtl number of the fluid.
Figure 5(a) and Figure 5(b) illustrate the influence of temperature difference parameter ( ) and time (t) on temperature profile. It is observed from Figure 5 that as the value of is increased there is a corresponding increase in the temperature of the fluid. In addition temperature increases as dimensionless time increases and finally reach it steady state value.
The variation of skin friction for different values of time and temperature difference parameter ( ) is presented in Figure 7(a), Figure 7(b) and Figure 8(a), Figure 8(b) . From these figures, it is observed that as increases, skin friction increases at the outer surface of the inner cylinder
Figure 5. Temperature profile
Figure 6. Temperature profile
Figure 7. The variation of skin friction with and t: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 8. The variation of skin friction with and t: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
and inner surface of the outer cylinder respectively. However, the values of skin friction are higher at the outer surface of the inner cylinder in comparison to the skin friction at the inner surface of the outer cylinder . In addition, as time increases skin friction increases and finally reached its steady state value.
From Figure 9(a), Figure 9(b) and Figure 10(a), Figure 10(b) it is evident that skin friction at the outer surface of inner cylinder as well as the inner surface of outer cylinder increases with increase in radiation parameter (R). Also as time increases the skin friction increases and achieved its steady state value. It is interesting to mention that skin friction is high for small value of . This is consistent with the fact that as Prandtl number decreases velocity increases.
Figure 11(a), Figure 11(b) and Figure 12(a), Figure 12(b) show the variation of the Nusselt number with and (t) respectively. As illustrated in Figure 11(a) and Figure 12(a), the Nusselt number decreases with increase in and t at the outer surface of inner cylinder , while the result is just converse at the inner surface of the outer cylinder . From these figures it is also clear that rate of heat transfer is higher for air in comparison to water. In fact the values of Nusselt number in Figure 11(a), Figure 11(b) are significantly higher in comparison with Figure 12(a), Figure 12(b) for small values of dimensionless time t.
Figure 9. The variation of skin friction with R and t: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 10. The variation of skin friction with R and t: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 11. The variation of the Nusselt number with and: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 12. The variation of the Nusselt number with and: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 13(a), Figure 13(b) and Figure 14(a), Figure 14(b) shows the variation of Nusselt number for varying values of R and time t respectively. In Figure 13(a) and Figure 14(a), the Nusselt number is significantly higher for small values of R and decreases with increase in dimensionless time (t) to zero with increase in R at the outer surface of inner cylinder , while at the inner surface of the outer cylinder Nusselt number increases with both and as illustrated in Figure 13(b) and Figure 14(b).
The effect of temperature difference, thermal radiation and Prandtl number on
Figure 13. The variation of the Nusselt number with R and: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
Figure 14. The variation of the Nusselt number with R and: (a) at the outer surface of the inner cylinder and (b) at the inner surface of the outer cylinder when .
transient-natural convection flow inside vertical annulus is analysed analytically as well as numerically. The expression for velocity, temperature, skin-friction and Nusselt number for the nonlinear partial differential equation are obtained analytically by the well-known perturbation series method under steady state operating condition. The numerical solution is obtained by implicit finite difference method for transient situation. The outcome of the result shows that:
1) Velocity and temperature increase with increase of thermal radiation and temperature difference parameter.
2) During transient state, the maximum value of velocity occurs at smaller radial distance from the outer surface of inner cylinder and then decreases to zero asymptotically.
3) Thermal radiation and temperature difference parameter enhances the skin friction at the outer surface of inner cylinder as well as inner surface of the outer cylinder.
4) Nusselt number is higher at small values of and R while decreases with increase in dimensionless time (t) to zero with increase in and R at the outer surface of inner cylinder .
5) Nusselt number increases with increase in and R and dimensionless time (t) at the inner surface of outer cylinder .
6) During course of this investigation, our results are found in good agreement between steady state and transient solution after some sufficiently large time.
The authors wish to thank the referees for their valuable comments to improve the quality of the paper.
radius of the inner cylinder
radius of the outer cylinder
acceleration due to gravity
temperature difference parameter
thermal radiation parameter
temperature of the outer surface of inner cylinder
vertical component of velocity
vertical co-ordinate, direction of flow
coefficient of thermal expansion
mean absorption coefficient
density of the fluid