Numerical Simulation of Wave Propagation, Breaking, and Setup on Steep Fringing Reefs

The prediction of wave transformation and associated hydrodynamics is essential in the design and construction of reef top structures on fringing reefs. To simulate the transformation process with better accuracy and time efficiency, a shock-capturing numerical model based on the extended Boussinesq equations suitable for rapidly varying topography with respect to wave transformation, breaking and runup, is established. A hybrid finite volume–finite difference scheme is used to discretize conservation form of the extended Boussinesq equations. The finite-volume method with a HLL Riemann solver is applied to the flux terms, while finite-difference discretization is applied to the remaining terms. The fourth-order MUSCL (Monotone Upstream-centered Schemes for Conservation Laws) scheme is employed to create interface variables, with in which the van-Leer limiter is adopted to improve computational accuracy on complex topography. Taking advantage of van-Leer limiter, a nested model is used to take account of both computational run time and accuracy. A modified eddy viscosity model is applied to better accommodate wave breaking on steep reef slopes. The established model is validated with laboratory measurements of regular and irregular wave transformation and breaking on steep fringing reefs. Results show the model can provide satisfactory predictions of wave height, mean water level and the generation of higher harmonics