A novel multidimensional solution reconstruction and edge-based limiting procedure for unstructured cell-centered finite volumes with application to shallow water dynamics

被引:34
作者
Delis, A. I. [1 ]
Nikolos, I. K. [2 ]
机构
[1] Tech Univ Crete, Dept Sci, Div Math, Khania 73100, Crete, Greece
[2] Tech Univ Crete, Dept Prod Engn & Management, Khania 73100, Crete, Greece
关键词
cell-centered finite volumes; linear reconstruction; limiting; shallow water equations; source terms; wet/dry fronts; HIGH-RESOLUTION SCHEMES; RIEMANN SOLVERS; SOURCE TERMS; RESIDUAL DISTRIBUTION; BOUNDARY-CONDITIONS; NUMERICAL TREATMENT; CONSERVATION-LAWS; WET/DRY FRONTS; UPWIND SCHEMES; EQUATIONS;
D O I
10.1002/fld.3674
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
On unstructured meshes, the cell-centered finite volume (CCFV) formulation, where the finite control volumes are the mesh elements themselves, is probably the most used formulation for numerically solving the two-dimensional nonlinear shallow water equations and hyperbolic conservation laws in general. Within this CCFV framework, second-order spatial accuracy is achieved with a Monotone Upstream-centered Schemes for Conservation Laws-type (MUSCL) linear reconstruction technique, where a novel edge-based multidimensional limiting procedure is derived for the control of the total variation of the reconstructed field. To this end, a relatively simple, but very effective modification to a reconstruction procedure for CCFV schemes, is introduced, which takes into account geometrical characteristics of computational triangular meshes. The proposed strategy is shown not to suffer from loss of accuracy on grids with poor connectivity. We apply this reconstruction in the development of a second-order well-balanced Godunov-type scheme for the simulation of unsteady two-dimensional flows over arbitrary topography with wetting and drying on triangular meshes. Although the proposed limited reconstruction is independent from the Riemann solver used, the well-known approximate Riemann solver of Roe is utilized to compute the numerical fluxes, whereas the GreenGauss divergence formulation for gradient computations is implemented. Two different stencils for the GreenGauss gradient computations are implemented and critically tested, in conjunction with the proposed limiting strategy, on various grid types, for smooth and nonsmooth flow conditions. Copyright (c) 2012 John Wiley & Sons, Ltd.
引用
收藏
页码:584 / 633
页数:50
相关论文
共 80 条
  • [1] Aftosmis, 1994, 940415 AIAA
  • [2] BEHAVIOR OF LINEAR RECONSTRUCTION TECHNIQUES ON UNSTRUCTURED MESHES
    AFTOSMIS, M
    GAITONDE, D
    TAVARES, TS
    [J]. AIAA JOURNAL, 1995, 33 (11) : 2038 - 2049
  • [3] [Anonymous], 2002, Cambridge Texts in Applied Mathematics, DOI [10.1017/CBO9780511791253, DOI 10.1017/CBO9780511791253]
  • [4] A well-balanced positivity preserving "second-order" scheme for shallow water flows on unstructured meshes
    Audusse, E
    Bristeau, MO
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2005, 206 (01) : 311 - 333
  • [5] Barth T., 2004, ENCY COMPUTATIONAL M
  • [6] Barth T.J., 1989, 27 AEROSPACE SCI M
  • [7] Batten P, 1996, INT J NUMER METH ENG, V39, P1821, DOI 10.1002/(SICI)1097-0207(19960615)39:11<1821::AID-NME929>3.0.CO
  • [8] 2-E
  • [9] Adaptive Godunov-based model for flood simulation
    Begnudelli, Lorenzo
    Sanders, Brett F.
    Bradford, Scott F.
    [J]. JOURNAL OF HYDRAULIC ENGINEERING, 2008, 134 (06) : 714 - 725
  • [10] Berger M, 2005, 43 AIAA AER SCI M EX