Abstract
A numerical method for the analysis of the wave breaking process capable of handling complex sea-beds and a free-surface was developed on the basis of the time-dependent, incompressible Navier-Stokes equations. The numerical method used the Boundary-Fitted Coordinate transformation (BFC) method consisting of the Boundary-Fitted Coordinates and Numerical Grid Generation techniques. Smoothing and relocation techniques were used to enhance accuracy and stability of the calculations. The accuracy of the present method was examined by performing calculations for the sloshing motions and plunging breaker on a uniform slope. A comparison of both the numerical method and the experimental results from LDV and image processing immediately before the breaker confirmed the validity of the numerical method.