A fast finite volume solver for multi-layered shallow water flows with mass exchange

JOURNAL OF COMPUTATIONAL PHYSICS

卷:272

Article

Audusse, Emmanuel^1,2  Benkhaldoun, Fayssal²  Sari, Saida²  Seaid, Mohammed³  Tassi, Pablo⁴ 

[1] INRIA Paris Rocquencourt, F-78153 Le Chesnay, France

[2] Univ Paris 13, LAGA, Sorbonne Paris Cite, F-93430 Villetaneuse, France

[3] Univ Durham, Sch Engn & Comp Sci, Durham DH1 3HP, England

[4] Lab Natl Hydraul & Environm, F-78401 Chatou, France

关键词: Multi-layer shallow water equations;    Hydrostatic flows;    Finite volume method;    Method of characteristics;    Mass exchange;    Wind-driven flows;   

DOI  :  10.1016/j.jcp.2014.04.026

来源: Elsevier

【 摘 要 】

A fast finite volume solver for hydrostatic multi-layered shallow water flows with mass exchange is investigated. In contrast to many models for multi-layered hydrostatic shallow water flows where the immiscible suppression is assumed, the present model allows for mass exchange between the layers. The multi-layered shallow water equations form a system of conservation laws with source terms for which the computation of the eigenvalues is not trivial. For most practical applications, complex eigenvalues may arise in the system and the multi-layered shallow water equations are not hyperbolic anymore. This property makes the application of conventional finite volume methods difficult or even impossible for those methods that require in their formulation the explicit computation of the eigenvalues. In the current study, we propose a finite volume method that avoids the solution of Riemann problems. At each time step, the method consists of two stages to update the new solution. In the first stage, the multi-layered shallow water equations are rewritten in a non-conservative form and the intermediate solutions are calculated using the method of characteristics. In the second stage, the numerical fluxes are reconstructed from the intermediate solutions in the first stage and used in the conservative form of the multi-layered shallow water equations. The proposed method is simple to implement, satisfies the conservation property and is suitable for multilayered shallow water equations on non-flat topography. The proposed finite volume solver is verified against several benchmark tests and it shows good agreement with analytical solutions of the incompressible hydrostatic Navier-Stokes equations. The method is conservative by construction and preserves the mass to the machine precision. The performance of the method is also demonstrated by comparing the results obtained using the proposed finite volume method to those obtained using the well-established kinetic method. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free