EGU General Assembly 2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.

A two-layer shallow flow model with two axes of integration with application to submarine avalanches and generated tsunamis

Enrique D. Fernandez-Nieto1, François Bouchut2, Juan M. Delgado-Sanchez1, Anne Mangeney3, and Gladys Narbona-Reina1
Enrique D. Fernandez-Nieto et al.
  • 1Universidad de Sevilla, Departamento de Matemática Aplicada I, Sevilla, Spain
  • 2Université Paris-Est, Laboratoire d'Analyse et de Mathématiques Appliquées (UMR 8050), CNRS, UPEM, UPEC, F-77454, Marne-la-Vallée, France
  • 3Institut de Physique du Globe de Paris, Seismology team, University Paris-Diderot, Sorbonne Paris Cité, 75238, Paris, France

There exits in the literature many approaches that has been used to model submarine avalanches (See [5]). These models are mainly based on the pioneer work of Savage and Hutter (SH) [4] that is a shallow water type model for aerial avalanches, which is written in local coordinates, in order to simulate the tangential velocity to the bottom. A depth-averaged SH model over a general bottom with curvature was introduced in [1]. An extension to submarine avalanches is developed in [2]. In this paper the same local coordinate system is used for the two layers. Nevertheless, using a local coordinates the model would prescribe the perturbation at the surface at a wrong placement. In [3] a bilayer depth-averaged model for submarine avalanches is presented with cartesian coordinates for the water layer and local coordinates for the avalanche. The drawback is that the seabed deformation is considered as an input data for the water layer equations, then no interaction between the two fluids are taken into account and it is necessary to do an interpolation of the granular surface at each time step of the numerical simulation. In this work we present firstly the details of the proposed model, a coupled two-layer shallow water system where we consider local coordinates for the granular layer and cartesian coordinates for the fluid one. The main difference with other models that adopt the same stragie is that any interpolation of the granular surface is required. Moreover, the velocity of the granular layer has an explicit influence on the mass and momentum conservation laws of the fluid layer. Secondly, several numerical tests will be presented.


[1] F. Bouchut, E.D. Fernández-Nieto, A. Mangeney, and P.Y. Lagrée. On new erosion models of Savage-Hutter type for avalanches. Acta Mechanica, 199(1):181--208, 2008.
[2] E.D. Fernández-Nieto, F. Bouchut, D. Bresch, M.J. Castro Díaz, and A. Mangeney. A new Savage-Hutter type model for submarine avalanches and generated tsunami. Journal of Computational Physics, 227(16):7720--7754, 2008.
[3] P.H. Heinrich, A. Piatanesi, and H. Hébert. Numerical modelling of tsunami generation and propagation from submarine slumps: the 1998 papua new guinea event. Geophysical Journal International, 145(1):97--111, 2001.
[4] S. B. Savage and K. Hutter. The dynamics of avalanches of granular materials from initiation to runout. part I: Analysis. Acta Mechanica, 86(1):201–223, 1991.
[5] S. Yavari-Ramshe and B. Ataie-Ashtiani. Numerical modeling of subaerial and submarine landslide-generated tsunami waves-recent advances and future challenges. Landslides, 13(6):1325–1368, 2016.

How to cite: Fernandez-Nieto, E. D., Bouchut, F., Delgado-Sanchez, J. M., Mangeney, A., and Narbona-Reina, G.: A two-layer shallow flow model with two axes of integration with application to submarine avalanches and generated tsunamis, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9485,, 2021.

Display materials

Display file