Adaptation of f-wave finite volume methods to the two-layer shallow-water equations in a moving vessel with a rigid-lid

Abstract

A numerical method is proposed to solve the two-layer inviscid, incompressible and immiscible 1D shallow-water equations in a moving vessel with a rigid-lid with different boundary conditions based on the high-resolution f-wave finite volume methods due to Bale et al. (2002). The method splits the jump in the fluxes and source terms including the pressure gradient at the rigid-lid into waves propagating away from each grid cell interface. For the influx–efflux boundary conditions the time dependent source terms are handled via a fractional step approach. In the linear case the numerical solutions are validated by comparison with the exact analytical solutions. Numerical solutions presented for the nonlinear case include shallow-water sloshing waves due to prescribed surge motion of the vessel.