Well-balanced hybrid compact-WENO scheme for shallow water equations

Abstract We investigate the performance of the high order well-balanced hybrid compact-weighted essentially non-oscillatory (WENO) finite difference scheme (Hybrid) for simulations of shallow water equations with source terms due to a non-flat bottom topography. The Hybrid scheme employs the nonlinear fifth order characteristic-wise WENO-Z finite difference scheme to capture high gradients and discontinuities in an essentially non-oscillatory manner, and the linear spectral-like sixth order compact finite difference scheme to resolve the fine scale structures in the smooth regions of the solution efficiently and accurately. The high order multi-resolution analysis is employed to identify the smoothness of the solution at each grid point. In this study, classical one- and two-dimensional simulations, including a long time two-dimensional dam-breaking problem with a non-flat bottom topography, are conducted to demonstrate the performance of the hybrid scheme in terms of the exact conservation property (C-property), good resolution and essentially non-oscillatory shock capturing of the smooth and discontinuous solutions respectively, and up to 2–3 times speedup factor over the well-balanced WENO-Z scheme.

[1]  Yulong Xing,et al.  Positivity-Preserving Well-Balanced Discontinuous Galerkin Methods for the Shallow Water Equations on Unstructured Triangular Meshes , 2013, J. Sci. Comput..

[2]  M. Hanif Chaudhry,et al.  Explicit Methods for 2‐D Transient Free Surface Flows , 1990 .

[3]  W. Don,et al.  High order Hybrid central-WENO finite difference scheme for conservation laws , 2007 .

[4]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[5]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[6]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[7]  R. LeVeque,et al.  Balancing Source Terms and Flux Gradientsin High-Resolution Godunov Methods : The Quasi-Steady Wave-Propogation AlgorithmRandall , 1998 .

[8]  Fayssal Benkhaldoun,et al.  Exact solutions to the Riemann problem of the shallow water equations with a bottom step , 2001 .

[9]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[10]  Peng Li,et al.  Hybrid Compact-WENO Finite Difference Scheme with Conjugate Fourier Shock Detection Algorithm for Hyperbolic Conservation Laws , 2016, SIAM J. Sci. Comput..

[11]  Gang Li,et al.  Hybrid Well-balanced WENO Schemes with Different Indicators for Shallow Water Equations , 2012, J. Sci. Comput..

[12]  Peng Li,et al.  Hybrid Compact-WENO Finite Difference Scheme For Detonation Waves Simulations , 2015 .

[13]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[14]  Guowei Yang,et al.  Hybrid finite compact‐WENO schemes for shock calculation , 2007 .

[15]  Sergio Pirozzoli,et al.  Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction , 2002 .

[16]  Alexandre Ern,et al.  A well‐balanced Runge–Kutta discontinuous Galerkin method for the shallow‐water equations with flooding and drying , 2008 .

[17]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[18]  Yulong Xing,et al.  A Survey of High Order Schemes for the Shallow Water Equations , 2014 .

[19]  Luka Sopta,et al.  Extension of ENO and WENO schemes to one-dimensional sediment transport equations , 2004 .

[20]  Yulong Xing,et al.  High order finite difference WENO schemes with the exact conservation property for the shallow water equations , 2005 .

[21]  P. Moin,et al.  A General Class of Commutative Filters for LES in Complex Geometries , 1998 .

[22]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[23]  Yuxin Ren,et al.  A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws , 2003 .

[24]  Wai-Sun Don,et al.  Multi-domain hybrid spectral-WENO methods for hyperbolic conservation laws , 2007, J. Comput. Phys..

[25]  Yulong Xing,et al.  Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations , 2010 .

[26]  Benedict D. Rogers,et al.  Mathematical balancing of flux gradient and source terms prior to using Roe's approximate Riemann solver , 2003 .

[27]  Gang Li,et al.  Hybrid weighted essentially non-oscillatory schemes with different indicators , 2010, J. Comput. Phys..

[28]  Luka Sopta,et al.  ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations , 2002 .

[29]  S. Balachandar,et al.  A Massively Parallel Multi-Block Hybrid Compact-WENO Scheme for Compressible Flows , 2009 .

[30]  Wai-Sun Don,et al.  High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws , 2011, J. Comput. Phys..

[31]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[32]  Luka Sopta,et al.  WENO schemes for balance laws with spatially varying flux , 2004 .