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 flat 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.
Mots clés : shallow water equations, Saint-Venant system, finite volume method, kinetic scheme, transport of pollutant, time discretization
@article{M2AN_2003__37_2_389_0, author = {Audusse, Emmanuel and Bristeau, Marie-Odile}, title = {Transport of pollutant in shallow water : a two time steps kinetic method}, journal = {ESAIM: Mod\'elisation math\'ematique et analyse num\'erique}, pages = {389--416}, publisher = {EDP-Sciences}, volume = {37}, number = {2}, year = {2003}, doi = {10.1051/m2an:2003034}, mrnumber = {1991208}, zbl = {1137.65392}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an:2003034/} }
TY - JOUR AU - Audusse, Emmanuel AU - Bristeau, Marie-Odile TI - Transport of pollutant in shallow water : a two time steps kinetic method JO - ESAIM: Modélisation mathématique et analyse numérique PY - 2003 SP - 389 EP - 416 VL - 37 IS - 2 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an:2003034/ DO - 10.1051/m2an:2003034 LA - en ID - M2AN_2003__37_2_389_0 ER -
%0 Journal Article %A Audusse, Emmanuel %A Bristeau, Marie-Odile %T Transport of pollutant in shallow water : a two time steps kinetic method %J ESAIM: Modélisation mathématique et analyse numérique %D 2003 %P 389-416 %V 37 %N 2 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an:2003034/ %R 10.1051/m2an:2003034 %G en %F M2AN_2003__37_2_389_0
Audusse, Emmanuel; Bristeau, Marie-Odile. Transport of pollutant in shallow water : a two time steps kinetic method. ESAIM: Modélisation mathématique et analyse numérique, Tome 37 (2003) no. 2, pp. 389-416. doi : 10.1051/m2an:2003034. http://www.numdam.org/articles/10.1051/m2an:2003034/
[1] Kinetic schemes for Saint-Venant equations with source terms on unstructured grids. INRIA Report, RR-3989 (2000), http://www.inria.fr/RRRT/RR-3989.html.
, and ,[2] Upwind methods for hyperbolic conservation laws with source terms. Comput. & Fluids 23 (1994) 1049-1071. | Zbl
and ,[3] Boundary conditions for the shallow water equations solved by kinetic schemes. INRIA Report, RR-4282 (2001), http://www.inria.fr/RRRT/RR-4282.html.
and ,[4] Transport of pollutant in shallow water using kinetic schemes. CEMRACS, ESAIM Proc. 10 (1999) 9-21, http://www.emath.fr/Maths/Proc/Vol.10.
and ,[5] Finite volume methods, Handbook of numerical analysis, Vol. VIII, P.G. Ciarlet and J.L. Lions Eds., Amsterdam, North-Holland (2000). | MR | Zbl
, and ,[6] Some approximate Godunov schemes to compute shallow water equations with topography. Comput. & Fluids 32 (2003) 479-513. | Zbl
, and ,[7] Derivation of viscous Saint-Venant system for laminar shallow water; Numerical validation. Discrete Contin. Dynam. Systems 1 (2001) 89-102. | Zbl
and ,[8] Numerical approximation of hyperbolic systems of conservation laws. Springer-Verlag, New York, Appl. Math. Sci. 118 (1996). | MR | Zbl
and ,[9] A well-balanced scheme designed for inhomogeneous scalar conservation laws. C. R. Acad. Sci. Paris Sér. I Math. 323 (1996) 543-546. | Zbl
and ,[10] Hydrodynamique des écoulements à surface libre, apport de la méthode des éléments finis. EDF (2001).
,[11] A steady state capturing method for hyperbolic system with geometrical source terms. ESAIM: M2AN 35 (2001) 631-646. | Numdam | Zbl
,[12] Numerical Methods for Conservation Laws. Second edition, ETH Zurich, Birkhauser, Lectures in Mathematics (1992). | MR | Zbl
,[13] Balancing source terms and flux gradients in high-resolution Godunov methods: the quasi-steady wave-propagation algorithm. J. Comput. Phys. 146 (1998) 346-365. | Zbl
,[14] Fonctionnement écologique de la Seine à l'aval de la station d'épuration d'Achères: données expérimentales et modélisation bidimensionnelle. Ph.D. Thesis, École des Mines de Paris, France (2001).
,[15] Kinetic formulations of conservation laws. Oxford University Press (2002). | MR | Zbl
,[16] A kinetic scheme for the Saint-Venant system with a source term. Calcolo 38 (2001) 201-231. | Zbl
and ,[17] Upwind differencing schemes for hyperbolic conservation laws with source terms, in Nonlinear Hyperbolic Problems, C. Carasso, P.A. Raviart and D. Serre Eds., Berlin, Springer-Verlag, Lecture Notes in Math. 1270 (1987) 41-51. | Zbl
,[18] Théorie du mouvement non permanent des eaux, avec application aux crues de rivières et à l'introduction des marées dans leur lit. C. R. Acad. Sci. Paris Sér. I Math. 73 (1871) 147-154. | JFM
,[19] The formation of breakers and bores. Comput. Appl. Math. 1 (1948). | MR | Zbl
,Cité par Sources :