View article

Many flows especially in geophysics involve turbulent boundary layers forming over rough surfaces with multiscale height distribution. Such surfaces pose special challenges for large-eddy simulation (LES) when the filter scale is such that only part of the roughness elements of the surface can be resolved. Here we consider LES of flows over rough surfaces with power-law height spectra Eh(k) ~ kβs (−3 ≤ βs < −1), as often encountered in natural terrains. The surface is decomposed into resolved and subgrid-scale height contributions. The effects of the unresolved small-scale height fluctuations are modelled using a local equilibrium wall model (log-law or Monin–Obukhov similarity), but the required hydrodynamic roughness length must be specified. It is expressed as the product of the subgrid-scale root-mean-square of the height distribution and an unknown dimensionless quantity, α, the roughness parameter …