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

TL;DR

This paper introduces a structure-preserving dual-field discretization method for 2D turbidity currents that conserves key physical quantities and accurately models dynamics with fewer degrees of freedom.

Contribution

The work presents a novel dual-field, structure-preserving discretization for turbidity currents that maintains energy, mass, enstrophy, and vorticity conservation, improving efficiency and accuracy.

Findings

Preserves energy balance with bounded residual.

Eliminates the need for solving fully nonlinear systems.

Achieves accurate dynamics with fewer degrees of freedom.

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.