Flux-gradient and source-term balancing for certain high resolution shock-capturing schemes

Abstract We present an extension of Marquina’s flux formula, as introduced in Fedkiw et al. [Fedkiw RP, Merriman B, Donat R, Osher S. The penultimate scheme for systems of conservation laws: finite difference ENO with Marquina’s flux splitting. In: Hafez M, editor. Progress in numerical solutions of partial differential equations, Arcachon, France; July 1998], for the shallow water system. We show that the use of two different Jacobians at cell interfaces prevents the scheme from satisfying the exact C -property [Bermudez A, Vazquez ME. Upwind methods for hyperbolic conservation laws with source terms. Comput Fluids 1994;23(8):1049–71] while the approximate C -property is satisfied for higher order versions of the scheme. The use of a single Jacobian in Marquina’s flux splitting formula leads to a numerical scheme satisfying the exact C -property, hence we propose a combined technique that uses Marquina’s two sided decomposition when the two adjacent states are not close and a single decomposition otherwise. Finally, we propose a special treatment at wet/dry fronts and situations of dry bed generation.

[1]  E. Toro.  Shock-Capturing Methods for Free-Surface Shallow Flows , 2001 .

[2]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[3]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[4]  Luka Sopta,et al.  Balanced finite volume WENO and central WENO schemes for the shallow water and the open-channel flow equations , 2004 .

[5]  José M. Corberán,et al.  Construction of second-order TVD schemes for nonhomogeneous hyperbolic conservation laws , 2001 .

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

[7]  Pilar García-Navarro,et al.  A numerical model for the flooding and drying of irregular domains , 2002 .

[8]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[9]  M. Berzins,et al.  An unstructured finite-volume algorithm for predicting flow in rivers and estuaries , 1998 .

[10]  Theodore E. Simos,et al.  A two-step method with phase-lag of order infinity for the numerical integration of second order periodic initial-value problem , 1991, Int. J. Comput. Math..

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

[12]  Graziano Guerra.  Well-posedness for a scalar conservation law with singular nonconservative source , 2004 .

[13]  Antonio Marquina,et al.  Capturing Shock Reflections , 1996 .

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

[15]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[16]  Ronald Fedkiw,et al.  The Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina’s Flux Splitting1 , 2010 .

[18]  On different flux splittings and flux functions in WENO schemes for balance laws , 2006 .

[19]  Philip L. Roe,et al.  Upwind differencing schemes for hyperbolic conservation laws with source terms , 1987 .

[20]  Yulong Xing,et al.  High order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source terms , 2006, J. Comput. Phys..

[21]  Pilar García-Navarro,et al.  Unsteady free surface flow simulation over complex topography with a multidimensional upwind technique , 2003 .

[22]  J. Greenberg,et al.  Analysis and Approximation of Conservation Laws with Source Terms , 1997 .

[23]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[24]  Laurent Gosse,et al.  A Well-Balanced Scheme Using Non-Conservative Products Designed for Hyperbolic Systems of Conservati , 2001 .

[25]  A. I. Delis,et al.  Relaxation schemes for the shallow water equations , 2003 .

[26]  L. Gosse.  A well-balanced flux-vector splitting scheme designed for hyperbolic systems of conservation laws with source terms☆ , 2000 .

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

[28]  D. Zhao,et al.  Finite‐Volume Two‐Dimensional Unsteady‐Flow Model for River Basins , 1994 .

[29]  Laurent Gosse,et al.  Localization effects and measure source terms in numerical schemes for balance laws , 2002, Math. Comput..

[30]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[31]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

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