[zrpanah]

I tried to model natural convection in a square enclosure with four circular cylinders with ANSYS-FLUENT and compare the results with previous studies. Four inner cylinders placed equal distance away from each other within an enclosure moving along the diagonal center lines. The walls of the square enclosure were maintained at a constant low temperature of Tc, whereas the cylinders was maintained at a constant high temperature of Th. for the chosen Ra number the results are expected to be symmetric with respect to the vertical and horizontal center lines at x = 0 and y = 0, respectively.
The issue is:
1. Simulation results does show symmetric temperature contours with symmetric values but, although the stream functions have the same symmetry shape they don't have symmetric values and the stream function values at the right side of the enclosure are smaller than that of the left side. (Picture of the results is attached).
2. Distribution of time averaged local Nusselt number along the surface of the cylinders does not match the expected pattern.(picture is attached)

I would appreciate any suggestions on the reason why these issues are occurred and how to solve them.