Two-dimensional simulation by regularization of free surface viscoplastic flows with Drucker-Prager yield stress and application to granular collapse

Christelle Lusso; Alexandre Ern; François Bouchut; Anne Mangeney; Maxime Farin; Olivier Roche

Two-dimensional simulation by regularization of free surface viscoplastic flows with Drucker-Prager yield stress and application to granular collapse

Christelle Lusso ¹ Alexandre Ern ^{1, 2} François Bouchut ³ Anne Mangeney ^{4, 5} Maxime Farin ⁴ Olivier Roche ⁶

1 CERMICS - Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique

2 SERENA - Simulation for the Environment: Reliable and Efficient Numerical Algorithms

Inria de Paris

3 LAMA - Laboratoire d'Analyse et de Mathématiques Appliquées

4 IPGP - Institut de Physique du Globe de Paris

5 ANGE - Numerical Analysis, Geophysics and Ecology

LJLL - Laboratoire Jacques-Louis Lions, Inria Paris-Rocquencourt

6 LMV - Laboratoire Magmas et Volcans

Abstract : This work is devoted to numerical modeling and simulation of granular flows relevant to geophysical flows such as avalanches and debris flows. We consider an incompressible viscoplastic fluid, described by a rheology with pressure-dependent yield stress, in a 2D setting with a free surface. We implement a regularization method to deal with the singularity of the rheological law, using a mixed finite element approximation of the momentum and incompressibility equations, and an arbitrary Lagrangian Eulerian (ALE) formulation for the displacement of the domain. The free surface is evolved by taking care of its deposition onto the bottom and of preventing it from folding over itself. Several tests are performed to assess the efficiency of our method. The first test is dedicated to verify its accuracy and cost on a one-dimensional simple shear plug flow. On this configuration we setup rules for the choice of the numerical parameters. The second test aims to compare the results of our numerical method to those predicted by an augmented Lagrangian formulation in the case of the collapse and spreading of a granular column over a horizontal rigid bed. Finally we show the reliability of our method by comparing numerical predictions to data from experiments of granular collapse of both trapezoidal and rectangular columns over horizontal rigid or erodible granular bed made of the same material. We compare the evolution of the free surface, the velocity profiles, and the static-flowing interface. The results show the ability of our method to deal numerically with the front behavior of granular collapses over an erodible bed.

Keywords : Drucker–Prager pressure-dependent yield stress Viscoplastic flows granular collapse ALE method regularization free surface

Document type :

Journal articles

Domain :

Mathematics [math] / Numerical Analysis [math.NA]

Sciences of the Universe [physics] / Earth Sciences / Geophysics [physics.geo-ph]

Complete list of metadatas

Cited literature [60 references] Display Hide Download

https://hal-upec-upem.archives-ouvertes.fr/hal-01133786
Contributor : François Bouchut <>
Submitted on : Friday, December 16, 2016 - 3:18:38 PM
Last modification on : Wednesday, September 4, 2019 - 1:52:04 PM

drucker-prager2D.pdf

Files produced by the author(s)

Two-dimensional simulation by regularization of free surface viscoplastic flows with Drucker-Prager yield stress and application to granular collapse

File

Identifiers

Collections

Citation

Export

Metrics

1536

412

Contact

Two-dimensional simulation by regularization of free surface viscoplastic flows with Drucker-Prager yield stress and application to granular collapse

File

Identifiers

Collections

Citation

Export

Share

Metrics

1536

412

Contact