Numerical Simulation of Natural Air Convection in Inclined Eccentric Hemispheres Enclosure

The subject of natural convection heat transfer is motivated by a wide range of applications in engineering technology. The hemispherical cavity is a part of basic geometries although it is not widely studied. The effect of inclinaison on natural convection fluid motions in the gap between two eccentric hemispheres is numerically studied. The inner hemisphere is subjected to a heat flux of a constant density and the outer one is maintened isothermal. The walls separating the two hemispheres are thermally adiabatic. Equations are formulated with vorticity and stream-functions variables. It is also assumed the fluid incompressible and obeys the approximation of Boussinesq. These equations are written by using bispherical coordinates system and solved by using a finite difference method. The results show the topology of flow is strongly dependent on the inclinaison because the flow can change from a unicellular regime to a multicellular regime by varying the inclination from 0 to π. By increasing the Rayleigh number ( Ra < < ), the flow intensifies. The results are shown in terms of streamlines and isotherms during their transient evolution.

spheres surfaces is of considerable interest in many engineering areas. The spheres or hemispheres can be used as constructing and insulated surfaces The rectangular or square cavity has been extensively treated [1] [2] given their simplicity of implementation. Sarr et al. [3] and Waheed [4] studied the cylindrical annuli.
Several analytic and numerical studies dealing with spherical sections have been presented by [5] [6]. However, very little works are done on natural convection between two hemispheres. Cabelli [7] examined the effect of inclinaison on the convective motion in hemispherical cavities. They showed that the axially symmetric motion which exists when the base of the hemisphere is horizontal becomes three-dimensional when the hemisphere is inclined. Shiina et al. [8] performed an experimental study of natural convection in a hemisphere that is differentially heated for nuclear applications. The heat transfer is empirically correlated by means of Nusselt-Rayleigh type. By studying theorical and experimental natural convection heat transfer from isothermal, Lewandowski et al. [9] presented their result in theorical consideration is presented in the Nusselt-Rayleigh correlation which the exponent of Rayleigh number is equal to a quarter. Saber et al. [10] numerically analysed the 2D convective in hemispherical cavities with planar inner surfaces. They observed that for differents boundaries temperatures conditions, the airflow in the cavities is monocellular and reaches steady state conditions for both cold and hot weather conditions. The numerical results are used to develop practical correlations for the Nusselt number. Bairi et al. [11] dealt with isothermal disks and dome faced downwards. Their results showed that the natural convection heat transfer is between 13% et 21% higher when the disk is maintened as compared to the case corresponding to imposed heat flux on the disk, for the same inclinaison angle and Rayleigh range. We note that all these hemispherical cavities mentioned have an enclosure formed by a dome and a horizontal plane. Little work has been done when the enclosure is delimited by two hemispheres. Thus, Alilat [12] studied the quantification of natural convection between two inclined concentric hemispheres for thermoregulation applications in electronics. His results showed that the inclinaison angle of the device affects only weakly the heat exchange and the natural convection heat transfers are also quantified by means of the Nusselt number. Recently, by studying natural convection between two eccentric hemispheres, Koita et al., [13] observed that the parameters of fluid depend strongly on the modified Rayleigh number as well as of the eccentricity. Also, the effects of eccentricity showed the center of rotation moves toward the top of enclosure when the eccentricity increases and the Nusselt number increases with the modified Rayleigh number. According to the authors' knowledge, transient natural convection between two eccentric inclined hemispheres has not been the subject of any study. The present numerical study is motivated by the interest to understand the effects of inclinaison, the eccentricity and the modified Rayleigh number on the convection motion over time between two hemispheres. The inner hemisphere is heated by a heat flux constant density and the outer one is maintened isothermal. The walls separating the two hemispheres are insulated. Open Journal of Applied Sciences

Mathematical Formulation
A schematic depiction of the problem is shown in Figure 1 contains a Newtonian fluid (air in this case) confined within enclosure delimited by two vertically eccentric hemispheres. Initially, the enclosure is at uniform temperature. From this time, the inner hemisphere of radius i R centered in i O is heated by a heat flux of constant density while the outer one of radius e R centered respectively on e O is maintened isothermal. The eccentricity e′ is defined as the algebric distance separating the two centers of hemispheres. The walls separating the two domes ( 0 θ = and θ = π ) are insulated. The enclosure is inclined with respect to the vertical by an angle ϕ from 0 to π. The two hemispheres being heated differently, a transient convection develops inside the enclosed space.
We considered the fluid Newtonian and the flow incompressible and bidimensional. Viscous dissipation, compressibility effects an thermal radiation are neglected. We further applied the Boussinesq approximation. An adapted conformal tranformation is used in order to reduce a curvilinear enclosure into a rectangular field. A bispherical coordinate system is chosen. The dimensionnal where a is the parameter of torus pole.
The inner copula is materialized by line of coordinate ( ) where Ψ and Ω are defined by: We defined the stream-function which intresically verifies incompressibility condition. His particularty is he has the dimension of a volume flow.
With φ is the stream-function which has a surface flow.
The dimensionless equation of stream-function is: These equations are subject to the following boundary conditions and initial: • On the inner hemisphere ( • On the outer hemisphere (

Numerical Procedure
Our physical domain is complex. the bispheric coordinate system we used allowed us to translate our curved walls into a parallelepiped grid. The implicit method alterning direction (A.D.I.) [15] is used for time integration with momentum and heat equations and the finite difference scheme is choosen for spatial integration. The system of linear equations obtained by using Thomas algorithm. The equations of stream-function is solved by using the Successive Over relaxation (S.O.R.) [16]. The convergence criterion as following is chosen for all variables dependant. It must be satisfied the following stop criterion: where m denotes the iteration and w stands for , ,T Ψ Ω .

Computation Conditions
Preliminary tests have been done on the influence on the mess and time step (see Table 1 and Table 2). We retain 51 × 51 for the grid system and 10 −5 for the time which constitutes a good compromise between a precision and an acceptable computation time.              tice that the value of the Nusselt number is almost the same for the angles π/3, π/2 and 2π/3. On the other hand, we note a difference for the value of π and 0, the value of Nusselt number for π is more important and its agrees with isotherms wich were more distorted than those of 0 ϕ = . The same remarks are observed for the variation of the average temperature of heated wall, namely that the values of the temperature for the angles π/3, π/2 and 2π/3 are equal and are the highest. The value of ϕ = π is the lower. The evolution of minimum of stream, it decrease with and the angle of inclinaison.    Likewise, the isotherms are much more deformed and more fill the cavity for this structure. When the eccentricity increases, we note a decrease in heat exchanges.

Conclusions
The numerical simulation of transient natural convection between two inclined hemispheres whose the inner hemisphere is subjected to a heat flux of constant density while the outer one is maintened isothermal with time is conducted. The bispherical coordinate and finite difference methods are used.
It emerges from the results presented that when the Rayleigh number is equal to 10 3 , heat transfers are dominated by conduction whatever the inclination and eccentricity. The convective transfers intensify as the Rayleigh number increases and there is a strongly destabilized regime for a Rayleigh number equal to 10 7 .
The simulations have also shown that the topology of the flow is strongly dependent on the inclinaison because the flow can change from a monocellular regime to a bicellular regime by varying the inclinaison. We remark also when the eccentricity increase, we note a decrease in heat exchanges. It would also be interesting to replace the Newtonian fluid with a complex fluid, for example, a fluid whose viscosity depends on the temperature.