摘要
In this work,we present a high-order discontinuous Galerkin method for the shallow water equations incorporating horizontal temperature gradients(also known as the Ripa model),which exactly maintains the lake at rest steady state.Herein,we propose original numerical fluxes defined on the basis of the hydrostatic reconstruction idea and a simple source term approximation.This novel approach allows us to achieve the well-balancing of the discontinuous Galerkin method without complication.Moreover,the proposed method retains genuinely high-order accuracy for smooth solutions and it shows good resolution for discontinuous solutions at the same time.Rigorous numerical analysis as well as extensive numerical results all verify the good performances of the proposed method.
基金
support of the Natural Science Foundation of China through Grants No.11771228
The author Qiang Niu is supported by the XJTLU research enhancement fund with No.REF-18-01-04 and the Key Programme Special Fund(KSF)in XJTLU with Nos.KSF-E-32,KSF-E-21 and KSF-P-02.