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

Emmanuel Audusse; Marie-Odile Bristeau

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2003)

  • Volume: 37, Issue: 2, page 389-416
  • ISSN: 0764-583X

Abstract

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

How to cite

top

Audusse, Emmanuel, and 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 37.2 (2003): 389-416. <http://eudml.org/doc/245541>.

@article{Audusse2003,
abstract = {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.},
author = {Audusse, Emmanuel, Bristeau, Marie-Odile},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {shallow water equations; Saint-Venant system; finite volume method; kinetic scheme; transport of pollutant; time discretization},
language = {eng},
number = {2},
pages = {389-416},
publisher = {EDP-Sciences},
title = {Transport of pollutant in shallow water : a two time steps kinetic method},
url = {http://eudml.org/doc/245541},
volume = {37},
year = {2003},
}

TY - JOUR
AU - Audusse, Emmanuel
AU - Bristeau, Marie-Odile
TI - Transport of pollutant in shallow water : a two time steps kinetic method
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2003
PB - EDP-Sciences
VL - 37
IS - 2
SP - 389
EP - 416
AB - 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.
LA - eng
KW - shallow water equations; Saint-Venant system; finite volume method; kinetic scheme; transport of pollutant; time discretization
UR - http://eudml.org/doc/245541
ER -

References

top
  1. [1] E. Audusse, M.O. Bristeau and B. Perthame, 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. [2] A. Bermudez and M.E. Vasquez, Upwind methods for hyperbolic conservation laws with source terms. Comput. & Fluids 23 (1994) 1049–1071. Zbl0816.76052
  3. [3] M.O. Bristeau and B. Coussin, 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. [4] M.O. Bristeau and B. Perthame, Transport of pollutant in shallow water using kinetic schemes. CEMRACS, ESAIM Proc. 10 (1999) 9–21, http://www.emath.fr/Maths/Proc/Vol.10. Zbl1151.76593
  5. [5] R. Eymard, T. Gallouet and R. Herbin, Finite volume methods, Handbook of numerical analysis, Vol. VIII, P.G. Ciarlet and J.L. Lions Eds., Amsterdam, North-Holland (2000). Zbl0981.65095MR1804748
  6. [6] T. Gallouet, J.M. Hérard and N. Seguin, Some approximate Godunov schemes to compute shallow water equations with topography. Comput. & Fluids 32 (2003) 479–513. Zbl1084.76540
  7. [7] J.F. Gerbeau and B. Perthame, Derivation of viscous Saint-Venant system for laminar shallow water; Numerical validation. Discrete Contin. Dynam. Systems 1 (2001) 89–102. Zbl0997.76023
  8. [8] E. Godlewski and P.A. Raviart, Numerical approximation of hyperbolic systems of conservation laws. Springer-Verlag, New York, Appl. Math. Sci. 118 (1996). Zbl0860.65075MR1410987
  9. [9] L. Gosse and A.Y. LeRoux, A well-balanced scheme designed for inhomogeneous scalar conservation laws. C. R. Acad. Sci. Paris Sér. I Math. 323 (1996) 543–546. Zbl0858.65091
  10. [10] J.M. Hervouet, Hydrodynamique des écoulements à surface libre, apport de la méthode des éléments finis. EDF (2001). 
  11. [11] S. Jin, A steady state capturing method for hyperbolic system with geometrical source terms. ESAIM: M2AN 35 (2001) 631–646. Zbl1001.35083
  12. [12] R.J. LeVêque, Numerical Methods for Conservation Laws. Second edition, ETH Zurich, Birkhauser, Lectures in Mathematics (1992). Zbl0723.65067MR1153252
  13. [13] R.J. LeVêque, Balancing source terms and flux gradients in high-resolution Godunov methods: the quasi-steady wave-propagation algorithm. J. Comput. Phys. 146 (1998) 346–365. Zbl0931.76059
  14. [14] L. Martin, 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. [15] B. Perthame, Kinetic formulations of conservation laws. Oxford University Press (2002). Zbl1030.35002MR2064166
  16. [16] B. Perthame and C. Simeoni, A kinetic scheme for the Saint-Venant system with a source term. Calcolo 38 (2001) 201–231. Zbl1008.65066
  17. [17] P.L. Roe, 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. Zbl0626.65086
  18. [18] A.J.C. de Saint-Venant, 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. Zbl03.0482.04JFM03.0482.04
  19. [19] J.J. Stoker, The formation of breakers and bores. Comput. Appl. Math. 1 (1948). Zbl0041.54602MR24307

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.