Flux globalization based well-balanced path-conservative central-upwind scheme for two-layer thermal rotating shallow water equations

Cao，Yangyang1; Kurganov，Alexander2; Liu，Yongle1,3; Zeitlin，Vladimir4

2023-02-01

10.1016/j.jcp.2022.111790

JOURNAL OF COMPUTATIONAL PHYSICS   Impact Factor & Quartile

0021-9991

1090-2716

We develop a flux globalization based well-balanced path-conservative central-upwind scheme for the two-layer thermal rotating shallow water (TRSW) equations, which arise in both oceanography and atmospheric sciences. As it occurs in all of the two-layer shallow water models, the studied two-layer TRSW system loses hyperbolicity due to the well-known Kelvin-Helmholtz type instabilities arising when the vertical velocity shear between the layers becomes sufficiently large. Therefore, prior to developing the numerical method, we discuss the hyperbolicity criterion. Besides the hyperbolicity loss, there are additional challenges in the development of numerical methods for the studied system. They are related to the presence of nonconservative terms modeling layer interaction and also to a quite complicated structure of steady-states solutions a good well-balanced numerical method should be able to exactly preserve. In order to treat the nonconservative product terms, we use the path-conservative technique, which is implemented within the flux globalization framework: the source and nonconservative terms are incorporated into the fluxes, which results in a quasi-conservative system, which is then numerically solved using the Riemann-problem-solver-free central-upwind scheme. A well-balanced property of the resulting scheme is ensured by performing piecewise linear reconstruction for the equilibrium rather than the conservative variables, by development of special quadratures required for the flux globalization procedure, and by switching off a part of the numerical diffusion when the computed solution is near or at thermo-geostrophic equilibria. The advantages and excellent performance of the proposed scheme are demonstrated on a number of numerical examples.

Flux globalization Path-conservative central-upwind schemes Two-layer thermal rotating shallow water equations Well-balanced schemes

English

First ; Corresponding

National Natural Science Foundation of China[12111530004];National Natural Science Foundation of China[12171226];Guangdong Provincial Key Laboratory Of Computational Science And Material Design[2019B030301001];

PHYSICS

2-s2.0-85142678329

Scopus

1.Department of Mathematics,Southern University of Science and Technology,Shenzhen,518055,China
2.Department of Mathematics,SUSTech International Center for Mathematics and Guangdong Provincial Key Laboratory of Computational Science and Material Design,Southern University of Science and Technology,Shenzhen,518055,China
3.Institute of Mathematics,University of Zürich,Zürich,8057,Switzerland
4.Laboratory of Dynamical Meteorology,Sorbonne University,Ecole Normale Supérieure,CNRS,Paris,24 rue Lhomond,75005,France

Cao，Yangyang,Kurganov，Alexander,Liu，Yongle,et al. Flux globalization based well-balanced path-conservative central-upwind scheme for two-layer thermal rotating shallow water equations[J]. JOURNAL OF COMPUTATIONAL PHYSICS,2023,474.

Cao，Yangyang,Kurganov，Alexander,Liu，Yongle,&Zeitlin，Vladimir.(2023).Flux globalization based well-balanced path-conservative central-upwind scheme for two-layer thermal rotating shallow water equations.JOURNAL OF COMPUTATIONAL PHYSICS,474.

Cao，Yangyang,et al."Flux globalization based well-balanced path-conservative central-upwind scheme for two-layer thermal rotating shallow water equations".JOURNAL OF COMPUTATIONAL PHYSICS 474(2023).