Multi-moment ADER-Taylor methods for systems of conservation laws with source terms in one dimension

A new integration method combining the ADER time discretization with a multi-moment finite-volume framework is introduced. ADER runtime is reduced by performing only one Cauchy-Kowalewski (C-K) procedure per cell per time step and by using the Differential Transform Method for high-order derivatives. Three methods are implemented: (1) single-moment WENO (WENO), (2) two-moment Hermite WENO (HWENO), and (3) entirely local multi-moment (MM-Loc). MM-Loc evolves all moments, sharing the locality of Galerkin methods yet with a constant time step during p-refinement.Five 1-D experiments validate the methods: (1) linear advection, (2) Burger's equation shock, (3) transient shallow-water (SW), (4) steady-state SW simulation, and (5) SW shock. WENO and HWENO methods showed expected polynomial h-refinement convergence and successfully limited oscillations for shock experiments. MM-Loc showed expected polynomial h-refinement and exponential p-refinement convergence for linear advection and showed sub-exponential (yet super-polynomial) convergence with p-refinement in the SW case.HWENO accuracy was generally equal to or better than a five-moment MM-Loc scheme. MM-Loc was less accurate than RKDG at lower refinements, but with greater h- and p-convergence, RKDG accuracy is eventually surpassed. The ADER time integrator of MM-Loc also proved more accurate with p-refinement at a CFL of unity than a semi-discrete RK analog of MM-Loc. Being faster in serial and requiring less frequent inter-node communication than Galerkin methods, the ADER-based MM-Loc and HWENO schemes can be spatially refined and have the same runtime, making them a competitive option for further investigation.

[1]  Michael Dumbser,et al.  A sub-cell based indicator for troubled zones in RKDG schemes and a novel class of hybrid RKDG+HWENO schemes , 2004, J. Comput. Phys..

[2]  Guy Capdeville,et al.  A Hermite upwind WENO scheme for solving hyperbolic conservation laws , 2008, J. Comput. Phys..

[3]  W. Collins,et al.  Description of the NCAR Community Atmosphere Model (CAM 3.0) , 2004 .

[4]  Eleuterio F. Toro,et al.  ADER: Arbitrary High Order Godunov Approach , 2002, J. Sci. Comput..

[5]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Transport Scheme on the Cubed Sphere , 2005 .

[6]  G. Russo,et al.  Central WENO schemes for hyperbolic systems of conservation laws , 1999 .

[7]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method II: Two dimensional case , 2005 .

[8]  Liu Wei,et al.  A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems , 2012, J. Comput. Phys..

[9]  Gabriella Puppo,et al.  A Fourth-Order Central WENO Scheme for Multidimensional Hyperbolic Systems of Conservation Laws , 2002, SIAM J. Sci. Comput..

[10]  Chi-Wang Shu,et al.  High Order ENO and WENO Schemes for Computational Fluid Dynamics , 1999 .

[11]  W. Pankiewicz.  Algorithms: Algorithm 337: calculation of a polynomial and its derivative values by Horner scheme , 1968, CACM.

[12]  Gerald Warnecke,et al.  Finite volume evolution Galerkin methods for nonlinear hyperbolic systems , 2002 .

[13]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

[14]  Feng Xiao,et al.  CIP/multi-moment finite volume method for Euler equations: A semi-Lagrangian characteristic formulation , 2007, J. Comput. Phys..

[15]  Guy Capdeville,et al.  A central WENO scheme for solving hyperbolic conservation laws on non-uniform meshes , 2008, J. Comput. Phys..

[16]  Clint Dawson,et al.  Semi discrete discontinuous Galerkin methods and stage-exceeding-order, strong-stability-preserving Runge-Kutta time discretizations , 2007, J. Comput. Phys..

[17]  Aydin Kurnaz,et al.  n-Dimensional differential transformation method for solving PDEs , 2005, Int. J. Comput. Math..

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

[19]  A. Gofen.  The ordinary differential equations and automatic differentiation unified , 2009 .

[20]  Cha'o-Kuang Chen,et al.  Application of differential transformation to eigenvalue problems , 1996 .

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

[22]  Michael Dumbser,et al.  ADER discontinuous Galerkin schemes for aeroacoustics , 2005 .

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