Leray-α simulations of wall-bounded turbulent flows

The Leray-α model reduces the range of active scales of the Navier-Stokes equations by smoothing the advective transport. Here we assess the potential of the Leray-α model in its standard formulation to simulate wall-bounded flows. Three flow cases are considered: plane channel flow at Reτ=590, Rayleigh-Bénard convection at Ra=107 and Pr = 1, and a side-heated vertical channel at Ra=5×106 and Pr = 0.7. The simulations are compared to results from a well resolved and coarse DNS. It is found that for all three flow cases, the variance in the velocity field increases as the filter width parameter a is increased, where a is connected to the filter width as αi=aΔxi, with Δxi the local grid size. Furthermore, the viscous and diffusive wall regions tend to thicken relative to the coarse DNS results as a function of a. In the cases where coarse DNS overpredicts wall gradients (as for Rayleigh-Bénard convection and the side-heated vertical channel), the thickening is beneficial. However, for the plane channel flow, coarse DNS underpredicts the wall-shear velocity, and increasing a only degrades the results. It is shown that buoyancy effects need to be included with care, because of the close relation of turbulent heat flux and the production of turbulent kinetic energy.