Transport of pollutant in shallow water a two time steps kinetic method

被引:37
作者
Audusse, E [1 ]
Bristeau, MO [1 ]
机构
[1] Inst Natl Rech Informat & Automat, F-78153 Le Chesnay, France
来源
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE | 2003年 / 37卷 / 02期
关键词
shallow water equations; Saint-Venant system; finite volume method; kinetic scheme; transport of pollutant; time discretization;
D O I
10.1051/m2an:2003034
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The aim of this paper is to present a finite volume kinetic method to compute the transport of a passive pollutant by a flow modeled by the shallow water equations using a new time discretization that allows large time steps for the pollutant computation. For the hydrodynamic part the kinetic solver ensures - even in the case of a non. at bottom - the preservation of the steady state of a lake at rest, the non-negativity of the water height and the existence of an entropy inequality. On an other hand the transport computation ensures the conservation of pollutant mass, a non-negativity property and a maximum principle for the concentration of pollutant and the preservation of discrete steady states associated with the lake at rest equilibrium. The interest of the developed method is to preserve these theoretical properties with a scheme that allows to disconnect the hydrodynamic time step - related to a classical CFL condition - and the transport one - related to a new CFL condition - and further the hydrodynamic calculation and the transport one. The CPU time is very reduced and we can easily solve different transport problems with the same hydrodynamic solution without large storage. Moreover the numerical results exhibit a better accuracy than with a classical method especially when using 1D or 2D regular grids.
引用
收藏
页码:389 / 416
页数:28
相关论文
共 19 条
  • [1] AUDUSSE E, 2000, RR3989 INRIA
  • [2] UPWIND METHODS FOR HYPERBOLIC CONSERVATION-LAWS WITH SOURCE TERMS
    BERMUDEZ, A
    VAZQUEZ, E
    [J]. COMPUTERS & FLUIDS, 1994, 23 (08) : 1049 - 1071
  • [3] Bristeau M.-O., 2001, 4248 RECH
  • [4] BRISTEAU MO, 1999, ESAIM P, V10, P9
  • [5] EYMARD R, 2000, HDB NUMERICAL ANAL, V8
  • [6] Gerbeau JF, 2001, DISCRETE CONT DYN-B, V1, P89
  • [7] GODLEWSKI E, 1996, APPL MATH SCI, V118
  • [8] Gosse L, 1996, CR ACAD SCI I-MATH, V323, P543
  • [9] Hervouet J. M., 2001, HYDRODYNAMIQUE ECOUL
  • [10] A steady-state capturing method for hyperbolic systems with geometrical source terms
    Jin, S
    [J]. ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 2001, 35 (04): : 631 - 645