A two-layer shallow flow model with two axes of integration, well-balanced discretization and application to submarine avalanches

J.M. Delgado-Sánchez (1) , François Bouchut (2) , E.D. Fernández-Nieto (1) , Anne Mangeney (3) , G. Narbona-Reina (1)

1 IMUS - Departamento de Matemática Aplicada I
2 LAMA - Laboratoire Analyse et de Mathématiques Appliquées
3 IPGP (UMR_7154) - Institut de Physique du Globe de Paris

We propose a two-layer model with two different axes of integration and a well-balanced finite volume method. The purpose is to study submarine avalanches and generated tsunamis by a depth-averaged model with different averaged directions for the fluid and the granular layers. Two-layer shallow depth-averaged models usually consider either Cartesian or local coordinates for both layers. However, the motion characteristics of the granular layer and the water wave are different: the granular flow velocity is mainly oriented downslope while water motion related to tsunami wave propagation is mostly horizontal. As a result, the shallow approximation and depth-averaging have to be imposed (i) in the direction normal to the topography for the granular flow and (ii) in the vertical direction for the water layer. To deal with this problem, we define a reference plane related to topography variations and use the associated local coordinates to derive the granular layer equations whereas Cartesian coordinates are used for the fluid layer. Depth-averaging is done orthogonally to that reference plane for the granular layer equations and in the vertical direction for the fluid layer equations. Then, a finite volume method is defined based on an extension of the hydrostatic reconstruction. The proposed method is exactly well-balanced for two kind of stationary solutions: the classical one, when both water and granular masses are at rest; the second one, when only the granular mass is at rest. Several tests are presented to get insight into the sensitivity of the granular flow, deposit and generated water waves to the choice of the coordinate systems. Our results show that even for moderate slopes (up to 30 •), strong relative errors on the avalanche dynamics and deposit (up to 60%) and on the generated water waves (up to 120%) are made when using Cartesian coordinates for both layers instead of an appropriate local coordinate system as proposed here.

Mots clés: Submarine avalanches; tsunami propagation; bilayer shallow water model; hydrostatic reconstruction; well-balanced finite volume method

Domaines: Géophysique [physics.geo-ph]; Analyse numérique [math.NA]

J.M. Delgado-Sánchez, François Bouchut, E.D. Fernández-Nieto, Anne Mangeney, G. Narbona-Reina. A two-layer shallow flow model with two axes of integration, well-balanced discretization and application to submarine avalanches. Journal of Computational Physics, 2020, 406, pp.109186. ⟨10.1016/j.jcp.2019.109186⟩. ⟨hal-02518084⟩