Structure-Preserving Nonstaggered Central Difference Schemes at Wet-Dry Fronts for the Shallow Water Equations

Abstract

We aim to propose a Riemann-problem-solver-free nonstaggered central difference scheme (NCDS), which is well-balanced and positivity-preserving for two-dimensional shallow water flows over nonflat bottom topography with wet-dry fronts. The well-balanced property of the NCDS for the two-dimensional shallow water equations at wet-dry fronts is still unclear. The NCDS is “generously multidimensional” and consists of three steps: a forward step, a corrector step, and a backward step. Each step needs to guarantee the water depth to be nonnegative and retain the lake at rest. The key ingredient of the NCDS is the backward step. It is very nontrivial and technical to obtain both the well-balanced and positivity-preserving properties when the computational domain contains wet-dry fronts for the NCDS. The traditional NCDS fails to retain the lake at rest when the computational domain contains wet-dry fronts. We introduce a novel backward step, which is well-balanced and positivity-preserving, to recover numerical solutions on nonstaggered cells. Finally, we show several numerical results of two-dimensional shallow water equations with wet-dry fronts to verify these properties.

Keywords

Two-dimensional shallow water equations / Modifications at wet-dry fronts / Nonstaggered central difference schemes (NCDSs) / Structure-preserving / 76M12 / 35L65 / 65L05

References

[1] Audusse E, Bouchut F, Bristeau MO, Klein R, Perthame B. A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows. SIAM J. Sci. Comput., 2004, 25(6): 2050-2065

[2] Bollermann A, Noelle S, Lukáčová-Medvidová M. Finite volume evolution Galerkin methods for the shallow water equations with dry beds. Commun. Comput. Phys., 2011, 10(2): 371-404

[3] Briggs MJ, Synolakis CE, Harkins GS, Green DR. Laboratory experiments of tsunami runup on a circular island. Pure Appl. Geophys., 1995, 144(3): 569-593

[4] Bryson S, Epshteyn Y, Kurganov A, Petrova G. Well-balanced positivity preserving central-upwind scheme on triangular grids for the Saint-Venant system. ESAIM Math. Model. Numer. Anal., 2011, 45(3): 423-446

[5] Chen G, Noelle S. A new hydrostatic reconstruction scheme based on subcell reconstructions. SIAM J. Numer. Anal., 2017, 55(2): 758-784

[6] de Saint-Venant A. Théorie du mouvement non-permanent des eaux avec application aux crues des rivières et à i’introduction des warées dans leur lit. C. R. Acad. Sci. Paris, 1871, 73(99): 148-154

[7] Dong J. A robust second-order surface reconstruction for shallow water flows with a discontinuous topography and a Manning friction. Adv. Comput. Math., 2020, 46(2): 1-33

[8] Dong J, Fang Li D. Exactly well-balanced positivity preserving nonstaggered central scheme for open-channel flows. Int. J. Numer. Meth. Fluids, 2021, 93(1): 273-292

[9] Dong J, Li DF. An effect non-staggered central scheme based on new hydrostatic reconstruction. Appl. Math. Comput., 2020, 372: 124992

[10] Dong J, Li DF. Well-balanced nonstaggered central schemes based on hydrostatic reconstruction for the shallow water equations with Coriolis forces and topography. Math. Methods Appl. Sci., 2021, 44(2): 1358-1376

[11] Dong J, Qian X, Song S. Adaptive physical-constraints-preserving unstaggered central schemes for shallow water equations on quadrilateral meshes. ESAIM Math. Model. Numer. Anal., 2022, 56(6): 2297-2338

[12] Gottlieb S, Shu C-W, Tadmor E. Strong stability-preserving high-order time discretization methods. SIAM Rev., 2001, 43(1): 89-112

[13] Jiang GS, Levy D, Lin CT, Osher S, Tadmor E. High-resolution nonoscillatory central schemes with nonstaggered grids for hyperbolic conservation laws. SIAM J. Numer. Anal., 1998, 35(6): 2147-2168

[14] Kurganov A, Petrova G. A second-order well-balanced positivity preserving central-upwind scheme for the Saint-Venant system. Commun. Math. Sci., 2007, 5(1): 133-160

[15] Liu X, Albright J, Epshteyn Y, Kurganov A. Well-balanced positivity preserving central-upwind scheme with a novel wet/dry reconstruction on triangular grids for the Saint-Venant system. J. Comput. Phys., 2018, 374: 213-236

[16] Noellea S, Pankratza N. Well-balanced finite volume schemes of arbitrary order of accuracy for shallow water flows. J. Comput. Phys., 2006, 213(2): 474-499

[17] Thacker WC. Some exact solutions to the nonlinear shallow-water wave equations. J. Fluid Mech., 1981, 107: 499-508

[18] Touma R. Well-balanced central schemes for systems of shallow water equations with wet and dry states. Appl. Math. Model., 2016, 40(4): 2929-2945

[19] Touma RG, Kanbar F. Well-balanced central schemes for two-dimensional systems of shallow water equations with wet and dry states. Appl. Math. Model., 2018, 62: 728-750

[20] Van Leer B. Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method. J. Comput. Phys., 1997, 32(1): 101-136

[21] Xing Y. Exactly well-balanced discontinuous Galerkin methods for the shallow water equations with moving water equilibrium. J. Comput. Phys., 2014, 257(2): 536-553

[22] Xing Y, Shu C-W. A new approach of high order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source. Commun. Comput. Phys., 2006, 1(1): 567-598

[23] Xing Y, Zhang X, Shu C-W. Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations. Adv. Water Resour., 2010, 33(12): 1476-1493

[24] Zhao J, Özgen I, Liang D, Hinkelmann R. Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations. Int. J. Numer. Method Fluids, 2018, 87(8): 401-436

[25] Zhou JG, Causon DM, Mingham CG, Ingram DM. The surface gradient method for the treatment of source terms in the shallow-water equations. J. Comput. Phys., 2001, 168(1): 1-25

RIGHTS & PERMISSIONS

Shanghai University