Investigating landslide-induced tsunamis using a low-dissipation Riemann weakly compressible smoothed particle hydrodynamics model

Abstract

In order to address the complex moving boundary problem and violent free-surface flows involved in the process of landslide tsunamis, this paper proposes a new low-dissipation Riemann-solver-based smoothed particle hydrodynamics (SPH) numerical model. In this model, the complex fluid-rigid slide interaction and friction effects on landslides are handled by an improved motion boundary approach based on a one-sided Riemann scheme. A weighted essentially non-oscillatory (WENO) reconstruction method is introduced to reduce excessive dissipation of the fluid caused by rigid landslide impact. The kernel gradient accuracy is assured by applying the corrected smoothed particle method (CSPM) into a continuity equation. Single-block and multi-block landslide tsunamis in two-dimensional (2D) or three-dimensional (3D) applications are simulated with different particle resolutions, and the results are compared to those of other numerical models and experimental results. The landslide trajectories obtained from the modified boundary model are in good agreements with the experimental values. Free-surface flows variations during the propagation are properly captured in the current simulations without introducing any artificial viscosity term, and the wave peaks exhibit deviations of −14%–10.5% from the experimental values, thus verifying the usability and robustness of the proposed Riemann-SPH model.