Non-hydrostatic 3D free surface layer-structured finite volume model for short wave propagation

In this paper a layer-structured finite volume model for non-hydrostatic 3D environmental free surface flow is presented and applied to several test cases, which involve the computation of gravity waves. The 3D unsteady momentum and mass conservation equations are solved in a collocated grid made of polyhedrons, which are built from a 2D horizontal unstructured mesh, by just adding several horizontal layers. The mesh built in such a way is unstructured in the horizontal plane, but structured in the vertical direction. This procedure simplifies the mesh generation and at the same time it produces a well-oriented mesh for stratified flows, which are common in environmental problems. The model reduces to a 2D depth-averaged shallow water model when one single layer is defined in the mesh. Pressure-velocity coupling is achieved by the Semi-Implicit Method for Pressure-Linked Equations algorithm, using Rhie-Chow interpolation to stabilize the pressure field. An attractive property of the model proposed is the ability to compute the propagation of short waves with a rather coarse vertical discretization. Several test cases are solved in order to show the capabilities and numerical stability of the model, including a rectangular free oscillating basin, a radially symmetric wave, short wave propagation over a ID bar, solitary wave runup on a vertical wall, and short wave refraction over a 2D shoal. In all the cases the numerical results are compared either with analytical or with experimental data. Copyright (C) 2008 John Wiley & Sons, Ltd.