Skip to Main content Skip to Navigation
Journal articles

Efficient numerical schemes for viscoplastic avalanches. Part 1: the 1D case

Abstract : This paper deals with the numerical resolution of a shallow water viscoplastic flow model. Viscoplastic materials are characterized by the existence of a yield stress: below a certain critical threshold in the imposed stress, there is no deformation and the material behaves like a rigid solid, but when that yield value is exceeded, the material flows like a fluid. In the context of avalanches, it means that after going down a slope, the material can stop and its free surface has a non-trivial shape, as opposed to the case of water (Newtonian fluid). The model involves variational inequalities associated with the yield threshold: finite-volume schemes are used together with duality methods (namely Augmented Lagrangian and Bermúdez-Moreno) to discretize the problem. To be able to accurately simulate the stopping behavior of the avalanche, new schemes need to be designed, involving the classical notion of well-balancing. In the present context, it needs to be extended to take into account the viscoplastic nature of the material as well as general bottoms with wet/dry fronts which are encountered in geophysical geometries. We derived such schemes and numerical experiments are presented to show their performances.
Complete list of metadata

Cited literature [40 references]  Display  Hide  Download
Contributor : Paul Vigneaux <>
Submitted on : Monday, October 9, 2017 - 5:39:43 PM
Last modification on : Tuesday, December 8, 2020 - 10:11:58 AM
Long-term archiving on: : Wednesday, January 10, 2018 - 3:27:48 PM


Files produced by the author(s)




Enrique D. Fernandez-Nieto, J.M. Gallardo, Paul Vigneaux. Efficient numerical schemes for viscoplastic avalanches. Part 1: the 1D case. Journal of Computational Physics, Elsevier, 2014, 264, pp.55-90. ⟨10.1016/⟩. ⟨hal-01011373⟩



Record views


Files downloads