A conservative, physically compatible, dual-field discretization for turbidity currents: with application to the lock-exchange problem

Abstract: In this work we present a structure preserving discretization for turbidity currents based on a mass-, energy-, enstrophy-, and vorticity-conserving formulation for 2D incompressible flows. This discretization exploits a dual-field formulation for the time evolution of the velocity and vorticity fields together with the transport equation for the particles. Due to its staggered time integration the resulting system of equations is quasi-linear, eliminating the need to solve for a fully nonlinear system of equations. It is shown that this discretization preserves the energy balance equation up to a bounded residual due to the staggering in time of the velocity and vorticity. This leads to a numerical scheme that does not introduce artificial energy dissipation. A comparison with literature results is presented showing that this approach can retrieve the dynamics of the system with a much smaller number of degrees of freedom.