ORIGINAL PAPERS

Quinpi: Integrating Conservation Laws with CWENO Implicit Methods

Expand
  • 1 Dipartimento di Matematica, Università La Sapienza, P. le Aldo Moro, 5, 00185 Roma, Italy;
    2 Dipartimento di Scienza e Alta Tecnologia, Università dell'Insubria, Via Valleggio, 11, 22100 Como, Italy

Received date: 2021-01-31

  Revised date: 2021-08-02

  Online published: 2023-03-08

Supported by

This work was partly supported by MIUR (Ministry of University and Research) PRIN2017 project number 2017KKJP4X, and Progetto di Ateneo Sapienza, number RM120172B41DBF3A.

Abstract

Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration is quite straightforward for first-order schemes. High order schemes instead also need to control spurious oscillations, which requires limiting in space and time also in the linear case. We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor, which is used both to set up the nonlinear weights of a standard high order space reconstruction, and to achieve limiting in time. In this preliminary work, we concentrate on the case of a third-order scheme, based on diagonally implicit Runge Kutta (DIRK) integration in time and central weighted essentially non-oscillatory (CWENO) reconstruction in space. The numerical tests involve linear and nonlinear scalar conservation laws.

Cite this article

G. Puppo, M. Semplice, G. Visconti . Quinpi: Integrating Conservation Laws with CWENO Implicit Methods[J]. Communications on Applied Mathematics and Computation, 2023 , 5(1) : 343 -369 . DOI: 10.1007/s42967-021-00171-0

References

1. Abbate, E., Iollo, A., Puppo, G.: An all-speed relaxation scheme for gases and compressible materials. J. Comput. Phys. 351, 1–24 (2017)
2. Alexander, R.: Diagonally implicit Runge-Kutta methods for stiff O.D.E.’s. SIAM J. Numer. Anal. 14(6), 1006–1021 (1977)
3. Arbogast, T., Huang, C., Zhao, X., King, D.N.: A third order, implicit, finite volume, adaptive Runge-Kutta WENO scheme for advection-diffusion equations. Comput. Methods Appl. Mech. Eng. (2020). https:// doi. org/ 10. 1016/j. cma. 2020. 113155
4. Balsara, D.S., Garain, S., Shu, C.-W.: An efficient class of WENO schemes with adaptive order. J. Comput. Phys. 326, 780–804 (2016)
5. Berger, M.J., LeVeque, R.J.: Adaptive mesh refinement using wave-propagation algorithms for hyperbolic systems. SIAM J. Numer. Anal. 35(6), 2298–2316 (1998). https:// doi. org/ 10. 1137/ S0036 14299 73159 74
6. Boscarino, S., Cho, S., Russo, G., Yun, S.: High order conservative semi-Lagrangian scheme for the BGK model of the Boltzmann equation. Commun. Comput. Phys. 29, 1–56 (2021)
7. Boscarino, S., Russo, G., Scandurra, L.: All Mach number second order semi-implicit scheme for the Euler equations of gas dynamics. J. Sci. Comput. 77(2), 850–884 (2018)
8. Carlini, E., Ferretti, R., Russo, G.: A weighted essentially nonoscillatory, large time-step scheme for Hamilton-Jacobi equations. SIAM J. Sci. Comput. 27, 1071–1091 (2005)
9. Castro-Dìaz, M.J., Semplice, M.: Third-and fourth-order well-balanced schemes for the shallow water equations based on the CWENO reconstruction. Int. J. Numer. Meth. Fluid 89(8), 304–325 (2019). https:// doi. org/ 10. 1002/ fld. 4700
10. Cravero, I., Puppo, G., Semplice, M., Visconti, G.: CWENO: uniformly accurate reconstructions for balance laws. Math. Comput. 87(312), 1689–1719 (2018). https:// doi. org/ 10. 1090/ mcom/ 3273
11. Cravero, I., Semplice, M.: On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes. J. Sci. Comput. 67, 1219–1246 (2016). https:// doi. org/ 10. 1007/ s10915-015-0123-3
12. Cravero, I., Semplice, M., Visconti, G.: Optimal definition of the nonlinear weights in multidimensional central WENOZ reconstructions. SIAM J. Numer. Anal. 57(5), 2328–2358 (2019). https:// doi. org/ 10. 1007/ s10915-015-0123-3
13. Degond, P., Tang, M.: All speed scheme for the low Mach number limit of the isentropic Euler equations. Commun. Comput. Phys. 10(1), 1–31 (2011)
14. Dellacherie, S.: Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number. J. Comput. Phys. 229(4), 978–1016 (2010)
15. Dimarco, G., Loubère, R., Vignal, M.H.: Study of a new asymptotic preserving scheme for the Euler system in the low Mach number limit. SIAM J. Sci. Comput. 39(5), A2099–A2128 (2017)
16. Dimarco, G., Loubère, R., Michel-Dansac, V., Vignal, M.H.: Second-order implicit-explicit total variation diminishing schemes for the Euler system in the low Mach regime. J. Comput. Phys. 372, 178–201 (2018)
17. Dimarco, G., Pareschi, L.: Numerical methods for kinetic equations. Acta Numer. 23, 369–520 (2014). https:// doi. org/ 10. 1017/ S0962 49291 40000 63
18. Dumbser, M., Balsara, D.S., Toro, E.F., Munz, C.D.: A unified framework for the construction of onestep finite volume and discontinuous Galerkin schemes on unstructured meshes. J. Comput. Phys. 227, 8209–8253 (2008). https:// doi. org/ 10. 1016/j. jcp. 2008. 05. 025
19. Dumbser, M., Boscheri, W., Semplice, M., Russo, G.: Central weighted ENO schemes for hyperbolic conservation laws on fixed and moving unstructered meshes. SIAM J. Sci. Comput. 39(6), A2564-A2591 (2017). https:// doi. org/ 10. 1137/ 17M11 11036
20. Duraisamy, K., Baeder, D.: Implicit scheme for hyperbolic conservation laws using non-oscillatory reconstruction in space and time. SIAM J. Sci. Comput. 29, 2607–2620 (2007)
21. Duraisamy, K., Baeder, D., Liu, J.: Concepts and application of time-limiters to high resolution schemes. J. Sci. Comput. 19, 139–162 (2003)
22. Forth, S.A.: A second order accurate space-time limited BDF scheme for the linear advection equation. In: Godunov Methods, pp. 335–342. Springer, Berlin (2001)
23. Gottlieb, S., Mullen, J.S., Ruuth, S.J.: A fifth order flux implicit WENO method. J. Sci. Comput. 27, 271–287 (2006). https:// doi. org/ 10. 1007/ s10915-005-9034-z
24. Gottlieb, S., Shu, C.-W., Tadmor, E.: Strong stability preserving high-order time discretization methods. SIAM Rev. 43, 73–85 (2001)
25. Groppi, M., Russo, G., Stracquadanio, G.: High order semi-Lagrangian methods for the BGK equation. Commun. Math. Sci. 14, 389–414 (2016)
26. Harten, A.: High resolution schemes for hyperbolic conservation laws. J. Comput. Phys. 49(3), 357-393 (1983). https:// doi. org/ 10. 1016/ 0021-9991(83) 90136-5
27. Harten, A.: On a class of high resolution total-variation-stable finite-difference schemes. SIAM J. Numer. Anal. 21, 1–23 (1984)
28. Jiang, G.S., Shu, C.-W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1996)
29. Jin, S., Xin, Z.: The relaxation schemes for systems of conservation laws in arbitrary space dimensions. Commun. Pure Appl. Math. 48(3), 235–276 (1995)
30. Ketcheson, D.I., MacDonald, C.B., Ruuth, S.J.: Spatially partitioned embedded Runge-Kutta methods. SIAM J. Numer. Anal. 51(5), 2887–2910 (2013). https:// doi. org/ 10. 1137/ 13090 6258
31. Kolb, O.: On the full and global accuracy of a compact third order WENO scheme. SIAM J. Numer. Anal. 52(5), 2335–2355 (2014). https:// doi. org/ 10. 1137/ 13094 7568
32. LeVeque, R.: Finite volume methods for hyperbolic problems. In: Cambridge Texts in Applied Mathematics. Cambridge University Press, Cambridge (2004)
33. Lemou, M., Mieussens, L.: A new asymptotic preserving scheme based on micro-macro formulation for linear kinetic equations in the diffusion limit. SIAM J. Sci. Comput. 31, 334–368 (2008)
34. Levy, D., Puppo, G., Russo, G.: Compact central WENO schemes for multidimensional conservation laws. SIAM J. Sci. Comput. 22(2), 656–672 (2000). https:// doi. org/ 10. 1137/ S1064 82759 93594 61
35. Nørsett, S.P., Wanner, G.: Perturbed collocation and Runge-Kutta methods. Numer. Math. 38, 193–208 (1981). https:// doi. org/ 10. 1007/ BF013 97089
36. Pieraccini, S., Puppo, G.: Microscopically implicit-macroscopically explicit schemes for the BGK equation. J. Comput. Phys. 231, 299–327 (2012)
37. Semplice, M., Coco, A., Russo, G.: Adaptive mesh refinement for hyperbolic systems based on thirdorder compact WENO reconstruction. J. Sci. Comput. 66, 692–724 (2016). https:// doi. org/ 10. 1007/ s10915-015-0038-z
38. Semplice, M., Visconti, G.: Efficient implementation of adaptive order reconstructions. J. Sci. Comput. 83, 6 (2020). https:// doi. org/ 10. 1007/ s10915-020-01156-6
39. Shu, C.-W.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations (Cetraro, 1997), Lecture Notes in Math., vol. 1697, pp. 325–432. Springer, Berlin (1998)
40. Sweby, P.: High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Anal. 21(5), 995–1011 (1984). https:// doi. org/ 10. 1137/ 07210 62
41. Tavelli, M., Dumbser, M.: A pressure-based semi-implicit space-time discontinuous Galerkin method on staggered unstructured meshes for the solution of the compressible Navier-Stokes equations at all Mach numbers. J. Comput. Phys. 341, 341–376 (2017). https:// doi. org/ 10. 1016/j. jcp. 2017. 03. 030
42. Zennaro, M.: Natural continuous extensions of Runge-Kutta methods. Math. Comput. 46, 119–133 (1986)
Options
Outlines

/