The breaking phenomenon of regular periodic waves generated by a numerical wave maker is simulated by finite-difference method which can cope with strong interface motions. The air and water flows are simultaneously solved in the time-marching solution procedure for the Navier-Stokes equation. A density-function technique is devised for the implemenation of the interface conditions. The accuracy is examined and applied to the simulation of two-dimensional breaking phenomena of periodic gravity waves.