Roughness effects in turbulent forced convection

We conducted direct numerical simulations (DNSs) of turbulent flow over three-dimensional sinusoidal roughness in a channel. A passive scalar is present in the flow with Prandtl number $Pr=0.7$, to study heat transfer by forced convection over this rough surface. The minimal channel is used to circumvent the high cost of simulating high Reynolds number flows, which enables a range of rough surfaces to be efficiently simulated. The near-wall temperature profile in the minimal channel agrees well with that of the conventional full-span channel, indicating it can be readily used for heat-transfer studies at a much reduced cost compared to conventional DNS. As the roughness Reynolds number, $k^+$, is increased, the Hama roughness function, $\Delta U^+$, increases in the transitionally rough regime before tending towards the fully rough asymptote of $\kappa_m^{-1}\log(k^+)+C$, where $C$ is a constant that depends on the particular roughness geometry and $\kappa_m\approx0.4$ is the von K\'arm\'an constant. In this fully rough regime, the skin-friction coefficient is constant with bulk Reynolds number, $Re_b$. Meanwhile, the temperature difference between smooth- and rough-wall flows, $\Delta\Theta^+$, appears to tend towards a constant value, $\Delta\Theta^+_{FR}$. This corresponds to the Stanton number (the temperature analogue of the skin-friction coefficient) monotonically decreasing with $Re_b$ in the fully rough regime. Using shifted logarithmic velocity and temperature profiles, the heat transfer law as described by the Stanton number in the fully rough regime can be derived once both the equivalent sand-grain roughness $k_s/k$ and the temperature difference $\Delta \Theta^+_{FR}$ are known. In meteorology, this corresponds to the ratio of momentum and heat transfer roughness lengths, $z_{0m}/z_{0h}$, being linearly proportional to $z_{0m}^+$, the momentum roughness length [continued]...