The computation of incompressible turbulent flow with two-equation closure models (k-ε and k-ω) is considered. The Cartesian staggered grid approach is generalized to general boundary-fitted coordinates. An accurate discretization on nonsmooth grids is presented. For higher-order monotone discretization of the equations for the turbulence quantities, flux-limited versions of the k-scheme are developed. In order to better assess the relative merits of explicit and implicit time discretization, a new approach to obtain von Neumann stability conditions is presented. A comparison is made between physical time scales for direct and large-eddy simulation, and stability restrictions on the time step for explicit schemes. Applications are presented for stationary turbulent flow computations with the k-ω model.
