A large time step 1D upwind explicit scheme (CFL > 1): Application to shallow water equations

It is possible to relax the Courant-Friedrichs-Lewy condition over the time step when using explicit schemes. This method, proposed by Leveque, provides accurate and correct solutions of non-sonic shocks. Rarefactions need some adjustments which are explored in the present work with scalar equation and systems of equations. The non-conservative terms that appear in systems of conservation laws introduce an extra difficulty in practical application. The way to deal with source terms is incorporated into the proposed procedure. The boundary treatment is analysed and a reflection wave technique is considered. In presence of strong discontinuities or important source terms, a strategy is proposed to control the stability of the method allowing the largest time step possible. The performance of the above scheme is evaluated to solve the homogeneous shallow water equations and the shallow water equations with source terms.

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

[2]  Ami Harten,et al.  On a large time-step high resolution scheme , 1986 .

[3]  Pilar García-Navarro,et al.  Flux difference splitting and the balancing of source terms and flux gradients , 2000 .

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

[5]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[6]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[7]  Javier Murillo,et al.  Extension of an explicit finite volume method to large time steps (CFL>1): application to shallow water flows , 2006 .

[8]  Matthew R. Norman,et al.  A low communication and large time step explicit finite-volume solver for non-hydrostatic atmospheric dynamics , 2011, J. Comput. Phys..

[9]  David L. George,et al.  Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation , 2008, J. Comput. Phys..

[10]  M. Vázquez-Cendón.  Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .

[11]  Chun-Hian Lee,et al.  A class of large time step Godunov schemes for hyperbolic conservation laws and applications , 2011, J. Comput. Phys..

[12]  Pilar García-Navarro,et al.  Efficient construction of high‐resolution TVD conservative schemes for equations with source terms: application to shallow water flows , 2001 .

[13]  Randall J. LeVeque,et al.  Large time step shock-capturing techniques for scalar conservation laws , 1981 .

[14]  Nancy Nichols,et al.  Analytic Benchmark Solutions for Open-Channel Flows , 1997 .

[15]  Randall J. LeVeque,et al.  A large Time Step Generalization of Godunov’s Method for Systems of Conservation Laws , 1985 .

[16]  Javier Murillo,et al.  The influence of source terms on stability, accuracy and conservation in two‐dimensional shallow flow simulation using triangular finite volumes , 2007 .

[17]  Javier Murillo,et al.  Weak solutions for partial differential equations with source terms: Application to the shallow water equations , 2010, J. Comput. Phys..

[18]  Giorgio Rosatti,et al.  The Riemann Problem for the one-dimensional, free-surface Shallow Water Equations with a bed step: Theoretical analysis and numerical simulations , 2010, J. Comput. Phys..