Simulation of Three-Dimensional Free-Surface Flows Using Two-Dimensional Multilayer Shallow Water Equations

We present an efficient and conservative Eulerian-Lagrangian method for solving two-dimensional hydrostatic multilayer shallow water flows with mass exchange between the vertical layers. The method consists of a projection finite volume method for the Eulerian stage and a method of characteristics to approximate the numerical fluxes for the Lagrangian stage. The proposed method is simple to implement, satisfies the conservation property and it can be used for multilayer shallow water equations on non-flat bathymetry including eddy viscosity and Coriolis forces. It offers a novel method of calculating stratified vertical velocities without the use of the Navier-Stokes equations. Numerical results are presented for several examples and the obtained results for a free-surface flow problem are in close agreement with the analytical solutions. We also test the performance of the proposed method for a test example of wind-driven flows with recirculation.