Abstract:This paper develops high-order well-balanced (WB) energy stable (ES) finite difference schemes for multi-layer (the number of layers $M\geqslant 2$) shallow water equations (SWEs) on both fixed and adaptive moving meshes, extending our previous works [20,51]. To obtain an energy inequality, the convexity of an energy function for an arbitrary $M$ is proved by finding recurrence relations of the leading principal minors or the quadratic forms of the Hessian matrix of the energy function with respect to the conservative variables, which is more involved than the single-layer case due to the coupling between the layers in the energy function. An important ingredient in developing high-order semi-discrete ES schemes is the construction of a two-point energy conservative (EC) numerical flux. In pursuit of the WB property, a sufficient condition for such EC fluxes is given with compatible discretizations of the source terms similar to the single-layer case. It can be decoupled into $M$ identities individually for each layer, making it convenient to construct a two-point EC flux for the multi-layer system. To suppress possible oscillations near discontinuities, WENO-based dissipation terms are added to the high-order WB EC fluxes, which gives semi-discrete high-order WB ES schemes. Fully-discrete schemes are obtained by employing high-order explicit SSP-RK methods and proved to preserve the lake at rest. The schemes are further extended to moving meshes based on a modified energy function for a reformulated system, relying on the techniques proposed in [51]. Numerical experiments are conducted for some two- and three-layer cases to validate the high-order accuracy, WB and ES properties, and high efficiency of the schemes, with a suitable amount of dissipation chosen by estimating the maximal wave speed due to the lack of an analytical expression for the eigenstructure of the multi-layer system.

High-order well-balanced central WENO scheme for pre-balanced shallow water equations

Well-balanced fifth-order finite volume WENO schemes with constant subtraction technique for shallow water equations

A new fifth-order finite difference well-balanced multi-resolution WENO scheme for solving shallow water equations

Well-balanced fifth-order finite difference Hermite WENO scheme for the shallow water equations

Well-balanced finite difference WENO-AO scheme for rotating shallow water equations with Coriolis force

High order well-balanced asymptotic preserving finite difference WENO schemes for the shallow water equations in all Froude numbers

A Consistent and Well‐Balanced Hybrid Weighted Essentially Non‐Oscillatory Scheme for Shallow Water Equations on Unstructured Meshes

Third and fourth order well-balanced schemes for the shallow water equations based on the CWENO reconstruction

Hybrid fifth-order unequal-sized weighted essentially non-oscillatory scheme for shallow water equations

A Mass Conservative, Well Balanced and Positivity-Preserving Central Scheme for Shallow Water Equations

High-order WENO finite-difference methods for hyperbolic nonconservative systems of Partial Differential Equations

Well-Balanced Central Schemes on Overlapping Cells with Constant Subtraction Techniques for the Saint-Venant Shallow Water System

Combined WENO schemes for increasing the accuracy of the numerical solution of conservation laws

Well-balanced, anti-diffusive non-oscillatory central difference (adNOC) scheme for the shallow water equations with wet-dry fronts

Central WENO schemes through a global average weight

High-order accurate well-balanced energy stable finite difference schemes for multi-layer shallow water equations on fixed and adaptive moving meshes

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

High-Order Oscillation-Eliminating Hermite WENO Method for Hyperbolic Conservation Laws

A Well-Balanced Reconstruction of Wet/Dry Fronts for the Shallow Water Equations