CMP (Chemical Mechanical Polishing) is to achieve adequate local and global planarization for future sub-micrometer VLSI requirements. In designing CMP, numerical computation is quite helpful in terms of reducing the amount of experimental works. Stresses on pad, concentration of particles and particle tracking are studied for design. In this research, the optimization of grooved pad shape of CMP is performed through numerical investigation of slurry flow in CMP process. The result indicates that the combination of sinusoidal groove and skewed pad is the most optimal shape among the twenty candidates. Useful information can be obtained in velocity, pressure, stress, concentration of particles and particles trajectories, etc.