A MOVING WATER EQUILIBRIA PRESERVING NONSTAGGERED CENTRAL SCHEME ACHIEVED VIA FLUX GLOBALIZATION FOR THE RIPA MODEL

  • Zhen LI ,
  • Min LIU ,
  • Dingfang LI
Expand
  • School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
Zhen LI, E-mail: zhen.li@whu.edu.cn; Min LIU, E-mail: liumin@whu.edu.cn

Received date: 2023-08-07

  Revised date: 2024-07-15

  Online published: 2024-12-06

Supported by

National Natural Science Foundation of China (51879194).

Abstract

In this paper, we propose a second-order moving-water equilibria preserving nonstaggered central scheme to solve the Ripa model via flux globalization. To maintain the moving-water steady states, we use the discrete source terms proposed by Britton et al. (J Sci Comput, 2020, 82(2): Art 30) by incorporating the expression of the source terms as a whole into the flux gradient, which directly avoids the discrete complexity of the source terms in order to maintain the well-balanced properties of the scheme. In addition, since the nonstaggered central scheme requires re-projecting the updated values of the nonstaggered cells onto the staggered cells, we modify the calculation of the global variables by constructing ghost cells and alternating the values of the global variables with the water depths obtained from the solution through the nonlinear relationship between the global flux and the water depth. In order to maintain the second-order accuracy of the scheme on the time scale, we incorporate a new Runge-Kutta type time discretization in the evolution of the numerical solution for the nonstaggered cells. Meanwhile, we introduce the "draining" time step technique to ensure that the water depth is positive and that it satisfies mass conservation. Numerical experiments verify that the scheme is well-balanced, positivity-preserving and robust.

Cite this article

Zhen LI , Min LIU , Dingfang LI . A MOVING WATER EQUILIBRIA PRESERVING NONSTAGGERED CENTRAL SCHEME ACHIEVED VIA FLUX GLOBALIZATION FOR THE RIPA MODEL[J]. Acta mathematica scientia, Series B, 2024 , 44(6) : 2307 -2340 . DOI: 10.1007/s10473-024-0615-z

References

[1] Nessyahu H, Tadmor E.Non-oscillatory central differencing for hyperbolic conservation laws. J Comput Phys, 1990, 87: 408-463
[2] Liu X, Tadmor E.Third order nonoscillatory central scheme for hyperbolic conservation laws. Numer Math, 1998, 79: 397-425
[3] Levy D, Puppo G, Russo G.Central WENO schemes for hyperbolic systems of conservation laws. M2AN Math Model Numer Anal, 1999, 33: 547-571
[4] Arminjon P, Stanescu D, Viallon M.A two-dimensional finite volume extension of the Lax-Friedrichs and Nessyahu-Tadmor schemes for compressible flow. Proceedings of the 6th International Symposium on Computational Fluid Dynamics, 1994, 4: 7-14
[5] Arminjon P, Viallon M, Madrane A.A finite volume extension of the Lax-Friedrichs and Nessyahu-Tadmor schemes for conservation laws on unstructured grids. Int J Comput Fluid D, 1998, 9: 1-22
[6] Arminjon P, Viallon M C.Convergence of a finite volume extension of the Nessyahu-Tadmor scheme on unstructured grids for a two-dimensional linear hyperbolic equation. SIAM J Numer Analy, 1999, 36: 738-771
[7] Jiang G S, Tadmor E.Nonoscillatory central schemes for multidimensional hyperbolic conservation laws. SIAM J Sci Comput, 1998, 19: 1892-1917
[8] Levy D, Puppo G, Russo G.Compact central weno schemes for multidimensional conservation laws. SIAM J Sci Comput, 2000, 22: 656-672
[9] Kurganov A, Levy D.Central-upwind schemes for the Saint-Venant system. M2AN Math Model Numer Anal, 2002, 36: 397-425
[10] Kurganov A, Petrova G.A second-order well-balanced positivity preserving central-upwind scheme for the Saint-Venant system. Commun Math Sci, 2007, 5: 133-160
[11] Kurganov A, Tadmor E.New high-resolution central schemes for nonlinear conservation laws and convection-diffusion equations. J Comput Phys, 2000, 160: 241-282
[12] Jiang G S, Levy D, Lin C T, et al.High-resolution nonoscillatory central schemes with nonstaggered grids for hyperbolic conservation laws. SIAM J Numer Anal, 1998, 35: 2147-2168
[13] Touma R.Central unstaggered finite volume schemes for hyperbolic systems: Applications to unsteady shallow water equations. Appl Math Comput, 2009, 213: 47-59
[14] Touma R, Klingenberg C.Well-balanced central finite volume methods for the Ripa system. Appl Numer Math, 2015, 97: 42-68
[15] Dellar P J.Common Hamiltonian structure of the shallow water equations with horizontal temperature gradients and magnetic fields. Phys Fluids, 2003, 15: 292-297
[16] Ripa P.Conservation laws for primitive equations models with inhomogeneous layers. Geophys Astro Fluid, 1993, 70: 85-111
[17] Ripa P.On improving a one-layer ocean model with thermodynamics. J Fluid Mech, 1995, 303: 169-201
[18] Russo G.Central schemes for conservation laws with application to shallow water equations// Trends and Applications of Mathematics to Mechanics. Milan: Springer, 2005: 225-246
[19] Touma R, Khankan S.Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems. Appl Math Comput, 2012, 218: 5948-5960
[20] Bollermann A, Chen G X, Kurganov A, Noelle S.A well-balanced reconstruction for wet/dry fronts for the shallow water equations. J Sci Comput, 2013, 56: 267-290
[21] Touma R.Well-balanced central schemes for systems of shallow water equations with wet and dry states. Appl Math Model, 2016, 40: 2929-2945
[22] 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: Art 35
[23] Dong J, Li D F.An effect non-staggered central scheme based on new hydrostatic reconstruction. Appl Math Comput, 2020, 372: 124992
[24] Kounadis G, Dougalis V.Galerkin finite element methods for the shallow water equations over variable bottom. J Comput Appl Math, 2020, 373: 112315
[25] Dong J, Li D F.A reliable second-order hydrostatic reconstruction for shallow water flows with the friction term and the bed source term. J Comput Appl Math, 2020, 376: 112871
[26] Britton J, Xing Y. High order still-water and moving-water equilibria preserving discontinuous Galerkin methods for the Ripa model. J Sci Comput, 2020, 82: Art 30
[27] Li J, Li G, Qian S, et al.A high-order well-balanced discontinuous Galerkin method based on the hydrostatic reconstruction for the Ripa model. Adv Appl Math Mech, 2020, 12: 1416-1437
[28] Chertock A, Kurganov A, Liu Y.Central-upwind schemes for the system of shallow water equations with horizontal temperature gradients. Numer Math, 2014, 127: 595-639
[29] Han X, Li G. Well-balanced finite difference WENO schemes for the Ripa model. Comput Fluids, 2016, 134-135: 1-10
[30] Qian S, Shao F, Li G.High order well-balanced discontinuous Galerkin methods for shallow water flow under temperature fields. Comput Appl Math, 2018, 37: 5775-5794
[31] Dong J, Li D F.Exactly well-balanced positivity preserving nonstaggered central scheme for open-channel flows. Int J Numer Meth Fl, 2021, 93: 273-292
[32] Wang J L, Zhou Z T, Wang Y S.Local structure-preserving algorithms for the Klein-Gordon-Zakharov equation. Acta Mathematica Scientia, 2023, 43B: 1211-1238
[33] Chen T W, Zhang Y Q.A strong solution of Navier-Stokes equations with a rotation effect for isentropic compressible fluids. Acta Mathematica Scientia, 2021, 41B: 1579-1605
[34] Li Z, Dong J, Luo Y, et al.Moving water equilibria preserving nonstaggered central scheme for open-channel flows. Math Method Appl Sci, 2023, 46(6): 7391-7412
[35] Bouchut F, Morales T.A subsonic-well-balanced reconstruction scheme for shallow water flows. SIAM J Numer Anal, 2010, 48: 1733-1758
[36] Cheng Y, Kurganov A.Moving-water equilibria preserving central-upwind schemes for the shallow water equations. Commun Math Sci, 2016, 14: 1643-1663
[37] Cheng Y, Chertock A, Herty M, et al.A new approach for designing moving-water equilibria preserving schemes for the shallow water equations. J Sci Comput, 2019, 80: 538-554
[38] Liu X, Chen X, Jin S, et al.Moving-water equilibria preserving partial relaxation scheme for the Saint-Venant system. SIAM J Sci Comput, 2020, 42: A2206-A2229
[39] Noelle S, Xing Y L, Shu C W.High-order well-balanced finite volume WENO schemes for shallow water equation with moving water. J Comput Phys, 2007, 226: 29-58
[40] Xing Y L, Shu C W, Noelle S.On the advantage of well-balanced schemes for moving-water equilibria of the shallow water equations. J Sci Comput, 2011, 48: 339-349
[41] Russo G, Khe A.High order well balanced schemes for systems of balance laws. Proc Symp Appl Math Amer Math, 2009, 67: 919-928
[42] Levy D, Puppo G, Russo G.A fourth-order central WENO scheme for multidimensional hyperbolic systems of conservation laws. SIAM J Sci Comput, 2002, 24: 480-506
[43] Castro M, López-García J, Parés C.High order exactly well-balanced numerical methods for shallow water systems. J Comput Phys, 2013, 246: 242-264
[44] Cao Y, Kurganov A, Liu Y, Xin R. Flux globalization based well-balanced path-conservative central-upwind schemes for shallow water models. J Sci Comput, 2022, 92: Art 69
[45] Xing Y L.Exactly well-balanced discontinuous Galerkin methods for the shallow water equations with moving water equilibrium. J Comput Phys, 2014, 257: 536-553
[46] Xing Y L, Shu C W.A survey of high order schemes for the shallow water equations. J Math Study, 2014, 47: 221-249
[47] Li Z, Dong J, Luo Y, et al.Moving water equilibria preserving nonstaggered central scheme for open-channel flows. Math Method Appl Sci, 2023, 46: 7391-7412
[48] Castro M, Parés C. Well-balanced high-order finite volume methods for systems of balance laws. J Sci Comput, 2020, 82: Art 48
[49] Chertock A, Herty M, Özcan, Ş N.Well-balanced central-upwind schemes for $2\times 2$ systems of balance laws// Kingenberg C, Westdickenberg M. Theory Numerics Applications of Hyperbolic Problems I. Cham: Springer, 2016: 345-361
[50] Chertock A, Cui S, Kurganov A, et al.Well-balanced schemes for the Euler equations with gravitation: conservative formulation using global fluxes. J Comput Phys, 2018, 358: 36-52
[51] Chertock A, Kurganov A, Liu X, et al. Well-balancing via flux globalization: Applications to shallow water equations with wet/dry fronts. J Sci Comput, 2022, 90: Art 9
[52] Gascón L, Corberán J.Construction of second-order TVD schemes for nonhomogeneous hyperbolic conservation laws. J Comput Phys, 2001, 172: 261-297
[53] Caselles V, Donat R, Haro G.Flux-gradient and source-term balancing for certain high resolution shock-capturing schemes. Comput Fluids, 2009, 38: 16-36
[54] Donat R, Martínez-Gavara A.Hybrid second order schemes for scalar balance laws. J Sci Comput, 2011, 48: 52-69
[55] Ciallella M, Torlo D, Ricchiuto M. Arbitrary high order WENO finite volume scheme with flux globalization for moving equilibria preservation. J Sci Comput, 2023, 96: Art 53
[56] Mantri Y, Öffner P, Ricchiuto M. Fully well balanced entropy controlled discontinuous Galerkin spectral element method for shallow water flows: global flux quadrature and cell entropy correction. J Comp Phys, 2024, 498: Art 112673
[57] Lee S, Wright N.Simple and efficient solution of the shallow water equations with source terms. Inter J Numer Meth Fl, 2010, 63: 313-340
[58] Berthon C, Bulteau S, Foucher F, et al.A very easy high-order well balanced reconstruction for hyperbolic systems with source terms. SIAM J Sci Comput, 2022, 44: A2506-A2535
[59] Gómez-Bueno I, Castro M, Parés C. High-order well-balanced methods for systems of balance laws: a control-based approach. Appl Math Comput, 2021, 394: Art 125820
[60] Gómez-Bueno I, Díaz M, Parés C, Russo G. Collocation methods for high-order well-balanced methods for systems of balance laws. Mathematics, 2021, 9: Art 1799
[61] Pareschi L, Puppo G, Russo G.Central runge-kutta schemes for conservation laws. SIAM J Sci Comput, 2005, 26: 979-999
[62] Van Leer B.Towards the ultimate conservative difference scheme. V. A second-order sequel to godunov's method. J Comput Phys, 1979, 32: 101-136
[63] Audusse E, Bouchut F, Bristeau M, et al.A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows, SIAM J Sci Comput, 2004, 25: 2050-2065
[64] Asher U, Ruuth S, Spiteri R.Implicit-explicit Runge-Kutta methods for time dependent partial differential equations. Appl Numer Math, 1997, 25: 151-167
[65] Bollermann A, Noelle S, Lukáčová-Medvid'ová M. Finite volume evolution Galerkin methods for the shallow water equations with dry beds. Commun Comput Phys, 2011, 10: 371-404
[66] Xing Y, Shu C.High order finite difference WENO schemes with the exact conservation property for the shallow water equations. J Comput Phys, 2005, 208: 206-227
[67] Khan A, Lai W.Modeling Shallow Water Flows Using the Discontinuous Galerkin Method. New York: CRC Press, 2014
[68] Vázquez-Cendón M.Improved treatment of source terms in upwind schemes for the shallow water equations in channels with irregular geometry. J Comput Phys, 1999, 148: 497-526
Options
Outlines

/