Classiﬁcation of Flow Modes for Natural Convection in a Square Enclosure with an Eccentric Circular Cylinder

: The present study investigated the natural convection for a hot circular cylinder embedded in a cold square enclosure. The numerical simulations are performed to solve a two-dimensional steady natural convection for three Rayleigh numbers of 10 3 , 10 4 and 10 5 at a ﬁxed Prandtl number of 0.7. This study considered the wide range of the inner cylinder positions to identify the eccentric effect of the cylinder on ﬂow and thermal structures. The present study classiﬁes the ﬂow structures according to the cylinder position. Finally, the present study provides the map for the ﬂow structures at each Rayleigh number ( Ra ). The Ra = 10 3 and 10 4 form the four modes of the ﬂow structures. These modes are classiﬁed by mainly the large circulation and inner vortices. When Ra = 10 5 , one mode that existed at Ra = 10 3 and 10 4 , disappears in the map of the ﬂow structures. The new three modes appear, resulting in total six modes of ﬂow structures at Ra = 10 5 . New modes at Ra = 10 5 are characterized by the top side secondary vortices. The corresponding isotherms are presented to explain the bifurcation of the ﬂow structure.

Shu and Zhu [19] numerically investigated the natural convection in a concentric annulus between a heated inner cylinder and a cold square enclosure. They found that the patterns of flow and thermal fields is significantly affected by the aspect ratio and the Rayleigh number (Ra).
Kim et al. [26] considered the variation of the vertical position of the inner hot cylinder in 10 3 ≤ Ra ≤ 10 6 . They reported that the dependence of the inner vortices and the primary circulation on the vertical position. Later, Yoon et al. [27] extended the re-search of Kim et al. [26] by considering the higher Ra of 10 7 , showing the occurrence of time dependence.
Hussain and Hussein [28] considered a uniform heat flux as a boundary condition of the inner cylinder alike the problem with Kim et al. [26]. They reported that the effect of a uniform heat flux condition on the Nusselt number in comparison with Kim et al. [26].
Lee et al. [36] considered the change of the cylinder location along the diagonal and horizontal directions in 10 3 ≤ Ra ≤ 10 6 . They revealed that the variation of the mean and local Nusselt numbers accompanying with the modification of the inner cells and the primary circulation according to the cylinder position.
As reviewed above, the natural convection for the internal cylinder in the cavity has been well established by many researches, resulting in greatly improving the understanding of this heat transfer, but based on the authors' literature survey, the previous studies considered the limited positions of the cylinder for the eccentric effect on the natural convection in the enclosure. Mainly three different directions of the vertical, horizontal, and diagonal direction for the position of the cylinder are dealt with the investigation of the eccentric effect.
Therefore, the main emphasis of the present study is to cover the wide range of the internal cylinder positions to identify the eccentric effect on the thermal and flow structures in the enclosure. The present study classifies the flow structures in the enclosure for the inner cylinder location. Finally, the present study provides the maps for the flow structures to reinforce the understanding of the eccentricity of the cylinder within the enclosure for three Rayleigh numbers. The isotherms are presented to explain the bifurcation of the flow structure. The distributions of the Nusselt numbers for the enclosure are presented as the contours in the polar plane embedded in the enclosure in order to match with the map of the flow structures according to the positions of the inner cylinder center. Figure 1 presents a schematic of the present system with the side length (L) of the square enclosure and the cylinder radius of R = 0.2L. In order to identify the inner cylinder location, the polar coordinates (r, θ) are introduced within the enclosure. The origin of the polar coordinate system is the enclosure center. The radial and circumferential positions of the center of the inner cylinder are in the range of 0 ≤ r ≤ 0.25L with intervals of ∆r = 0.05L and −90 • ≤ θ ≤ 90 • with intervals of ∆θ = 10 • , respectively. To resolve the diagonal positions, θ = 45 • and θ = −45 • are also considered. The solutions of the present problems are steady and symmetric about x = 0. Thus, we consider the cylinder positions in the right half of the enclosure. In total 10 6 positions of the inner cylinder are considered in this study.

Numerical Details
The constant low temperature (T c ) and high temperature (T h ) are put on the enclosure walls and the cylinder surface, respectively. The present study does not consider the radiation which is assumed to be negligible. The Boussinesq approximation is used for the density. The constant assumption is used for the rest of properties. Thus, under the present assumptions, the following continuity, momentum and energy equations satisfying the unsteady incompressible viscous flow and the thermal fields are the dimensionless governing equations: where, the superscript * represents the dimensionless variable. In the above equations, u i are the components of velocity corresponding to the Cartesian coordinates of x i , p is the pressure, t is the time and T is the temperature. The dimensionless variables in the above equations are defined as: where, α and ρ denote the thermal diffusivity and the density, respectively. The non-dimensionalization results in dimensionless parameters of Rayleigh number, and Prandtl number, Pr = ν/α, where β, g and ν are the expansion coefficient, the gravitational acceleration and the kinematic viscosity, respectively. The value of Pr is 0.7 which corresponds to the air. The present study considers three values of Ra of 10 3 , 10 4 and 10 5 .
For the discretization of the space, a second-order central difference scheme is utilized to structured grids, based on the finite volume method (FVM). A second-order implicit scheme is employed for temporal discretization. The Crank-Nicolson scheme and the second-order Adams-Bashforth are adopted for the diffusion and convection terms, respectively. A two-step, time-split scheme is used for the time advance of the flow and thermal fields. The solution of Poisson equation satisfies the continuity equation. Further details of this scheme are can be referred from Kim and Moin [48] and Zang et al. [49]. We derived the governing equations using a generalized curvilinear coordinate system.
The boundary conditions of the no-penetration and no-slip are imposed on the walls for the flow. As the wall boundaries for the energy equation, the hot and cold isothermal conditions are defined for the internal cylinder and the enclosure, respectively.
After achieving the velocity and temperature, the local and surface-averaged Nusselt numbers are obtained by following definitions; where n and A are the normal direction to the wall and the surface area. Figure 2 shows the present grid system. The present computations use the grid numbers of 201 × 201 in the horizontal (x) and vertical (y) directions. The non-uniform grid distribution is generated by using the algebraic function to resolve the high gradients near the walls. Particularly, the dense resolution of the grid is maintained by tuning the grid points to cover the variation of the cylinder position. The test of the grid dependence is performed by using coarser and much finer grids of 102(x) × 102(y), 303(x) × 303(y) and 402(x) × 402(y). The difference of Nusselt numbers between grid numbers becomes smaller, which is less than 0.5%, 0.1% and 0.07% for the coarse, medium and fine grids. Finally, the grid resolution of 201(x) × 201 (y) is used in all commutations for the computational efficiency. For the validation of the present numerical methods, the surface-averaged Nusselt numbers of the inner cylinder for different values of Ra are compared with those of Park et al. [18], Kim et al. [26] and Lee et al. [36], as presented in Table 1. The present results agree well with the previous researches. In addition, in order to identify the reliability of the present numerical methods to achieve the solutions of the eccentricity between the enclosure and the cylinder, the previous results [26,40] for total surface-averaged Nu of the square enclosure (Nu en ) and the inner circular cylinder (Nu cy ) for the vertical variation are considered. The present results agree well with the previous results [26,40], as presented in Figure 3.

Map and Patterns of Flow and Thermal Fields
The map is introduced to show the patterns of thermal and flow fields in the (r, θ) plane which covers the inner cylinder positions. In order to reinforce the understanding of the eccentricity of the cylinder within the enclosure, the (r, θ) plane for the map is inserted in the enclosure, as shown in Figure 4. In general, the Mode I covers the circular formation with the enclosure center and extends to the upper and lower narrow areas long the vertical center line. At the enclosure center of r = 0 and θ = 0 • , the isothermal lines are almost linearly distributed along the radial direction. The corresponding streamlines reveal two inner vortices embedded in two symmetric circulations, as presented in Figure 5a. Therefore, Mode I reveals two overall rotating primary eddies and each rotating eddy has two inner vortices (TPE_TIV).
When the internal cylinder positioned below the horizontal center of y = 0 and moves to lower boundary of Mode I, the area between the bottom wall and the cylinder becomes narrow. Thus, the size of lower inner vortex is reduced gradually, but the upper inner vortex is getting bigger. The corresponding isotherms is getting denser in the region between the bottom wall and the inner cylinder, as shown in Figure   In Mode I, in order to investigate the effect of the circumferential position on the thermal and flow fields, the variation of the circumferential direction is considered from the vertical center line (θ = 90 • and θ = −90 • ) to the horizontal center line (θ = −0 • ). When θ approaches counterclockwise from θ = −90 • to θ = 0 • within Mode I, the two inner vortices in the left and right primary counter-clockwise eddies become smaller and larger, respectively, as shown in Figure 5c for θ = −30 • . Consequently, when θ approaches counterclockwise from θ = −90 • to θ = 0 • , the distance between the right wall and inner cylinder enclosure becomes shorter. Thus, the gradient of isotherms augments, as the primary eddy becomes smaller in the right side. When θ approaches clockwise from θ = 90 • to θ = 0 • θ = 0 • within Mode I, the thermal and flow fields undergo the symmetric pattern about the horizontal center line at y = 0 to those for the counterclockwise rotation from θ = −90 • to θ = 0 • .
When the inner cylinder gets away from the boundary of Mode I for (r = 0.1L, θ = 90 • ) and move downward and upward, two inner vortices are merged into the primary eddy at (r = 0.15L, θ = 90 • ) in Figure 6a and (r = 0.15L, θ = −90 • ) which is symmetry to (r = 0.15L, θ = 90 • ). Therefore, Mode II reveals two overall rotating primary eddies (TPE). The area of Mode II become wider in the circumferential direction with increasing the redial position of the inner cylinder center. The critical circumferential locations for the radial positions of r = 0.2L and r = 0.25L are θ = 80 • and θ = 70 • , respectively, as shown in Figure 4 for the Map. Thus, as the radial position increases, the critical circumferential location, where the flow structure changes, decreases. The effect of the cylinder radial position is that the core of the primary circulation moves toward the same direction of the variation of the internal cylinder location because the enclosure can secure spaces opposite the inner cylinder moving enough to enlarge the core of the circulation. When r increases within Mode II, isotherms become denser in between the cylinder and the enclosure wall facing to the cylinder, otherwise isotherms are getting sparser in the opposite region.
In Mode II, the effect of the circumferential position on the flow and thermal fields is almost the same as Mode I. When the inner cylinder rotates clockwise from θ = 90 • to θ = 70 • within Mode II, the space in the direction of the inner cylinder rotating becomes small, as shown in Figure 6b. As a result, in this space, the primary circulations are squeezed and small. The isotherms form dense and coarse gradients in the cylinder side and opposite regions, respectively. The counter-clockwise rotation from θ = −90 • to θ = −70 • within Mode II undergoes the same variation as the clockwise rotation from θ = 90 • to θ = 70 • .
When the inner cylinder gets away from the boundary of Mode I for (r = 0.1L, θ = 90 • ) and rotates clockwise from θ = 90 • to θ = 80 • , in the left side of the enclosure, two secondary inner eddies at (r = 0.1L, θ = 90 • ) is merged into the primary eddy at (r = 0.1L, θ = 80 • ), as shown in Figure 7a. When θ continuously turn clockwise from θ = 80 • to θ = 0 • , the center of the merged inner vortex in the left circulation rotates clockwise, as shown in Figure 7a,b for θ = 80 • to θ = 80 • , respectively. Simultaneously, in the right primary eddy, the upper inner vortex becomes bigger and the lower inner vortex squeezes into the clockwise direction.
Therefore, Mode III reveals two overall rotating primary eddies which have two inner vortices in the primary circulation in the direction of the inner cylinder rotation and one inner vortex in the opposite primary circulation (TPE_TOIV). When the inner cylinder gets away from the boundary of Mode II holding TPE, the primary circulation in the inner cylinder rotating direction contains additional inner vortex. When the inner cylinder circumferentially approaches to the horizontal center (θ = 0 • ), this additional inner vortex becomes bigger and the lower inner vortex squeezes into the clockwise direction. When θ approaches counter-clockwise from θ = −90 • with the negative sign in the lower half of the enclosure to θ = 0 • within Mode III, the flow and thermal fields undergo the symmetric pattern about the horizontal center line to those for the clockwise rotation from θ = 90 • to θ = 0 • . When the position of the inner cylinder changes along the radial direction, gets away from the boundary of Mode II and enters into Mode IV, the primary eddy in the direction of the cylinder moving are divided in Figure 8 for r = 0.1L and θ = 0 • . Also, when the inner cylinder rotates along the circumferential direction, gets away from the boundary of Mode III and enters into Mode IV, two inner vortices in the primary eddy are divided in the direction of the inner cylinder rotation. Therefore, Mode IV reveals one overall rotating primary eddy and two divided vortices (OPE_TDV). The isothermal lines and streamlines at Ra = 10 4 have almost the same formations as those at Ra = 10 4 . For each mode, the effects of the radial and circumferential locations of the internal cylinder center on the thermal and flow fields at Ra = 10 4 are very similar to those at Ra = 10 4 . Thus, the map of the flow structures for Ra = 10 4 consists with the same four modes of Mode I, Mode II, Mode II and Mode IV corresponding to TPE_TIV, TPE, TPE_TOIV, and OPE_TDV, respectively, as Ra = 10 3 . However, the buoyance in the direction of the gravity gives a greater effect on both fields, because of the augmentation of the convection. Thus, in contrast to Ra = 10 3 having a four-fold symmetry, the map for Ra = 10 4 is a two-fold symmetry x = 0, as presented in Figure 9.
In general, the Ra = 10 4 case forms thinner thermal boundary layer on the cylinder bottom than the upper side, and a smaller and weaker inner lower vortex than the upper one, comparing to Ra = 10 3 . The regions of Modes I, II and IV for Ra = 10 4 are larger than the regions of corresponding Modes for Ra = 10 3 , as shown in Figures 4 and 9, respectively.  As the cylinder is positioned downward from the enclosure at r = 0 along the vertical center line at θ = −90 • , two internal vortices merges at r = 0.1L. Namely, Mode I changes into Mode II, which is earlier than r = 0.15L at Ra = 10 3 . In Mode II, as r increases at θ = −90 • , isotherms becomes distorted and the thermal plumes on the cylinder top is getting bigger.
As the cylinder is located upward at θ = 90 • , the occurrence of the bifurcation from Mode I to Mode II at r = 0.15L is later than that for the downward movement of the cylinder, due to the stronger convective flow in the space between the top wall and the cylinder.
In Mode II formation, when the cylinder is positioned further upward, the two overall rotating primary eddies are more obligated than the downward movement of the cylinder, as observed in Figure 10b  At r = 0.2L and r = 0.25L, the flow structures undergo three modes by the effect of the circumferential eccentricity of the cylinder. At r = 0.2L, when the cylinder rotates clockwise from θ = 90 • to θ = 0 • , Mode II changes into Mode III at θ = 70 • and then the bifurcation from Mode III to Mode IV occurs at θ = 0 • , as presented in Figure 12a-c, respectively. When the internal cylinder rotates counter-clockwise from θ = −90 • to θ = 0 • , Mode II changes into Mode III at θ = −70 • and then the bifurcation from Mode III to Mode IV occurs at θ = −20 • , as presented in Figure 12d-f, respectively. At r = 0.25L, when the inner cylinder rotates, the bifurcation from Mode II to Mode III occurs later than that at r = 0.2L, resulting in that the size of Mode II at r = 0.25L is wider than that at r = 0.2L, as observed in Figure 9 for the Map. In comparison with Ra = 10 3 , the region for Mode IV enlarges when Ra = 10 4 .

Ra = 10 5
When Ra = 10 5 , the Mode I with TPE_TIV, existed at Ra = 10 4 and 10 4 , disappears in the map of the flow structures, as shown in Figure 13. The new three modes of V, VI and VII appear. The Mode I is divided into about the four modes of II, III, V, VI and VII.
The center of the enclosure is occupied by Mode II, since a stronger buoyancy makes merging the lower inner vortex into the upper vortex of Mode I which is existed at Ra = 10 3 and 10 4 . Therefore, the centers of the overall rotating eddies locate at the upper half of the enclosure and the upward thermal plume is originated from the internal cylinder in Figure 14a.
The radial eccentricity by the cylinder downward movement along the vertical centerline from (r = 0, θ = −90 • ), the Mode II keep maintains, as shown Figure 14b for r = 0.25L, respectively. As the cylinder moves downward from the enclosure center, the circulation is getting larger in the upper half of the enclosure and the upward thermal plume is more elongated, since wider area between the bottom wall and the cylinder secures a stronger convection in the upper half of the enclosure.   Figure 13. When the cylinder moves upward from the enclosure center, the mode changes from Mode II to Mode V and then recovers again Mode II, which is clearly identified in Figure 14c,d for r = 0.05L and r = 0.15L, respectively.
In the range of 0.05L ≤ r ≤ 0.1L corresponding to Mode V in the map, a fair of secondary vortices appears in the cylinder upper surface. Therefore, Mode V presents two overall rotating primary eddies and secondary two vortices (TPE_STV). The corresponding isothermal lines exhibit that one upward plume in Mode II is divided into two plumes from the cylinder and one downward plume from the enclosure top wall, respectively, as presented in Figure 14c.
As the cylinder moves further upward from r = 0.1L to r = 0.25L, the area between the enclosure top wall and cylinder is not enough to activate the secondary two vortices in this Ra = 10 5 . In other words, the buoyance force is not enough to overcome the conduction and to keep the convection to maintain the secondary two vortices in this narrow space at this Ra = 10 5 . Thus, the secondary two vortices disappear, as shown in Figure 14d for r = 0.15L. Mainly, the isotherms are linearly distributed in this narrow region between the enclosure top wall and cylinder. As a result, Mode V changes into Mode II.
The circumferential eccentricity by the clockwise rotation of the cylinder from (r = 0.05L, θ = 90 • ) to the horizontal center line (r = 0.05L, θ = 0 • ) leads to two bifurcations from Mode V to Mode VI and from Mode VI to Mode III, as shown in Figure 13. When θ = 70 • at r = 0.05L, the upper secondary vortices still appear, but these secondary vortices become asymmetric, since the left upper secondary vortex is squeezed and smaller by the augmentation of the left primary eddy due to the enlargement of the left upper region by the clockwise rotation of the cylinder. Particularly, the right primary eddy contains additional tiny lower inner vortex by the compression of the clockwise rotation of the cylinder, as presented in Figure 15a. As a result, Mode VI forms secondary two vortices and two overall rotating primary eddies which have two inner vortices in a large circulation in the direction of the inner cylinder rotation and one inner vortex in the opposite primary circulation. Finally, Mode VI is named as TPE_TOIV_STV. When the inner cylinder clockwise rotates from θ = 70 • to θ = 50 • at r = 0.05L, the upper secondary vortices becomes smaller but the lower inner vortex in the right primary eddy becomes bigger at θ = 60 • in Figure 15b, and then the upper secondary vortices disappears at θ = 50 • in Figure 15c. Therefore, at θ = 50 • , the bifurcation from Mode VI to Mode III occurs. The Mode III maintains up to θ = 0 • . At r = 0.1L, the circumferential eccentricity by the cylinder clockwise rotating from θ = 90 • to the horizontal center line of θ = 0 • leads to three bifurcations from Mode V to Mode VI, from Mode VI to Mode VII, and from Mode VII to Mode III, as shown in Figure 13 for the Map.
At r = 0.1L, when the inner cylinder clockwise rotates from θ = 90 • to θ = 60 • , first bifurcation from Mode V to Mode VI occurs, as shown in Figure 16a. When θ = 45 • , the left upper secondary vortex disappears and simultaneously the right upper secondary vortex becomes smaller, as shown in Figure 16b. Thus, new flow structure forms in Mode VII where one secondary upper vortex in company with TPE_TOIV. The Mode VII is defined as TPE_TOIV_SOV. When the inner cylinder reaches to θ = 40 • , the right upper secondary vortex also disappears, as shown in Figure 16c, resulting in Mode III which maintains to θ = 0 • . At r = 0.15L, the clockwise rotation of the inner cylinder from θ = 90 • to θ = 0 • derives three bifurcations from Mode II to Mode III, from Mode III to Mode VII, and from Mode VII to Mode III, as shown in Figure 14 for the Map. At this radial location, Mode II at θ = 90 • changes into Mode III at θ = 80 • , and the lower inner vortex of the primary eddy in the rotational direction becomes bigger at θ = 70 • , as shown in Figure 17a for θ = 80 • .
When θ = 60 • at r = 0.15L, the secondary tiny vortex appears on the cylinder upper surface and it exists up to θ = 45 • under the main flow structure of Mode III, as shown in Figure 17b for θ = 60 • . The primary circulation in the opposite direction to the rotation of the cylinder strongly turns around the inner cylinder surface in the opposite side circulation, as the inner cylinder rotates clockwise from θ = 90 • .
When θ = 40 • at r = 0.15L, the secondary tiny vortex disappears by the augmentation in size and strength of the upper inner vortex of the primary eddy in the opposite direction to the rotation of the inner cylinder, resulting in Mode III which is preserved up to θ = 0 • , as shown in Figure 17c at θ = 40 • . Therefore, at r = 0.15L, the third bifurcation from Mode VII to Mode III occurs at θ = 40 • . At r = 0.2L, when the inner cylinder rotates clockwise from θ = 90 • to θ = 0 • and counterclockwise from θ = −90 • to θ = 0 • , single bifurcation from Mode II to Mode III occurs at θ = 80 • and θ = −60 • , respectively, as shown in Figure 13 for the Map.
At r = 0.25L, when the inner cylinder rotates clockwise from θ = 90 • to θ = 0 • , two bifurcations occur, as shown in Figure 13. First bifurcation from Mode II to Mode III occurs at θ = 70 • , as shown in Figure 18a for θ = 70 • . The second bifurcation from Mode III to Mode IV occurs at θ = 20 • and Mode IV preserves up to θ = 0 • , as shown in Figure 18b for θ = 20 • . At r = 0.25L, the counter-clockwise rotation from θ = −90 • to θ = 0 • provides two bifurcations, as shown in Figure 13. First bifurcation from Mode II to Mode III occurs at θ = −60 • , as shown in Figure 19a. The second bifurcation from Mode III to Mode IV occurs at θ = −30 • , as shown in Figure 19b. The circumferential eccentricity by the counter-clockwise rotation from θ = −90 • to θ = 0 • in the lower half leads to almost the same distribution of isotherms, regardless of the radial position. At θ = −90 • , the isotherms form upward and large plume appears on the cylinder, which forms sharper thermal gradient near the enclosure top. The dense isotherms form over the lower surface of the cylinder due to the increased buoyance with increasing Rayleigh number. The stagnant area under the cylinder becomes wider with augmenting r, resulting in very coarse isotherms. When the inner cylinder rotates counterclockwise from θ = −90 • to θ = 0 • , the upward thermal plume inclines to the cylinder rotational direction.

Nusselt Number
The distributions of the surface-averaged Nusselt numbers ( Nu ) for each walls of the enclosure are presented as the contours in the (r, θ) plane in the (r, θ) plane embedded in the enclosure in order to match with the map of the flow structures according to the positions of the cylinder center, as presented in Figure 20.
For Ra = 10 3 , the gap between the surface of the enclosure wall and the cylinder is predominant to the change of Nu , as shown in Figure 20a-c for the top, bottom, right walls of the enclosure, respectively. The Nu for the bottom and top walls are dependent on nearly the vertical direction, as shown in Figure 20a,b, respectively, due to the strong conduction at low Ra = 10 3 . Because the upward movement of the cylinder forms denser and coarser isotherms near the top and bottom walls, respectively, the larger and smaller values of Nu for the top and bottom walls appear at θ = 90 • and θ = −90 • for r = 0.25L. The cylinder downward movement derives the opposite patterns of Nu for the bottom and top walls to the upward movement.
At this low Ra = 10 3 , the right wall exhibits the symmetric distribution of Nu about the center line of x = 0, observed in Figure 20c. The values of Nu for the right wall is dependent on nearly the horizontal direction. Therefore, the distribution of Nu is composed of almost the vertical lines. However, the careful observation identifies the very weak curvature of the contours which means the existence of the weak convection at low Ra = 10 3 . The magnitude of Nu for the right wall is governed by the distance between the enclosure wall and the cylinder surface, like that for the bottom and top walls.
For Ra = 10 4 , the distributions of Nu for the enclosure walls are almost the same as Ra = 10 3 , which is consistent with the similar maps for the thermal and flow structures for Ra = 10 4 and Ra = 10 3 . Therefore, the Nu for the bottom and top walls depend on mainly the vertical distance and form the horizontal lines based contours, as presented in Figure 20d,e. The Nu for the right wall is dependent on nearly the horizontal direction. Thus, the distribution of Nu is composed of almost the vertical lines, as shown in Figure 20f. However, the increase of the convection with augmenting Ra from 10 3 and Ra = 10 4 contributes the contours of Nu to be more curvature. As Ra becomes large to Ra = 10 5 , the distributions of Nu for the walls of the enclosure considerably change in comparison with those for lower Ra = 10 3 and Ra = 10 4 , which can be clarified in Figure 20g-i,a-c and d-f for Ra = 10 5 , Ra = 10 3 and Ra = 10 4 , respectively.
As discussed in the thermal and flow structures, when Ra = 10 5 , the large upward thermal plume form the cylinder upper surface characterizes these fields due to the increase of the buoyance. Therefore, the magnitudes of Nu for the top wall becomes much bigger than those for lower Ra = 10 3 and Ra = 10 4 .
As the cylinder moves upward along the centerline from (r = 0.25L, θ = −90 • ) which is the nearest cylinder position to the bottom wall, the Nu for the top wall increases, decreases and increases again. The increase of Nu according the upward movement of the cylinder from (r = 0.25L, θ = −90 • ) is attributed to the combined effects of the augmentation of the conduction due to the reduction gap between the cylinder and the top wall and a stronger convection generating the large upward thermal plume. As the cylinder keeps moving upward, the value of Nu decreases by the slightly weaken convection due to a smaller area between the cylinder and the top wall. Also, the downward plume from the top wall of the enclosure, corresponding to the secondary vortices over the upper surface of the inner cylinder, contributes to the decrease of Nu for the top wall. Finally, as the internal cylinder continuously moves upward, the Nu increases and reaches to the maximum.
When the cylinder becomes far away, the upward thermal plume inclines and forms coarse isotherms near the top wall. Therefore, this inclined thermal plume contributes the value of Nu for the top wall to the reduction. The small Nu distribute around the outer boundary of the position of the inner cylinder, as presented in Figure 20g.
As the cylinder locates in the lower half of the enclosure for Ra = 10 5 , the distribution of Nu for the bottom wall is similar to that for lower Ra = 10 3 and Ra = 10 4 . Namely, when inner cylinder moves upward from (r = 0.25L, θ = −90 • ) to the enclosure center (r = 0, θ = 0 • ), the Nu decreases, as shown in Figure 20h. Otherwise, when the inner cylinder locates in the upper half for Ra = 10 5 , the variation of Nu is very weak with a small value, since the thermal stratification appears in the area under the cylinder, as early observed in the isotherms.
As the cylinder locates in the right half of the enclosure for Ra = 10 5 , the Nu for the right wall depends on mainly the gap between the right wall and cylinder, resulting in the same increasing pattern of Nu with closer to the right wall as lower Ra = 10 3 and Ra = 10 4 , as observed in Figure 20i. Otherwise, as the cylinder locates in the left half of the enclosure, the coarse isotherms form the opposite corners. Hence, the small Nu s for the right wall appear as the cylinder positions near the left bound of the position of the cylinder. Figure 21 shows the total surfaces-averaged Nusselt number of enclosure, Nu en , in the present ranges of the radial and circumferential locations of the inner cylinder center for different Ra. The Ra = 10 3 and 10 4 form about the same distribution of Nu en , which is almost the four-fold symmetry about x = 0 and y = 0, as presented in Figure 21a,b, respectively. As discussed in the thermal and flow fields, because the gap is predominant to the heat transfer, Nu en has large values at the horizontal and vertical boundaries of the cylinder position. In other words, the magnitude of Nu en is mainly as a function of the radial direction at these low Ra of 10 3 and 10 4 . For the validation in Figure 3, Nu en has a parabolic profile according to r, regardless of θ.
As Ra increases to 10 5 , the location of the minimum of Nu en moves diagonally upward from the center of the enclosure (x = 0, y = 0) where the minimum of Nu en occurs at Ra = 10 3 and 10 4 in Figure 21c. In the regions showing the minimum of Nu en , the inclined rising thermal plumes on the upper surface of the cylinder form the secondary vortices which accompany the down welling plume from the top wall of the enclosure. This down welling plume derives the coarse isotherms near the enclosure top wall, which contributes to forming regions containing the minimum of Nu en . The value of Nu en near the bottom are bigger than that those near the top wall, because the stagnant region at the top wall becomes wider than that near the bottom.

Conclusions
The present study researched the natural convection in a cold square enclosure with a hot circular cylinder. The numerical simulations are performed to solve a 2-D unsteady natural convection for Ra of Ra = 10 3 , 10 4 and 10 5 at Pr = 0.7.
This study covered the wide range of the cylinder positions to identify the eccentric effect of the inner cylinder on flow and thermal structures in the square enclosure. The present study classified the flow structures according to the cylinder center in the enclosure. The isotherms were presented to explain the bifurcation of the flow structure.
The maps for the flow structures were provided to reinforce the understanding of the eccentricity of the cylinder within the enclosure for different Rayleigh numbers. The Rayleigh numbers of 10 3 and 10 4 form the four modes of the flow structure. These modes are classified by mainly the large circulation and inner vortices. Therefore, the maps for these Rayleigh numbers consist of four modes.
Mode I reveals two overall rotating primary eddies and each rotating eddy has two inner vortices (TPE_TIV). When the inner cylinder gets away from the boundary of Mode I, the two inner vortices are merged. Mode II reveals two overall rotating primary eddies (TPE). When the inner cylinder gets away from the boundary of Mode II, the primary circulation in the inner cylinder rotating direction contains additional inner vortex. Therefore, Mode III reveals two overall rotating primary eddies which have two inner vortices in the primary circulation in the direction of the inner cylinder rotation and one inner vortex in the opposite primary circulation (TPE_TOIV). When the inner cylinder gets away from the boundary of Mode III, two inner eddies embedded in the primary eddy in the direction of the inner cylinder rotation are divided. Therefore, Mode IV reveals one overall rotating primary eddy and two divided vortices (OPE_TDV).
When Ra = 10 5 , the Mode I with TPE_TIV, existed at Ra = 10 3 and 10 4 , disappears in the map of the flow structures. The new three modes of V, VI and VII appear. Mode V presents two overall rotating primary eddies and secondary two vortices (TPE_STV). Mode VI forms secondary two vortices and two overall rotating primary eddies which have two inner vortices in the primary circulation in the direction of the inner cylinder rotation and one inner vortex in the opposite primary circulation. Finally, Mode VI is named as TPE_TOIV_STV. Mode VII has one secondary upper vortex in company with TPE_TOIV. The Mode VII is defined as TPE_TOIV_SOV. At Ra = 10 5 , these new modes are characterized by the top side secondary vortices.