In this paper, a modification of PISO algorithm based on standard k-ε turbulence model was proposed. The numerical technique used in this research is finite volume method, hybrid scheme for discretizing convection term, Euler implicit scheme for discretizing time term, and non-staggered grid. The basic idea of the modification of PISO algorithm is to perform an additional corrector stage for turbulence kinetic energy and dissipation rate to correct the inconsistence of flow and turbulence. In order to validate this algorithm, simulation of flow around a square cylinder (Re=3000) was performed in two-dimensional case. The results obtained from the proposed scheme show better agreement with those from the experiment than using original PISO algorithm in coherent velocity field.