An Energy-Stable High-Order Central Difference Scheme for the Two-Dimensional Shallow Water Equations

被引:0
作者
Matthew Brown
Margot Gerritsen
机构
[1] University of Auckland,Department of Engineering Science, School of Engineering
[2] Stanford University,Department of Petroleum Engineering, Green Earth Sciences Building
来源
Journal of Scientific Computing | 2006年 / 28卷
关键词
Shallow water equations; entropy-splitting; semi-discrete energy estimate; central finite difference scheme;
D O I
暂无
中图分类号
学科分类号
摘要
An energy-stable high-order central finite difference scheme is derived for the two-dimensional shallow water equations. The scheme is mathematically formulated using the semi-discrete energy method for initial boundary value problems described in Olsson (1995, Math. Comput. 64, 1035–1065): after symmetrizing the equations via a change to entropy variables, the flux derivatives are entropy-split enabling the formulation of a semi-discrete energy estimate. We show experimentally that the entropy-splitting improves the stability properties of the fully discretized equations. Thus, the dependence on numerical dissipation to keep the scheme stable for long term time integrations is reduced relative to the original unsplit form, thereby decreasing non-physical damping of solutions. The numerical dissipation term used with the entropy-split equations is in a form which preserves the semi-discrete energy estimate. A random one-dimensional dam break calculation is performed showing that the shock speed is computed correctly for this particular case, however it is an open question whether the correct shock speed will be computed in general
引用
收藏
页码:1 / 30
页数:29
相关论文
共 29 条
[1]  
Bova S.W.(1996)An entropy variable formulation and applications for the two dimensional shallow water equations Int. J. Numerical Methods Fluids 24 29-46
[2]  
Carey G.F.(1996)A symmetric formulation and SUPG scheme for the shallow water equations Adv. Water Res. 19 123-131
[3]  
Bova S.W.(1998)Designing an efficient solution strategy for fluid flows. I. A stable high-order finite difference scheme and sharp shock resolution for the Euler equations J. Comput. Phys. 129 245-262
[4]  
Carey G.F.(1998)Designing an efficient solution strategy for fluid flows. II. Stable high-order central finite difference schemes on composite adaptive grids with sharp shock resolution J. Comput. Phys. 147 293-317
[5]  
Gerritsen M.(1995)Fourth–order difference methods for hyperbolic IBVPs J. Comput. Phys. 117 300-164
[6]  
Olsson P.(1984)On the symmetric form of systems of conservation laws with entropy J. Comput. Phys. 49 151-112
[7]  
Gerritsen M.(1998)A symmetric formulation for computing transient shallow water flows Comput. Methods Appl. Mech. Eng. 163 111-234
[8]  
Olsson P.(1986)A new finite element formulation for computational fluid dynamics: 1. Symmetric forms of the compressible Euler and Navier-Stokes equations and the second law of thermodynamics Comput. Methods Appl. Mech. Eng. 54 223-210
[9]  
Gustafsson B.(1990)Theoretical and practical aspects of some initial boundary value problems in fluid dynamics J. Comput. Phys 86 187-446
[10]  
Olsson P.(1978)Summation-by-parts, projections and stability SIAM J. Appl. Math. 35 419-1065