This paper incorperates Bingham and bi-viscosity rheology models with the Navier–Stokes solver to simulate the dynamics and kinematics processes of slumps for tsunami generation. The rheology models are integrated into a computational fluid dynamics code, Splash3D, to solve the incompressible Navier–Stokes equations with volume of fluid surface tracking algorithm. The change between un-yield and yield phases of the slide material is controlled by the yield stress and yield strain rate in Bingham and bi-viscosity models, respectively. The integrated model is carefully validated by the theoretical results and laboratory data with good agreements. This validated model is then used to simulate the benchmark problem of the failure of the gypsum tailings dam in East Texas in 1966. The accuracy of predicted flood distances simulated by both models is about 73% of the observation data. To improve the prediction, a fixed large viscosity is introduced to describe the un-yield behavior of tailings material. The yield strain rate is obtained by comparing the simulated inundation boundary to the field data. This modified bi-viscosity model improves not only the accuracy of the spreading distance to about 97% but also the accuracy of the spreading width. The un-yield region in the modified bi-viscosity model is sturdier than that described in the Bingham model. However, once the tailing material yields, the material returns to the Bingham property. This model can be used to simulate landslide tsunamis.

Energy dissipation mechamism is the key to study tsunami hazard mitigation. Numerical method is adopted to study the interaction between bores and square cylinders. The model solves the three-dimensional Navier–Stokes equations with Large-Eddy Simulation turbulence model. The Volume-of-fluid (VOF) method is used to track the complex free surface. We focus the investigation on the effect of cylinder height on the flow field. The results show that the turbulence diffusion is the main mechanism for energy dissipation. The flow patterns are significantly different within and beyond the cylinder array. The taller cylinders cause smaller velocity magnitude in the downstream area. In addition, a larger value of velocity magnitude and vorticity near the bottom is identified in the tall-cylinder case. These unique featuers make different dissipation rates.

This paper incorporates the Bingham rheology model with the Navier–Stokes solver to simulate the tsunamis excited by a slump-type landslide. The slump is modeled as the Bingham material, in which the rheological properties changing from the un-yield phase to yield phase is taken into account. The volume of fluid method is used to track the interfaces between three materials: air, water, and slump. The developed model is validated by the laboratory data of the benchmark landslide tsunami problem. A series of rheological properties analyses is performed to identify the parameter sensitivity to the tsunami generation. The results show that the yield stress plays a more important role than the yield viscosity in terms of the slump kinematics and tsunami generation. Moreover, the scale effect is investigated under the criterion of Froude number similarity and Bingham number similarity. With the same Froude number and Bingham number, the result from the laboratory scale can be applied to the field scale. If the slump material collected in the field is used in the laboratory experiments, only the result of the maximum wave height can be used, and significant errors in slump shape and moving speed are expected.