The wave diffraction around a submerged circular cylinder was simulated numerically with two dimensional waves. A special wave maker motion mode was determined by analysis in the frequency domain and simulation in the time domain to generate Stokes waves with second order accuracy. Numerical simulation employed the time stepping boundary integral equation method. Particular consideration is paid to the compatibility of the free surface and the body surface conditions at the intersection point, the ope...