“…As in [7,8], discretization of equations is carried out on a staggered grid to prevent mismatch between the velocity and pressure fields defined by simultaneous iterations and the artificial compressibility method. In the preliminary work stages, effects of the grid size and time step, the grid rarefaction/refinement, the length/width of the computational domain, and boundary conditions are studied.…”