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.
@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: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
volume = {37},
year = {2003},
pages = {389-416},
doi = {10.1051/m2an:2003034},
mrnumber = {1991208},
zbl = {1137.65392},
language = {en},
url = {http://dml.mathdoc.fr/item/M2AN_2003__37_2_389_0}
}
Audusse, Emmanuel; Bristeau, Marie-Odile. Transport of pollutant in shallow water : a two time steps kinetic method. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 37 (2003) pp. 389-416. doi : 10.1051/m2an:2003034. http://gdmltest.u-ga.fr/item/M2AN_2003__37_2_389_0/
[1] , and, 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.
[2] and, Upwind methods for hyperbolic conservation laws with source terms. Comput. & Fluids 23 (1994) 1049-1071. | Zbl 0816.76052
[3] and, Boundary conditions for the shallow water equations solved by kinetic schemes. INRIA Report, RR-4282 (2001), http://www.inria.fr/RRRT/RR-4282.html.
[4] and, Transport of pollutant in shallow water using kinetic schemes. CEMRACS, ESAIM Proc. 10 (1999) 9-21, http://www.emath.fr/Maths/Proc/Vol.10. | Zbl pre01614316
[5] , and, Finite volume methods, Handbook of numerical analysis, Vol. VIII, P.G. Ciarlet and J.L. Lions Eds., Amsterdam, North-Holland (2000). | MR 1804748 | Zbl 0981.65095
[6] , and, Some approximate Godunov schemes to compute shallow water equations with topography. Comput. & Fluids 32 (2003) 479-513. | Zbl 1084.76540
[7] and, Derivation of viscous Saint-Venant system for laminar shallow water; Numerical validation. Discrete Contin. Dynam. Systems 1 (2001) 89-102. | Zbl 0997.76023
[8] and, Numerical approximation of hyperbolic systems of conservation laws. Springer-Verlag, New York, Appl. Math. Sci. 118 (1996). | MR 1410987 | Zbl 0860.65075
[9] and, A well-balanced scheme designed for inhomogeneous scalar conservation laws. C. R. Acad. Sci. Paris Sér. I Math. 323 (1996) 543-546. | Zbl 0858.65091
[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 1001.35083
[12] , Numerical Methods for Conservation Laws. Second edition, ETH Zurich, Birkhauser, Lectures in Mathematics (1992). | MR 1153252 | Zbl 0723.65067
[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 0931.76059
[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 2064166 | Zbl 1030.35002
[16] and, A kinetic scheme for the Saint-Venant system with a source term. Calcolo 38 (2001) 201-231. | Zbl 1008.65066
[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 0626.65086
[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 03.0482.04
[19] , The formation of breakers and bores. Comput. Appl. Math. 1 (1948). | MR 24307 | Zbl 0041.54602