<p>This paper proposes a robust structure-preserving dissipation-reducing nonstaggered central scheme for the multi-dimensional two-layer shallow water equations with wet-dry fronts. Since the two-layer shallow water equations are conditionally hyperbolic, computing Riemann solutions for numerical fluxes across cell interfaces is challenging. This issue is circumvented by the proposed scheme, which is genuinely Riemann-problem-solver-free, meaning that the numerical fluxes are independent of characteristic speed information. The source term and non-conservative product arising from the momentum exchange, are discretized using auxiliary variables. One key advantage of the scheme is its ability to preserve stationary solutions, even in the presence of wet-dry fronts in the computational domain. The auxiliary variables are interpreted as point values along paths in the phase space, which are used to discretize the non-conservative product, reinterpreted as a Borel measure. We prove that the scheme is well-balanced and positivity-preserving. Specifically, we establish a Lax-Wendroff type convergence theorem, in which the dissipation-reducing parameter plays a crucial and nontrivial role. Finally, we apply the scheme to several classical problems of the two-layer shallow water equations, with numerical results confirming its efficiency and robustness.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

A Robust Structure-Preserving Nonstaggered Central Scheme for Multidimensional Two-Layer Shallow Water Equations with Wetting and Drying Transitions

  • Xingpo Qiu,
  • Jian Dong,
  • Xu Qian,
  • Zige Wei

摘要

This paper proposes a robust structure-preserving dissipation-reducing nonstaggered central scheme for the multi-dimensional two-layer shallow water equations with wet-dry fronts. Since the two-layer shallow water equations are conditionally hyperbolic, computing Riemann solutions for numerical fluxes across cell interfaces is challenging. This issue is circumvented by the proposed scheme, which is genuinely Riemann-problem-solver-free, meaning that the numerical fluxes are independent of characteristic speed information. The source term and non-conservative product arising from the momentum exchange, are discretized using auxiliary variables. One key advantage of the scheme is its ability to preserve stationary solutions, even in the presence of wet-dry fronts in the computational domain. The auxiliary variables are interpreted as point values along paths in the phase space, which are used to discretize the non-conservative product, reinterpreted as a Borel measure. We prove that the scheme is well-balanced and positivity-preserving. Specifically, we establish a Lax-Wendroff type convergence theorem, in which the dissipation-reducing parameter plays a crucial and nontrivial role. Finally, we apply the scheme to several classical problems of the two-layer shallow water equations, with numerical results confirming its efficiency and robustness.