Jacobian-free explicit multiderivative Runge-Kutta methods for hyperbolic conservation laws

Based on the recent development of Jacobian-free Lax-Wendroff (LW) approaches for solving hyperbolic conservation laws [Zorio, Baeza and Mulet, Journal of Scientific Computing 71:246-273, 2017], [Carrillo and Parés, Journal of Scientific Computing 80:1832-1866, 2019], a novel collection of explicit Jacobian-free multistage multiderivative solvers for hyperbolic conservation laws is presented in this work. In contrast to Taylor time-integration methods, multiderivative RungeKutta (MDRK) techniques achieve higher-order of consistency not only through the excessive addition of higher temporal derivatives, but also through the addition of Runge-Kutta-type stages. This adds more flexibility to the time integration in such a way that more stable and more efficient schemes could be identified. The novel method permits the practical application of MDRK schemes. In their original form, they are difficult to utilize as higher-order flux derivatives have to be computed analytically. Here we overcome this by adopting a Jacobian-free approximation of those derivatives. In this paper, we analyze the novel method with respect to order of consistency and stability. We show that the linear CFL number varies significantly with the number of derivatives used. Results are verified numerically on several representative testcases.

[1]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[2]  Angelika Fruehauf,et al.  A First Course In Numerical Analysis , 2016 .

[3]  Tsuyoshi Murata,et al.  {m , 1934, ACML.

[4]  Michael Dumbser,et al.  The discontinuous Galerkin method with Lax-Wendroff type time discretizations , 2005 .

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

[6]  Alexander Jaust,et al.  Implicit Multiderivative Collocation Solvers for Linear Partial Differential Equations with Discontinuous Galerkin Spatial Discretizations , 2017, Journal of Scientific Computing.

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

[8]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[9]  Alexander Jaust,et al.  Implicit Multistage Two-Derivative Discontinuous Galerkin Schemes for Viscous Conservation Laws , 2015, J. Sci. Comput..

[10]  Oleg A. Yakimenko,et al.  Symbolic Math Toolbox , 2010 .

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

[12]  Wei Guo,et al.  A New Lax–Wendroff Discontinuous Galerkin Method with Superconvergence , 2015, J. Sci. Comput..

[13]  Chi-Wang Shu,et al.  Finite Difference WENO Schemes with Lax-Wendroff-Type Time Discretizations , 2002, SIAM J. Sci. Comput..

[14]  Zachary Grant,et al.  Explicit Strong Stability Preserving Multistage Two-Derivative Time-Stepping Schemes , 2015, J. Sci. Comput..

[15]  With Invariant Submanifolds,et al.  Systems of Conservation Laws , 2009 .

[16]  Zachary J. Grant,et al.  HIGH ORDER UNCONDITIONALLY STRONG STABILITY PRESERVING MULTI-DERIVATIVE IMPLICIT AND IMEX RUNGE–KUTTA METHODS WITH ASYMPTOTIC PRESERVING PROPERTIES , 2021 .

[17]  Eleuterio F. Toro,et al.  ADER schemes for three-dimensional non-linear hyperbolic systems , 2005 .

[18]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[19]  Jianxian Qiu,et al.  Development and Comparison of Numerical Fluxes for LWDG Methods , 2008 .

[20]  Michael Dumbser,et al.  Efficient Implementation of ADER Discontinuous Galerkin Schemes for a Scalable Hyperbolic PDE Engine , 2018, Axioms.

[21]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations, Second Edition , 2004 .

[22]  William E. Schiesser,et al.  Linear and nonlinear waves , 2009, Scholarpedia.

[23]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[24]  Editors , 2003 .

[25]  Turgut Özis,et al.  Derivation of three-derivative Runge-Kutta methods , 2016, Numerical Algorithms.

[26]  Jiequan Li,et al.  A Two-Stage Fourth Order Time-Accurate Discretization for Lax-Wendroff Type Flow Solvers I. Hyperbolic Conservation Laws , 2015, SIAM J. Sci. Comput..

[27]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[28]  Andrew J. Christlieb,et al.  High-Order Multiderivative Time Integrators for Hyperbolic Conservation Laws , 2013, J. Sci. Comput..

[29]  高等学校計算数学学報編輯委員会編.  高等学校計算数学学報 = Numerical mathematics , 1979 .

[30]  M. A. Wolfe.  A first course in numerical analysis , 1972 .

[31]  Robert P. K. Chan,et al.  On explicit two-derivative Runge-Kutta methods , 2010, Numerical Algorithms.

[32]  Jianxian Qiu,et al.  Simulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes , 2011, J. Sci. Comput..

[33]  Claus-Dieter Munz,et al.  ADER: A High-Order Approach for Linear Hyperbolic Systems in 2D , 2002, J. Sci. Comput..

[34]  P. Mulet,et al.  On approximate implicit Taylor methods for ordinary differential equations , 2020, Computational and Applied Mathematics.

[35]  Antonio Baeza,et al.  An Approximate Lax–Wendroff-Type Procedure for High Order Accurate Schemes for Hyperbolic Conservation Laws , 2017, J. Sci. Comput..

[36]  Hugo Carrillo,et al.  Compact Approximate Taylor Methods for Systems of Conservation Laws , 2019, Journal of Scientific Computing.

[37]  David Zorío,et al.  Lax-Wendroff Approximate Taylor Methods with Fast and Optimized Weighted Essentially Non-oscillatory Reconstructions , 2020, J. Sci. Comput..