An unstructured finite-volume method for coupled models of suspended sediment and bed load transport in shallow-water flows

The aim of this work is to develop a well-balanced finite-volume method for the accurate numerical solution of the equations governing suspended sediment and bed load transport in two-dimensional shallow-water flows. The modelling system consists of three coupled model components: (i) the shallow-water equations for the hydrodynamical model; (ii) a transport equation for the dispersion of suspended sediments; and (iii) an Exner equation for the morphodynamics. These coupled models form a hyperbolic system of conservation laws with source terms. The proposed finite-volume method consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The gradient fluxes are discretized using a modified Roe's scheme using the sign of the Jacobian matrix in the coupled system. A well-balanced discretization is used for the treatment of source terms. In this paper, we also employ an adaptive procedure in the finite-volume method by monitoring the concentration of suspended sediments in the computational domain during its transport process. The method uses unstructured meshes and incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep sediment concentrations and bed load gradients that may form in the approximate solutions. Details are given on the implementation of the method, and numerical results are presented for two idealized test cases, which demonstrate the accuracy and robustness of the method and its applicability in predicting dam-break flows over erodible sediment beds. The method is also applied to a sediment transport problem in the Nador lagoon

Schéma SRNHS Analyse et Application d'un schéma aux volumes finis dédié aux systèmes non homogènes

International audienceThis article is devoted to the analysis, and improvement of a finite volume scheme proposed recently for a class of non homogeneous systems. We consider those for which the corressponding Riemann problem admits a selfsimilar solution. Some important examples of such problems are Shallow Water problems with irregular topography and two phase flows. The stability analysis of the considered scheme, in the homogeneous scalar case, leads to a new formulation which has a naturel extension to non homogeneous systems. Comparative numerical experiments for Shallow Water equations with sourec term, and a two phase problem (Ransom faucet) are presented to validate the scheme.Cet article concerne l'analyse et l'application, d'un schéma proposé récemment por une classe de systèmes non homogènes. Nous considérons ceux pour lesquels le problème de Riemann correpondant admet une solution autosimilaire. Deux exemples importants de tels problèmes sont l'écoulement d'eau peu profonde au-dessus d'un fond non plat et les problèmes diphasiques. l'analyse de stabilité du schéma, dans le cas scalaire homogène, amène à une nouvelle écriture qui a une extension naturelle pour le cas non homogène. Des expériences numériques comparatives pour des équations de saint-Venant avec topographie variable, et un problème diphasique (Robinet de Ransom) sont présentés pour évaluer l'efficacité du schéma

A non homogeneous Riemann Solver for shallow water and two phase flows

In this work we consider a two steps finite volume scheme, recently developed to solve nonhomogeneous systems. The first step of the scheme depends on a diffusion control parameter which we modulate, using the limiters theory. Results on Shallow water equations and two phase flows are presented

Finite volume characteristic flux scheme for transonic flow problems

This work deals with the numerical solution of internal transonic flow problems. Currently we expect 1D and 2D inviscid flow of perfect gas modelled by the Euler equations. We use finite volume characteristic flux scheme, called VFFC scheme, which can be viewed as a generalization of Roe scheme. The dissipation matrix is computed analytically or numerically. We present numerical results for 1D shock tube problem computed by the second order method, where the spatial accuracy is due to linear reconstruction with the minmod limiter and temporal discretization is done using explicit three stage Runge-Kutta method. Further numerical results for 2D transonic flow in GAMM channel have been achieved by the first order method on structured quadrilateral as well as unstructured triangular meshes

A sign matrix based scheme for non-homogeneous PDE's with an analysis of the convergence stagnation phenomenon

International audienceThis work is devoted to the analysis of a finite volume method recently proposed for the numerical computation of a class of non homogenous systems of partial differencial equations of interest in fluid dynamics. The stability analysis of the proposed scheme leads to the introduction of the sign matrix of the flux jacobian. It appears that this formulation is equivalent to the VFRoe scheme introduced in the homogeneous case and has a natural extension here to non homogeneous sys- tems. Comparative numerical experiments for the Shallow Water and Euler equa- tions with source terms, and a model problem of two phase flow (Ransom faucet) are presented to validate the scheme. The numerical results present a convergence stagnation phenomenon for certain forms of the source term, notably when it is singular. Convergence stagnation has been also shown in the past for other numerical schemes. This issue is addressed in a specific section where an explanation is given with the help of a linear model equation, and a cure is demonstrated

Finite volume methods for unidirectional dispersive wave models

We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular we consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves and their various interactions.Comment: 25 pages, 12 figures, 51 references. Other authors papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh

Tsunami generation by dynamic displacement of sea bed due to dip-slip faulting

In classical tsunami-generation techniques, one neglects the dynamic sea bed displacement resulting from fracturing of a seismic fault. The present study takes into account these dynamic effects. Earth's crust is assumed to be a Kelvin-Voigt material. The seismic source is assumed to be a dislocation in a viscoelastic medium. The fluid motion is described by the classical nonlinear shallow water equations (NSWE) with time-dependent bathymetry. The viscoelastodynamic equations are solved by a finite-element method and the NSWE by a finite-volume scheme. A comparison between static and dynamic tsunami-generation approaches is performed. The results of the numerical computations show differences between the two approaches and the dynamic effects could explain the complicated shapes of tsunami wave trains.Comment: 16 pages, 10 figures, Accepted to Mathematics and Computers in Simulation. Other author's papers can be downloaded at http://www.cmla.ens-cachan.fr/~dutyk