数学物理学报(英文版) ›› 2024, Vol. 44 ›› Issue (6): 2307-2340.doi: 10.1007/s10473-024-0615-z

• • 上一篇    下一篇

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

Zhen LI, Min LIU, Dingfang LI   

  1. School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
  • 收稿日期:2023-08-07 修回日期:2024-07-15 发布日期:2024-12-06
  • 通讯作者: † Dingfang LI, E-mail: di@whu.edu.cn
  • 作者简介:Zhen LI, E-mail: zhen.li@whu.edu.cn; Min LIU, E-mail: liumin@whu.edu.cn
  • 基金资助:
    National Natural Science Foundation of China (51879194).

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

Zhen LI, Min LIU, Dingfang LI   

  1. School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
  • Received:2023-08-07 Revised:2024-07-15 Published:2024-12-06
  • Contact: † Dingfang LI, E-mail: di@whu.edu.cn
  • About author:Zhen LI, E-mail: zhen.li@whu.edu.cn; Min LIU, E-mail: liumin@whu.edu.cn
  • Supported by:
    National Natural Science Foundation of China (51879194).

摘要: 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.

关键词: Ripa model, moving-water steady states, nonstaggered central scheme, flux globalization, Runge-Kutta solvers

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.

Key words: Ripa model, moving-water steady states, nonstaggered central scheme, flux globalization, Runge-Kutta solvers

中图分类号: 

  • 76M12