We propose a three-dimensional non-hydrostatic shock-capturing numerical model for the simulation of wave propagation, transformation and breaking, which is based on an original integral formulation of the contravariant Navier-Stokes equations, devoid of Christoffel symbols, in general time-dependent curvilinear coordinates. A coordinate transformation maps the time-varying irregular physical domain that reproduces the complex geometries of coastal regions to a fixed uniform computational one. The advancing of the solution is performed by a second-order accurate strong stability preserving Runge-Kutta fractional-step method in which, at every stage of the method, a predictor velocity field is obtained by the shock-capturing scheme and a corrector velocity field is added to the previous one, to produce a non-hydrostatic divergence-free velocity field and update the water depth. The corrector velocity field is obtained by numerically solving a Poisson equation, expressed in integral contravariant form, by a multigrid technique which uses a four-colour Zebra Gauss-Seidel line-by-line method as smoother. Several test cases are used to verify the dispersion and shockcapturing properties of the proposed model in time-dependent curvilinear grids.

In this paper a new finite-volume non-hydrostatic and shock-capturing three-dimensional model for the simulation of wave-structure interaction and hydrodynamic phenomena (wave refraction, diffraction, shoaling and breaking) is proposed. The model is based on an integral formulation of the Navier-Stokes equations which are solved on a time dependent coordinate system: a coordinate transformation maps the varying coordinates in the physical domain to a uniform transformed space. The equations of motion are discretized by means of a finite-volume shock-capturing numerical procedure based on high order WENO reconstructions. The solution procedure for the equations of motion uses a third order accurate Runge-Kutta (SSPRK) fractional-step method and applies a pressure corrector formulation in order to obtain a divergence-free velocity field at each stage. The proposed model is validated against several benchmark test cases

In this paper a non-hydrostatic and shock-capturing model for the simulation of wave-structure interaction and hydrodynamic phenomena (wave refraction, diffraction, shoaling and breaking) is proposed. The model is based on an integral formulation of the motion equations which are solved on a time dependent curvilinear coordinate system: a coordinate transformation maps the varying coordinates in the physical domain to a uniform transformed space. The motion equations are discretized by means of a shock-capturing numerical procedure based on high order WENO reconstructions. The solution procedure for the motion equations uses a five stage four order accurate Runge-Kutta (SSPRK) fractionalstep method and applies a pressure corrector formulation in order to obtain a divergence-free velocity field at each stage. The proposed model is validated against two benchmark test cases and is used in order to analyse the effects that a submerged breakwater produces on wave motion.

Submerged shore-parallel breakwaters for coastal defence are a good compromise between the need to mitigate the effects of waves on the coast and the ambition to ensure the preservation of the landscape and water quality. In this work we simulate, in a fully three-dimensional form, the hydrodynamic effects induced by submerged breakwaters on incident wave trains with different wave height. The proposed three-dimensional non-hydrostatic finite-volume model is based on an integral form of the Navier-Stokes equations in σ-coordinates and is able to simulate the shocks in the numerical solution related to the wave breaking. The obtained numerical results show that the hydrodynamic phenomena produced by wave-structure interaction have features of three-dimensionality (undertow), that are locally important, and emphasize the need to use a non-hydrostatic fully-three-dimensional approach.