“…The grid is uniform along the streamwise and spanwise directions. In the vertical direction, to increase the grid resolution near the free surface and the wall, the mesh is stretched with a transformation used in [23,25], so that the minimum vertical grid sizes about Dz þ min ¼ 0:1, 0.06 and 0.03 are used for Pr = 1, 10 and 100, respectively. According to Bergant and Tiselj [31], the grid system used is ensured to be capable of resolving all the essential scales of turbulence, with the smallest resolved temperature scales equal to the smallest scales of the velocity field for high-Pr calculations.…”