skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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.

Authors:
 [1]
  1. Unitée de Mathématiques Pures et Appliquées, Ecole Normale Supérieure de Lyon, 46 allée d'Italie, 69364 Lyon Cedex 07 (France)
Publication Date:
OSTI Identifier:
22314861
Resource Type:
Journal Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 264; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0021-9991
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; DEFORMATION; FLOW MODELS; FLUIDS; LAGRANGIAN FUNCTION; PERFORMANCE; PLASTICITY; RESOLUTION; SHAPE; SOLIDS; STRESSES; SURFACES; VARIATIONAL METHODS; VISCOUS FLOW; WATER

Citation Formats

Fernández-Nieto, Enrique D., E-mail: edofer@us.es, Gallardo, José M., E-mail: jmgallardo@uma.es, and Vigneaux, Paul. Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case. United States: N. p., 2014. Web. doi:10.1016/J.JCP.2014.01.026.
Fernández-Nieto, Enrique D., E-mail: edofer@us.es, Gallardo, José M., E-mail: jmgallardo@uma.es, & Vigneaux, Paul. Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case. United States. https://doi.org/10.1016/J.JCP.2014.01.026
Fernández-Nieto, Enrique D., E-mail: edofer@us.es, Gallardo, José M., E-mail: jmgallardo@uma.es, and Vigneaux, Paul. 2014. "Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case". United States. https://doi.org/10.1016/J.JCP.2014.01.026.
@article{osti_22314861,
title = {Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case},
author = {Fernández-Nieto, Enrique D., E-mail: edofer@us.es and Gallardo, José M., E-mail: jmgallardo@uma.es and Vigneaux, Paul},
abstractNote = {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.},
doi = {10.1016/J.JCP.2014.01.026},
url = {https://www.osti.gov/biblio/22314861}, journal = {Journal of Computational Physics},
issn = {0021-9991},
number = ,
volume = 264,
place = {United States},
year = {Thu May 01 00:00:00 EDT 2014},
month = {Thu May 01 00:00:00 EDT 2014}
}