Computation of turbulent flow in general domains