I released a new set of boundary conditions for simulating rough surfaces in OpenFOAM, using fine meshes with near wall distances around y+=1 or lower.
kAupoixLowReWallFunctionkLeeLowReWallFunctionomegaAupoixWallFunctionomegaKnoppWallFunctionomegaLeeWallFunction
GitHub Repository, with source code and binaries: