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

Slah Sahmim (1) , Fayssal Benkhaldoun (1) , Francisco Alcrudo (2)

1 LAGA - Laboratoire Analyse, Géométrie et Applications
2 CPS - CPS

This 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.

Mots clés: Finite volumes; Riemann Problems; SRNH scheme; Source terms; Shallow Water equations; Euler equations; Two-phase flow

Domaines: Analyse numérique [math.NA]

Slah Sahmim, Fayssal Benkhaldoun, Francisco Alcrudo. A sign matrix based scheme for non-homogeneous PDE's with an analysis of the convergence stagnation phenomenon. Journal of Computational Physics, 2007, 226 (2), pp.1753-1783. ⟨10.1016/j.jcp.2007.06.017⟩. ⟨hal-00114522⟩