# 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.
##### MSC:
 65-XX Numerical analysis 15-XX Linear and multilinear algebra; matrix theory
##### Software:
Eigen; PVM ; HLLE
Full Text:
##### References:
  Roe, P., Approximate Riemann solvers, parameter vectors, and difference schemes, J. Comput. Phys., 43, 2, 357-372 (1981) · Zbl 0474.65066  Toumi, I., A weak formulation of Roe’s approximate Riemann solver, J. Comput. Phys., 102, 2, 360-373 (1992) · Zbl 0783.65068  Parés, C., Numerical methods for nonconservative hyperbolic systems: a theoretical framework., SIAM J. Numer. Anal., 44, 1, 300-321 (2006) · Zbl 1130.65089  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  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  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  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  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  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)  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  G. Guennebaud, B. Jacob, et al., Eigen v3, 2010, (http://eigen.tuxfamily.org).  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  Schijf, J.; Schönfled, J., Theoretical considerations on the motion of salt and fresh water (1953), IAHR  Struchtrup, H., Macroscopic transport equations for rarefied gas flows: approximation methods in kinetic theory, Interaction of Mechanics and Mathematics (2006), Springer Berlin Heidelberg  Torrilhon, M., Modeling nonequilibrium gas flow based on moment equations, Annu. Rev. Fluid Mech., 48, 1, 429-458 (2016) · Zbl 1356.76297  Grad, H., On the kinetic theory of rarefied gases, Commun. Pure Appl. Math., 2, 4, 331-407 (1949) · Zbl 0037.13104  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  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  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  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  Fan, Y.; Li, R., Globally hyperbolic moment system by generalized Hermite expansion, Sci. Sin. Math., 45(10), 10, 1635-1676 (2015)  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  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  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  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.