For time advancement, a second order accurate implicit scheme is employed, which uses an algebraic multi-grid method and a constant time step of 1 107 s. A second order bounded central-difference scheme is used for spatial discretization to guarantee boundedness and high numerical accuracy [19]. The normalizedvariable diagram value is computed based on the local flow variables to determine the amount of upwind-blending to stabilize the central difference scheme when required [20]. Fureby and Grinstein [21] show that this kind of discretization schemes are suitable for implicit LES calculations.