High order schemes for hyperbolic problems using globally continuous approximation and avoiding mass matrices

Abgrall, Rémi (2017). High order schemes for hyperbolic problems using globally continuous approximation and avoiding mass matrices. Journal of Scientific Computing, 73(2-3):461-494.

Abstract

When integrating unsteady problems using globally continuous representation of the solution, as for continuous finite element methods, one faces the problem of inverting a mass matrix. In some cases, one has to recompute this mass matrix at each time steps. In some other methods that are not directly formulated by standard variational principles, it is not clear how to write an invertible mass matrix. Hence, in this paper, we show how to avoid this problem for hyperbolic systems, and we also detail the conditions under which this is possible. Analysis and simulation support our conclusions, namely that it is possible to avoid inverting mass matrices without sacrificing the accuracy of the scheme. This paper is an extension of Abgrall et al. (in: Karasözen B, Manguoglu M, Tezer-Sezgin M, Goktepe S, Ugur O (eds) Numerical mathematics and advanced applications ENUMATH 2015. Lecture notes in computational sciences and engineering, vol 112, Springer, Berlin, 2016) and Ricchiuto and Abgrall (J Comput Phys 229(16):5653–5691, 2010). © 2017 Springer Science+Business Media, LLC

Abstract

When integrating unsteady problems using globally continuous representation of the solution, as for continuous finite element methods, one faces the problem of inverting a mass matrix. In some cases, one has to recompute this mass matrix at each time steps. In some other methods that are not directly formulated by standard variational principles, it is not clear how to write an invertible mass matrix. Hence, in this paper, we show how to avoid this problem for hyperbolic systems, and we also detail the conditions under which this is possible. Analysis and simulation support our conclusions, namely that it is possible to avoid inverting mass matrices without sacrificing the accuracy of the scheme. This paper is an extension of Abgrall et al. (in: Karasözen B, Manguoglu M, Tezer-Sezgin M, Goktepe S, Ugur O (eds) Numerical mathematics and advanced applications ENUMATH 2015. Lecture notes in computational sciences and engineering, vol 112, Springer, Berlin, 2016) and Ricchiuto and Abgrall (J Comput Phys 229(16):5653–5691, 2010). © 2017 Springer Science+Business Media, LLC

Statistics

Citations

Dimensions.ai Metrics
25 citations in Web of Science®
30 citations in Scopus®

Altmetrics

Detailed statistics

Item Type: Journal Article, refereed, original work 07 Faculty of Science > Institute of Mathematics 510 Mathematics Physical Sciences > Software Physical Sciences > Theoretical Computer Science Physical Sciences > Numerical Analysis Physical Sciences > General Engineering Physical Sciences > Computational Theory and Mathematics Physical Sciences > Computational Mathematics Physical Sciences > Applied Mathematics English 18 July 2017 09 Oct 2017 16:38 26 Jan 2022 13:41 Springer 0885-7474 Closed https://doi.org/10.1007/s10915-017-0498-4 : FunderSNSF: Grant ID200021_153604: Project TitleHigh fidelity simulation for compressible material