Numerical modeling of the propagation and morphological changes of turbidity currents using a cost-saving strategy of solution updating

Abstract Existing layer-averaged numerical models for turbidity currents have mostly adopted the global minimum time step (GMiTS) for solution updating, which confines their computational efficiency and limits their attractiveness for field applications. This paper presents a highly efficient layer-averaged numerical model for turbidity currents by implementing the combined approach of the local graded-time-step (LGTS) and the global maximum-time-step (GMaTS). The governing equations are solved for unstructured triangular meshes by the shock-capturing finite volume method along with a set of well-balanced evaluations of the numerical flux and geometrical slope source terms. The quantitative accuracy of the model, given reasonably estimated empirical and model parameters (e.g., bed friction, water entrainment, sediment deposition and erosion coefficients), is demonstrated by comparing the numerical solutions against laboratory data of the current front positions and deposition profiles, as well as field data of the current front positions. The improved computational efficiency is demonstrated by comparing the computational cost of the present model against that of a traditional model that uses a GMiTS. For the present simulated cases, the maximum reduction of the computational cost is approximately 80% (e.g., a simulation that cost 1 h before will only require 12 min with the new model).