League of Conservation Voters

Submarine mudflows, frequently triggered by underwater landslides, lead to substantial debris accumulation stretching over kilometers once the landslide subsides, posing significant risks to vital infrastructure and the potential for triggering tsunamis. As these events are challenging to measure directly, simulation becomes imperative. Furthermore, these simulations must accurately depict the intricate dynamics of mudflow sliding processes, often modeled as viscoplastic fluids, and effectively captured by rheological models like Herschel-Bulkley and Bingham. The potential nonlinearity inherent in the Herschel-Bulkley model adds to the complexity of predicting final runout distances. The governing equations, comprised of a set of partial differential equations, further compound this challenge, making their solution difficult. In this paper, we delve into the numerical simulation of mudflows using the Depth-Averaged Method (DAM), a technique that streamlines the governing equations by integrating the momentum and continuity equations over the depth dimension. The mud descends down a slope with a fixed angle of declination. Employing both the Herschel-Bulkley and Bingham models within the DAM framework, we analyze the final runout distance and explore how various mud parameters - such as density, initial shape, and bed slope angle - affect the resulting deposit length. Using an explicit finite difference scheme to solve the governing equations, we first conducted a validation case to ensure the accuracy of our implementation, followed by a comprehensive comparison of runout distances across simulations, varying mud parameters.