Instability in Leapfrog and Forward-Backward Schemes, Part II-Numerical Simulations of Dam Break

Abstract

Following Sun’s approach [17], Shuman smoothing instead of conventional diffusion terms is used in a simple two-time step semi-implicit finite volume scheme to simulate dam break. When the Courant number is less than one, the absolute value of amplification factor of the 1D linearized shallow-water equations is 1 in this new scheme. Compared with the characteristic-based semi-Lagrangian schemes and the Riemann solver, this scheme produces excellent results of free water depth and speed of the shock. Numerical simulations show that the water inside the dam initially moves away radially until water almost depletes near the center; then the water moves back to the center and forms a vertical water column there. This paper proves that Shuman smoothing can be used not only in the linearized shallow-water equations discussed in Sun [17] but also in the nonlinear wave equations to control instability around shocks.

Keywords

Riemann solver, shallow-water eqyatuibsm senu-implicit scheme, Eigenvalue, characteristic-base, semi-Lagranian, shock

Date of this Version

2010

DOI

http://dx.doi.org/10.1016/j.compfluid.2010.11.002,

Volume

45

Issue

1

Pages

70-76

Link Out to Full Text

http://www.sciencedirect.com/science/article/pii/S0045793010003063

Link to Full Text

Share

COinS