Asymmetric secondary flows above geometrically symmetric surface roughness

The mean flow behaviour of a turbulent boundary layer over rough walls is expected to exhibit symmetries that govern the flow dynamics. In particular, when roughness elements are arranged in a spanwise symmetric manner, the mean flow above them should also exhibit spanwise symmetry. This symmetrical consideration has garnered substantial empirical support. We conduct direct numerical simulations (DNS) of flow over aligned cube arrays to test such symmetry considerations further. We vary the surface coverage density from 0.25% to 6.25%, and employ an averaging time of about 100 large-eddy turnover times, which is longer than the typical averaging time in prior DNS studies of rough-wall boundary layers. The results suggest the presence of spanwise asymmetry in the mean flow. Specifically, we observe the development of a prominent secondary vortex on one side of the cubical roughness, accompanied by a relatively smaller secondary vortex on the other side. This asymmetry becomes most pronounced when the surface coverage density is approximately 0.59%, and diminishes as the coverage density approaches either a low or a high value. We also establish that this mean flow asymmetry is robust across variations in the domain size, the initial condition, and the placement of the cubes in the spanwise direction.