zbMATH — the first resource for mathematics

On the efficient implementation of PVM methods and simple Riemann solvers. Application to the Roe method for large hyperbolic systems. (English) Zbl 07329288
Summary: Polynomial Viscosity Matrix (PVM) methods can be considered as approximations of the Roe method in which the absolute value of the Roe matrix appearing in the numerical viscosity is replaced by the evaluation of the Roe matrix at a chosen polynomial that approximates the absolute value function. They are in principle cheaper than the Roe method since the computation and the inversion of the eigenvector matrix is not necessary. In this article, an efficient implementation of the PVM based on polynomials that interpolate the absolute value function at some points is presented. This implementation is based on the Newton form of the polynomials. Moreover, many numerical methods based on simple Riemann solvers (SRS) may be interpreted as PVM methods and thus this implementation can be also applied to them: the close relation between PVM methods and simple Riemann solvers is revisited here and new shorter proofs based on the classical interpolation theory are given. In particular, Roe method can be interpreted both as a SRS and as a PVM method so that the new implementation can be used. This implementation, that avoids the computation and the inversion of the eigenvector matrix, is called Newton Roe method. Newton Roe method yields the same numerical results of the standard Roe method, with less runtime for large PDE systems. Numerical results for two-layer Shallow Water Equations and Quadrature-Based Moment Equations show a significant speedup if the number of equations is large enough.
65-XX Numerical analysis
15-XX Linear and multilinear algebra; matrix theory
Eigen; PVM ; HLLE
Full Text: DOI
[1] Roe, P., Approximate Riemann solvers, parameter vectors, and difference schemes, J. Comput. Phys., 43, 2, 357-372 (1981) · Zbl 0474.65066
[2] Toumi, I., A weak formulation of Roe’s approximate Riemann solver, J. Comput. Phys., 102, 2, 360-373 (1992) · Zbl 0783.65068
[3] Parés, C., Numerical methods for nonconservative hyperbolic systems: a theoretical framework., SIAM J. Numer. Anal., 44, 1, 300-321 (2006) · Zbl 1130.65089
[4] Harten, A.; Lax, P. D.; Leer, B. V., On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev., 25, 1, 35-61 (1983) · Zbl 0565.65051
[5] Castro, M.; Fernández-Nieto, E., A class of computationally fast first order finite volume solvers: PVM methods, SIAM J. Sci. Comput., 34, 4, A2173-A2196 (2012) · Zbl 1253.65167
[6] Morales de Luna, T.; Castro, M. J.; Parés, C., Relation between PVM schemes and simple Riemann solvers, Numer. Methods Partial Differ. Equ., 30, 4, 1315-1341 (2014) · Zbl 1297.65102
[7] Castro, M.; Macías, J.; Parés, C., A Q-scheme for a class of systems of coupled conservation laws with source term. Application to a two-layer 1-D shallow water system, ESAIM, 35, 1, 107-127 (2001) · Zbl 1094.76046
[8] Koellermeier, J., Derivation and Numerical Solution of Hyperbolic Moment Equations for Rarefied Gas Flows (2017), RWTH Aachen University: RWTH Aachen University Aachen, dissertation · Zbl 1408.76454
[9] Krvavica, N.; Tuhtan, M.; Jelenić, G., Analytical implementation of roe solver for two-layer shallow water equations with accurate treatment for loss of hyperbolicity, Adv. Water Resour., 122, 187-205 (2018)
[10] Parés, C.; Castro, M., On the well-balance property of Roe’s method for nonconservative hyperbolic systems. applications to shallow-water systems, ESAIM: M2AN, 38, 5, 821-852 (2004) · Zbl 1130.76325
[11] G. Guennebaud, B. Jacob, et al., Eigen v3, 2010, (http://eigen.tuxfamily.org).
[12] Fernández-Nieto, E.; Castro, M.; Parés, C., On an intermediate field capturing Riemann solver based on a parabolic viscosity matrix for the two-layer shallow water system, J. Sci. Comput., 48, 1-3, 117-140 (2011) · Zbl 1421.76163
[13] Schijf, J.; Schönfled, J., Theoretical considerations on the motion of salt and fresh water (1953), IAHR
[14] Struchtrup, H., Macroscopic transport equations for rarefied gas flows: approximation methods in kinetic theory, Interaction of Mechanics and Mathematics (2006), Springer Berlin Heidelberg
[15] Torrilhon, M., Modeling nonequilibrium gas flow based on moment equations, Annu. Rev. Fluid Mech., 48, 1, 429-458 (2016) · Zbl 1356.76297
[16] Grad, H., On the kinetic theory of rarefied gases, Commun. Pure Appl. Math., 2, 4, 331-407 (1949) · Zbl 0037.13104
[17] Cai, Z.; Fan, Y.; Li, R., Globally hyperbolic regularization of Grad’s moment system, Commun. Pure Appl. Math., 67, 3, 464-518 (2014) · Zbl 1307.35182
[18] Fan, Y.; Koellermeier, J.; Li, J.; Li, R.; Torrilhon, M., Model reduction of kinetic equations by operator projection, J. Stat. Phys., 162, 2, 457-486 (2016) · Zbl 1335.35175
[19] Koellermeier, J.; Schaerer, R. P.; Torrilhon, M., A framework for hyperbolic approximation of kinetic equations using quadrature-based projection methods, Kinet. Relat. Models, 7, 3, 531-549 (2014) · Zbl 1304.35475
[20] Koellermeier, J.; Torrilhon, M., Numerical solution of hyperbolic moment models for the Boltzmann equation, Eur. J. Mech. B. Fluids, 64, 41-46 (2017) · Zbl 1408.76454
[21] Fan, Y.; Li, R., Globally hyperbolic moment system by generalized Hermite expansion, Sci. Sin. Math., 45(10), 10, 1635-1676 (2015)
[22] Koellermeier, J.; Torrilhon, M., Numerical study of partially conservative moment equations in kinetic theory, Commun. Comput. Phys., 21(04), 4, 981-1011 (2017) · Zbl 1373.82079
[23] Bhatnagar, P. L.; Gross, E. P.; Krook, M., A model for collision processes in gases. 1. Small amplitude processes in charged and neutral one-component systems, Phys. Rev., 94, 511-525 (1954) · Zbl 0055.23609
[24] Cai, Z.; Fan, Y.; Li, R., Globally hyperbolic regularization of Grad’s moment system in one dimensional space, Commun. Math. Sci., 11(2), 547-571 (2013) · Zbl 1301.35083
[25] Koellermeier, J.; Torrilhon, M., Simplified hyperbolic moment equations, Proceedings of the 16th International Conference on Hyperbolic Problems (2016) · Zbl 1408.76454
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.