2D well-balanced augmented ADER schemes for the Shallow Water Equations with bed elevation and extension to the rotating frame

Abstract In this work, an arbitrary order augmented WENO-ADER scheme for the resolution of the 2D Shallow Water Equations (SWE) with geometric source term is presented and its application to other shallow water models involving non-geometric sources is explored. This scheme is based in the 1D Augmented Roe Linearized-ADER (ARL-ADER) scheme, presented by the authors in a previous work and motivated by a suitable compromise between accuracy and computational cost. It can be regarded as an arbitrary order version of the Augmented Roe solver, which accounts for the contribution of continuous and discontinuous geometric source terms at cell interfaces in the resolution of the Derivative Riemann Problem (DRP). The main novelty of this work is the extension of the ARL-ADER scheme to 2 dimensions, which involves the design of a particular procedure for the integration of the source term with arbitrary order that ensures an exact balance between flux fluctuations and sources. This procedure makes the scheme preserve equilibrium solutions with machine precision and capture the transient waves accurately. The scheme is applied to the SWE with bed variation and is extended to handle non-geometric source terms such as the Coriolis source term. When considering the SWE with bed variation and Coriolis, the most relevant equilibrium states are the still water at rest and the geostrophic equilibrium. The traditional well-balanced property is extended to satisfy the geostrophic equilibrium. This is achieved by means of a geometric reinterpretation of the Coriolis source term. By doing this, the formulation of the source terms is unified leading to a single geometric source regarded as an apparent topography. The numerical scheme is tested for a broad variety of situations, including some cases where the first order scheme ruins the solution.

[1]  G. Kesserwani,et al.  Application of a second‐order Runge–Kutta discontinuous Galerkin scheme for the shallow water equations with source terms , 2008 .

[2]  Michael Dumbser,et al.  Comparison of solvers for the generalized Riemann problem for hyperbolic systems with source terms , 2012, J. Comput. Phys..

[3]  Adrián Navas-Montilla,et al.  Asymptotically and exactly energy balanced augmented flux-ADER schemes with application to hyperbolic conservation laws with geometric source terms , 2016, J. Comput. Phys..

[4]  Michael Dumbser,et al.  Space–time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting , 2014, 1412.0081.

[5]  Daniel Y. Le Roux,et al.  Analysis of Numerically Induced Oscillations in 2D Finite-Element Shallow-Water Models Part I: Inertia-Gravity Waves , 2007, SIAM J. Sci. Comput..

[6]  Javier Murillo,et al.  Augmented Roe's approaches for Riemann problems including source terms: Definition of stability region with application to the shallow water equations with rigid and deformable bed , 2012 .

[7]  Eleuterio F. Toro,et al.  Towards Very High Order Godunov Schemes , 2001 .

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

[9]  Rupert Klein,et al.  Preservation of the Discrete Geostrophic Equilibrium in Shallow Water Flows , 2011 .

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

[11]  Eleuterio F. Toro,et al.  Solvers for the high-order Riemann problem for hyperbolic balance laws , 2008, J. Comput. Phys..

[12]  Pilar García-Navarro,et al.  Implicit schemes with large time step for non‐linear equations: application to river flow hydraulics , 2004 .

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

[14]  Philippe G. LeFloch,et al.  A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime , 2011, J. Comput. Phys..

[15]  Adrián Navas-Montilla,et al.  A comprehensive explanation and exercise of the source terms in hyperbolic systems using Roe type solutions. Application to the 1D-2D shallow water equations , 2016 .

[16]  Eitan Tadmor,et al.  Well-balanced and energy stable schemes for the shallow water equations with discontinuous topography , 2011, J. Comput. Phys..

[17]  B. Sanders,et al.  Discretization of Integral Equations Describing Flow in Nonprismatic Channels with Uneven Beds , 2003 .

[18]  Vladimir Zeitlin,et al.  Frontal geostrophic adjustment and nonlinear wave phenomena in one-dimensional rotating shallow water. Part 2. High-resolution numerical simulations , 2004, Journal of Fluid Mechanics.

[19]  Vladimir Zeitlin,et al.  A robust well-balanced scheme for multi-layer shallow water equations , 2010 .

[20]  Eleuterio F. Toro,et al.  Derivative Riemann solvers for systems of conservation laws and ADER methods , 2006, J. Comput. Phys..

[21]  Mária Lukácová-Medvid'ová,et al.  Well-balanced schemes for the shallow water equations with Coriolis forces , 2018, Numerische Mathematik.

[22]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[23]  Emmanuel Audusse,et al.  Conservative discretization of Coriolis force in a finite volume framework , 2009, J. Comput. Phys..

[24]  A. Leroux,et al.  A well-balanced numerical scheme for the approximation of the shallow-water equations with topography: the resonance phenomenon , 2004 .

[25]  Manuel Jesús Castro Díaz,et al.  Finite Volume Simulation of the Geostrophic Adjustment in a Rotating Shallow-Water System , 2008, SIAM J. Sci. Comput..

[26]  Eleuterio F. Toro,et al.  ADER schemes for the shallow water equations in channel with irregular bottom elevation , 2006, J. Comput. Phys..

[27]  Valerio Caleffi,et al.  Well balancing of the SWE schemes for moving-water steady flows , 2017, J. Comput. Phys..

[28]  Nicolas G. Wright,et al.  Simple and efficient solution of the shallow water equations with source terms , 2010 .

[29]  Eleuterio F. Toro,et al.  Implicit, semi-analytical solution of the generalized Riemann problem for stiff hyperbolic balance laws , 2015, J. Comput. Phys..

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

[31]  Michael Dumbser,et al.  Well-Balanced High-Order Centred Schemes for Non-Conservative Hyperbolic Systems. Applications to Shallow Water Equations with Fixed and Mobile Bed , 2009 .

[32]  P. García-Navarro,et al.  On numerical treatment of the source terms in the shallow water equations , 2000 .

[33]  Yulong Xing,et al.  Exactly well-balanced discontinuous Galerkin methods for the shallow water equations with moving water equilibrium , 2014, J. Comput. Phys..

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

[35]  Ashley F. Emery,et al.  An Evaluation of Several Differencing Methods for Inviscid Fluid Flow Problems , 1968 .

[36]  L. Polvani,et al.  Nonlinear geostrophic adjustment, cyclone/anticyclone asymmetry, and potential vorticity rearrangement , 2000 .

[37]  E. Toro,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes - V. Local time stepping and p-adaptivity , 2007 .

[38]  Eleuterio F. Toro,et al.  ADER scheme on unstructured meshes for shallow water: simulation of tsunami waves , 2012 .

[39]  Javier Murillo,et al.  Augmented versions of the HLL and HLLC Riemann solvers including source terms in one and two dimensions for shallow flow applications , 2012, J. Comput. Phys..

[40]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[41]  Luca Solari,et al.  Conservative Scheme for Numerical Modeling of Flow in Natural Geometry , 2008 .

[42]  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 .

[43]  Adrián Navas-Montilla,et al.  Energy balanced numerical schemes with very high order. The Augmented Roe Flux ADER scheme. Application to the shallow water equations , 2015, J. Comput. Phys..

[44]  Eleuterio F. Toro,et al.  Exact solution of the Riemann problem for the shallow water equations with discontinuous bottom geometry , 2008, J. Comput. Phys..

[45]  W. Melville,et al.  Kelvin Fronts on the Equatorial Thermocline , 2000 .

[46]  Eleuterio F. Toro,et al.  A sub-cell WENO reconstruction method for spatial derivatives in the ADER scheme , 2013, J. Comput. Phys..

[47]  Javier Murillo,et al.  Energy balance numerical schemes for shallow water equations with discontinuous topography , 2013, J. Comput. Phys..

[48]  Valerio Caleffi,et al.  A comparison between bottom-discontinuity numerical treatments in the DG framework , 2015 .

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

[50]  Jostein R. Natvig,et al.  High-order well-balanced finite-volume schemes for barotropic flows: Development and numerical comparisons , 2014, 1412.3609.

[51]  M. J. Castro,et al.  ADER schemes on unstructured meshes for nonconservative hyperbolic systems: Applications to geophysical flows , 2009 .

[52]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[53]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[54]  Manuel Jesús Castro Díaz,et al.  High order exactly well-balanced numerical methods for shallow water systems , 2013, J. Comput. Phys..

[55]  Numerical simulation of Rossby wave in shallow water , 2013 .

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

[57]  Javier Murillo,et al.  Accurate numerical modeling of 1D flow in channels with arbitrary shape. Application of the energy balanced property , 2014, J. Comput. Phys..

[58]  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 .

[59]  Michael Dumbser,et al.  Fast high order ADER schemes for linear hyperbolic equations , 2004 .

[60]  Eleuterio F. Toro,et al.  ADER schemes for scalar non-linear hyperbolic conservation laws with source terms in three-space dimensions , 2005 .

[61]  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..

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

[63]  M. Thanh,et al.  The Riemann problem for the shallow water equations with discontinuous topography , 2007, 0712.3778.

[64]  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..

[65]  E. Toro,et al.  Solution of the generalized Riemann problem for advection–reaction equations , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[66]  Yulong Xing,et al.  High-order well-balanced finite volume WENO schemes for shallow water equation with moving water , 2007, J. Comput. Phys..

[67]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .