ORIGINAL PAPER

Enforcing Strong Stability of Explicit Runge-Kutta Methods with Superviscosity

Expand
  • 1 Department of Mathematics, The Ohio State University, Columbus, OH 43210, USA;
    2 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA

Received date: 2019-12-26

  Revised date: 2020-07-17

  Online published: 2021-11-25

Supported by

Research supported by NSF Grants DMS-1719410 and DMS-2010107, and by AFOSR Grant FA9550-20-1-0055.

Abstract

A time discretization method is called strongly stable (or monotone), if the norm of its numerical solution is nonincreasing. Although this property is desirable in various of contexts, many explicit Runge-Kutta (RK) methods may fail to preserve it. In this paper, we enforce strong stability by modifying the method with superviscosity, which is a numerical technique commonly used in spectral methods. Our main focus is on strong stability under the inner-product norm for linear problems with possibly non-normal operators. We propose two approaches for stabilization: the modifed method and the fltering method. The modifed method is achieved by modifying the semi-negative operator with a high order superviscosity term; the fltering method is to post-process the solution by solving a difusive or dispersive problem with small superviscosity. For linear problems, most explicit RK methods can be stabilized with either approach without accuracy degeneration. Furthermore, we prove a sharp bound (up to an equal sign) on difusive superviscosity for ensuring strong stability. For nonlinear problems, a fltering method is investigated. Numerical examples with linear non-normal ordinary diferential equation systems and for discontinuous Galerkin approximations of conservation laws are performed to validate our analysis and to test the performance.

Cite this article

Zheng Sun, Chi-Wang Shu . Enforcing Strong Stability of Explicit Runge-Kutta Methods with Superviscosity[J]. Communications on Applied Mathematics and Computation, 2021 , 3(4) : 671 -700 . DOI: 10.1007/s42967-020-00098-y

References

1. Abgrall, R.: A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes. J. Comput. Phys. 372, 640–666 (2018)
2. Butcher, J.C.: Numerical Methods for Ordinary Diferential Equations. John Wiley and Sons, Ltd., Chichester (2016)
3. Chen, G.Q., Du, Q., Tadmor, E.: Spectral viscosity approximations to multidimensional scalar conservation laws. Math. Comput. 61(204), 629–643 (1993)
4. Chen, T., Shu, C.-W.: Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws. J. Comput. Phys. 345, 427–461 (2017)
5. Chen, T., Shu, C.-W.: Review of entropy stable discontinuous Galerkin methods for systems of conservation laws on unstructured simplex meshes. CSIAM Trans. Appl. Math. 1, 1–52 (2020)
6. Cockburn, B., Hou, S., Shu, C.-W.: The Runge-Kutta local projection discontinuous Galerkin fnite element method for conservation laws. IV: the multidimensional case. Math. Comput. 54(190), 545–581 (1990)
7. Cockburn, B., Lin, S.-Y., Shu, C.-W.: TVB Runge-Kutta local projection discontinuous Galerkin fnite element method for conservation laws Ⅲ: one-dimensional systems. J. Comput. Phys. 84(1), 90–113 (1989)
8. Cockburn, B., Shu, C.-W.: TVB Runge-Kutta local projection discontinuous Galerkin fnite element method for conservation laws Ⅱ: general framework. Math. Comput. 52(186), 411–435 (1989)
9. Cockburn, B., Shu, C.-W.: The Runge-Kutta local projection P1-discontinuous-Galerkin fnite element method for scalar conservation laws. ESAIM Math. Model. Numer. Anal. 25(3), 337–361 (1991)
10. Cockburn, B., Shu, C.-W.: The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems. J. Comput. Phys. 141(2), 199–224 (1998)
11. Cockburn, B., Shu, C.-W.: Runge-Kutta discontinuous Galerkin methods for convection-dominated problems. J. Sci. Comput. 16(3), 173–261 (2001)
12. Fehlberg, E.: Klassische Runge-Kutta-formeln vierter und niedrigerer ordnung mit schrittweiten-kontrolle und ihre anwendung auf waermeleitungsprobleme. Computing 6(1/2), 61–71 (1970)
13. Fisher, T.C., Carpenter, M.H.: High-order entropy stable fnite diference schemes for nonlinear conservation laws: fnite domains. J. Comput. Phys. 252, 518–557 (2013)
14. Fjordholm, U.S., Mishra, S., Tadmor, E.: Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM J. Numer. Anal. 50(2), 544–573 (2012)
15. Gottlieb, D., Hesthaven, J.S.: Spectral methods for hyperbolic problems. J. Comput. Appl. Math. 128(1/2), 83–131 (2001)
16. Gottlieb, S., Shu, C.-W., Tadmor, E.: Strong stability-preserving high-order time discretization methods. SIAM Rev. 43(1), 89–112 (2001)
17. Guermond, J.-L., Prudhomme, S.: Mathematical analysis of a spectral hyperviscosity LES model for the simulation of turbulent fows. ESAIM: Math. Model. Numer. Anal. 37(6), 893–908 (2003)
18. Guo, H., Yang, Y.: Bound-preserving discontinuous Galerkin method for compressible miscible displacement in porous media. SIAM J. Sci. Comput. 39(5), A1969–A1990 (2017)
19. Gustafsson, B., Kreiss, H.-O., Oliger, J.: Time Dependent Problems and Diference Methods. John Wiley and Sons, Inc., New Jersey (1995)
20. Hesthaven, J.S., Gottlieb, S., Gottlieb, D.: Spectral Methods for Time-Dependent Problems. Cambridge University Press, Cambridge (2007)
21. Iserles, A.: A First Course in the Numerical Analysis of Diferential Equations. Cambridge University Press, Cambridge (2009)
22. Karamanos, G., Karniadakis, G.E.: A spectral vanishing viscosity method for large-eddy simulations. J. Comput. Phys. 163(1), 22–50 (2000)
23. Ketcheson, D.I.: Relaxation Runge-Kutta methods: conservation and stability for inner-product norms. SIAM J. Numer. Anal. 57(6), 2850–2870 (2019)
24. Kraaijevanger, J.F.B.M.: Contractivity of Runge-Kutta methods. BIT Numeri. Math. 31(3), 482–528 (1991)
25. Levy, D., Tadmor, E.: From semidiscrete to fully discrete: stability of Runge-Kutta schemes by the energy method. SIAM Rev. 40(1), 40–73 (1998)
26. Lozano, C.: Entropy production by explicit Runge-Kutta schemes. J. Sci. Comput. 76(1), 521–564 (2018)
27. Lozano, C.: Entropy production by implicit Runge-Kutta schemes. J. Sci. Comput. 79(3), 1832–1853 (2019)
28. Maday, Y., Kaber, S.M.O., Tadmor, E.: Legendre pseudospectral viscosity method for nonlinear conservation laws. SIAM J. Numer. Anal. 30(2), 321–342 (1993)
29. Maday, Y., Tadmor, E.: Analysis of the spectral vanishing viscosity method for periodic conservation laws. SIAM J. Numer. Anal. 26(4), 854–870 (1989)
30. Öfner, P., Glaubitz, J., Ranocha, H.: Analysis of artifcial dissipation of explicit and implicit time-integration methods. Int. J. Numer. Anal. Model. 17(3), 332–349 (2020)
31. Passot, T., Pouquet, A.: Hyperviscosity for compressible fows using spectral methods. J. Comput. Phys. 75(2), 300–313 (1988)
32. Qin, T., Shu, C.-W., Yang, Y.: Bound-preserving discontinuous Galerkin methods for relativistic hydrodynamics. J. Comput. Phys. 315, 323–347 (2016)
33. Ranocha, H.: On strong stability of explicit Runge-Kutta methods for nonlinear semibounded operators. IMA J. Numer. Anal. (2020)
34. Ranocha, H., Öfner, P.: L2 stability of explicit Runge-Kutta schemes. J. Sci. Comput. 75(2), 1040–1056 (2018)
35. Ranocha, H., Sayyari, M., Dalcin, L., Parsani, M., Ketcheson, D.I.: Relaxation Runge-Kutta methods: fully-discrete explicit entropy-stable schemes for the Euler and Navier-Stokes equations. SIAM J. Sci. Comput. 42(2), A612–A638 (2020)
36. Ranocha, H., Ketcheson, D.I.: Energy stability of explicit Runge-Kutta methods for non-autonomous or nonlinear problems. SIAM J. Numer. Anal. 58(6), 3382–3405 (2020)
37. Spijker, M.: Contractivity in the numerical solution of initial value problems. Numerische Mathematik 42(3), 271–290 (1983)
38. Sun, Z., Carrillo, J.A., Shu, C.-W.: A discontinuous Galerkin method for nonlinear parabolic equations and gradient fow problems with interaction potentials. J. Comput. Phys. 352, 76–104 (2018)
39. Sun, Z., Carrillo, J.A., Shu, C.-W.: An entropy stable high-order discontinuous Galerkin method for crossdifusion gradient fow systems. Kinetic Related Models 12(4), 885–908 (2019)
40. Sun, Z., Shu, C.-W.: Stability analysis and error estimates of Lax-Wendrof discontinuous Galerkin methods for linear conservation laws. ESAIM Math. Model. Numeri. Anal. 51(3), 1063–1087 (2017)
41. Sun, Z., Shu, C.-W.: Stability of the fourth order Runge-Kutta method for time-dependent partial diferential equations. Ann. Math. Sci. Appl. 2(2), 255–284 (2017)
42. Sun, Z., Shu, C.-W.: Strong stability of explicit Runge-Kutta time discretizations. SIAM J. Numer. Anal. 57(3), 1158–1182 (2019)
43. Tadmor, E.: The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Math. Comput. 49(179), 91–103 (1987)
44. Tadmor, E.: Convergence of spectral methods for nonlinear conservation laws. SIAM J. Numer. Anal. 26(1), 30–44 (1989)
45. Tadmor, E.: Super viscosity and spectral approximations of nonlinear conservation laws. Numerical methods for fuid dynamics IV. In: Baines, M.J., Morton, K.W. (eds) Proceedings of the 1992 Conference on Numerical Methods for Fluid Dynamics, Clarendon Press, Oxford, 69–82 (1993)
46. Tadmor, E.: From semidiscrete to fully discrete: stability of Runge-Kutta schemes by the energy method. Ⅱ. In: Estep, D., Tavener, S. (eds.) Collected Lectures on the Preservation of Stability under Discretization, Lecture Notes from Colorado State University Conference, Fort Collins, CO, 2001 Proceedings in Applied Mathematics, SIAM, 109, 25–49 (2002)
47. Tadmor, E.: Burgers’ equation with vanishing hyper-viscosity. Commun. Math. Sci. 2(2), 317–324 (2004)
48. Tadmor, E., Waagan, K.: Adaptive spectral viscosity for hyperbolic conservation laws. SIAM J. Sci. Comput. 34(2), A993–A1009 (2012)
49. Wu, K., Tang, H.: High-order accurate physical-constraints-preserving fnite diference WENO schemes for special relativistic hydrodynamics. J. Comput. Phys. 298, 539–564 (2015)
50. Xing, Y., Zhang, X., Shu, C.-W.: Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations. Adv. Water Resour. 33(12), 1476–1493 (2010)
51. Xu, Y., Zhang, Q., Shu, C.-W., Wang, H.: The L2-norm stability analysis of Runge-Kutta discontinuous Galerkin methods for linear hyperbolic equations. SIAM J. Numer. Anal. 57(4), 1574–1601 (2019)
52. Xu, Y., Meng, X., Shu, C.-W., Zhang, Q.: Superconvergence analysis of the Runge-Kutta discontinuous Galerkin methods for a linear hyperbolic equation. J. Sci. Comput. 84, 23 (2020)
53. Zhang, Q., Gao, F.: A fully-discrete local discontinuous Galerkin method for convection-dominated Sobolev equation. J. Sci. Comput. 51, 107–134 (2012)
54. Zhang, Q., Shu, C.-W.: Stability analysis and a priori error estimates of the third order explicit RungeKutta discontinuous Galerkin method for scalar conservation laws. SIAM J. Numer. Anal. 48(3), 1038– 1063 (2010)
55. Zhang, X., Shu, C.-W.: On maximum-principle-satisfying high order schemes for scalar conservation laws. J. Comput. Phys. 229(9), 3091–3120 (2010)
56. Zhang, X., Shu, C.-W.: On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. J. Comput. Phys. 229(23), 8918–8934 (2010)
Options
Outlines

/