No Arabic abstract
We develop well-balanced central schemes on overlapping cells for the Saint-Venant shallow water system and its variants. The main challenge in deriving the schemes is related to the fact that the Saint-Venant system contains a geometric source term due to nonflat bottom topography and therefore a delicate balance between the flux gradients and source terms has to be preserved. We propose a constant subtraction technique, which helps one to ensure a well-balanced property of the schemes, while maintaining arbitrary high-order of accuracy. Hierarchical reconstruction limiting procedure is applied to eliminate spurious oscillations without using characteristic decomposition. Extensive one- and two-dimensional numerical simulations are conducted to verify the well-balanced property, high-order of accuracy, and non-oscillatory high-resolution for both smooth and nonsmooth solutions.
This note aims at demonstrating the advantage of moving-water well-balanced schemes over still-water well-balanced schemes for the shallow water equations. We concentrate on numerical examples with solutions near a moving-water equilibrium. For such examples, still-water well-balanced methods are not capable of capturing the small perturbations of the moving-water equilibrium and may generate significant spurious oscillations, unless an extremely refined mesh is used. On the other hand, moving- water well-balanced methods perform well in these tests. The numerical examples in this note clearly demonstrate the importance of utilizing moving-water well-balanced methods for solutions near a moving-water equilibrium.
The standard multilayer Saint-Venant system consists in introducing fluid layers that are advected by the interfacial velocities. As a consequence there is no mass exchanges between these layers and each layer is described by its height and its average velocity. Here we introduce another multilayer system with mass exchanges between the neighborhing layers where the unknowns are a total height of water and an average velocity per layer. We derive it from Navier-Stokes system with an hydrostatic pressure and prove energy and hyperbolicity properties of the model. We also give a kinetic interpretation leading to effective numerical schemes with positivity and energy properties. Numerical tests show the versatility of the approach and its ability to compute recirculation cases with wind forcing.
In this paper, high order semi-implicit well-balanced and asymptotic preserving finite difference WENO schemes are proposed for the shallow water equations with a non-flat bottom topography. We consider the Froude number ranging from O(1) to 0, which in the zero Froude limit becomes the lake equations for balanced flow without gravity waves. We apply a well-balanced finite difference WENO reconstruction, coupled with a stiffly accurate implicit-explicit (IMEX) Runge-Kutta time discretization. The resulting semi-implicit scheme can be shown to be well-balanced, asymptotic preserving (AP) and asymptotically accurate (AA) at the same time. Both one- and two-dimensional numerical results are provided to demonstrate the high order accuracy, AP property and good performance of the proposed methods in capturing small perturbations of steady state solutions.
In this paper, we develop a well-balanced oscillation-free discontinuous Galerkin (OFDG) method for solving the shallow water equations with a non-flat bottom topography. One notable feature of the constructed scheme is the well-balanced property, which preserves exactly the hydrostatic equilibrium solutions up to machine error. Another feature is the non-oscillatory property, which is very important in the numerical simulation when there exist some shock discontinuities. To control the spurious oscillations, we construct an OFDG method with an extra damping term to the existing well-balanced DG schemes proposed in [Y. Xing and C.-W. Shu, CICP, 1(2006), 100-134.]. With a careful construction of the damping term, the proposed method achieves both the well-balanced property and non-oscillatory property simultaneously without compromising any order of accuracy. We also present a detailed procedure for the construction and a theoretical analysis for the preservation of the well-balancedness property. Extensive numerical experiments including one- and two-dimensional space demonstrate that the proposed methods possess the desired properties without sacrificing any order of accuracy.
The blood flow model maintains the steady state solutions, in which the flux gradients are non-zero but exactly balanced by the source term. In this paper, we design high order finite difference weighted non-oscillatory (WENO) schemes to this model with such well-balanced property and at the same time keeping genuine high order accuracy. Rigorous theoretical analysis as well as extensive numerical results all indicate that the resulting schemes verify high order accuracy, maintain the well-balanced property, and keep good resolution for smooth and discontinuous solutions.