Using piecewise parabolic reconstruction of physical variables in the Rusanov solver. I. The special relativistic hydrodynamics equations

I. M. Kulikov

Sibirskij žurnal industrialʹnoj matematiki, Tome 26 (2023) no. 4, pp. 49-64.

Résumé

The Rusanov solver for solving hydrodynamic equations is one of the most robust schemes in the class of Riemann solvers. For special relativistic hydrodynamics, the robustness condition of the scheme is the most important property, especially for sufficiently high values of the Lorentz factor. At the same time, the Rusanov solver is known to be very dissipative. It is proposed to use a piecewise parabolic representation of physical variables to reduce the dissipation of the Rusanov scheme. Using this approach has made it possible to obtain a scheme with the same dissipative properties as Roe-type schemes and the family of Harten—Lax—van Leer schemes. Using the problem of the decay of a relativistic hydrodynamic discontinuity, it is shown that the present author's version of the Rusanov scheme is advantageous in terms of reproducing a contact discontinuity. The scheme is verified on classical problems of discontinuity decay and on the problem of the interaction of two relativistic jets in the three-dimensional formulation.

Détail

Keywords: computational astrophysics, relativistic hydrodynamics, numerical method, Rusanov solver, piecewise parabolic method.

Bibliographie

[1] Chechetkin V. M., D'yachenko V. F., Ginzburg S. L., Paleichik V. V., Fimin N. N., Sudarikov A. L., “On the generation mechanism of hard cosmic gamma-ray emission from AGN jets”, Astron. Rep, 53 (2009), 501–509 | DOI

[2] Barkov M. V., Bisnovatyi-Kogan G. S., “Interaction of a cosmological gamma-ray burst with a dense molecular cloud and the formation of jets”, Astron. Rep, 49 (2005), 24–35 | DOI

[3] Sotomayor P., Romero G., “Nonthermal radiation from the central region of super-accreting active galactic nuclei”, Astron. Astrophys, 664 (2022), A178 | DOI

[4] Sokolov V. V., Bisnovatyi-Kogan G. S., Kurt V. G., Gnedin Yu. N., Baryshev Yu. V., “Observational constraints on the angular and spectral distributions of photons in gamma-ray burst sources”, Astron. Rep, 50 (2006), 612–625 | DOI

[5] Komissarov S. S., “Simulations of the axisymmetric magnetospheres of neutron stars”, Mon. Not. R. Astron. Soc, 367 (2006), 19–31 | DOI

[6] Tutukov A. V., Bogomazov A. I., “Radio pulsars in close binaries with neutron stars”, Astron. Rep, 52 (2008), 390–402 | DOI

[7] Istomin Ya. N., Komberg B. V., “Gamma-ray bursts as a result of the interaction of a shock from a supernova and a neutron-star companion”, Astron. Rep, 46 (2002), 908–917 | DOI

[8] Fateev V. F., Davlatov R. A., “Space-Based Gravitational-Wave Detectors: Development of Ground-Breaking Technologies for Future Space-Based Gravitational Gradiometers”, Astron. Rep, 63 (2019), 699–709 | DOI

[9] Belyaev V. S., Bisnovatyi-Kogan G. S., Gromov A. I., Zagreev B. V., Lobanov A. V., Matafonov A. P., Moiseenko S. G., Toropina O. D., “Numerical Simulations of Magnetized Astrophysical Jets and Comparison with Laboratory Laser Experiments”, Astron. Rep, 62 (2018), 162–182 | DOI

[10] Krauz V. I., Mitrofanov K. N., Kharrasov A. M., Il'ichev I. V., Myalton V. V., Anan'ev S. S., Beskin V. S., “Laboratory Modeling of the Rotation of Jets Ejected from Young Stellar Objects at Studies the Azimuthal Structure of an Axial Jet at the PF-3 Facility”, Astron. Rep, 65 (2021), 26–44 | DOI

[11] Kulikov I., “A new code for the numerical simulation of relativistic flows on supercomputers by means of a low-dissipation scheme”, Comput. Phys. Commun, 257 (2020), 107532 | DOI | MR | Zbl

[12] Rusanov V. V., “The calculation of the interaction of non-stationary shock waves with barriers”, Comput. Math. Math. Phys, 1 (1961), 267–279 | MR

[13] Mohamed K., Benkhaldoun F., “A modified Rusanov scheme for shallow water equations with topography and two phase flows”, Eur. Phys. J. Plus, 131 (2016), 207 | DOI

[14] Wu C., Walker R., Dawson J., “A three dimensional MHD model of the Earth's magnetosphere”, Geophys. Res. Lett., 8:5 (1981), 523–526 | DOI

[15] Mohamed K., Abdelrahman M., “The modified Rusanov scheme for solving the ultra-relativistic Euler equations”, Eur. J. Mech. B Fluids, 90 (2021), 89–98 | DOI | MR | Zbl

[16] Mohammadian S., Moghaddam A., Sahaf A., “On the performance of HLL, HLLC, and Rusanov solvers for hyperbolic traffic models”, Comput. Fluids, 231 (2021), 105161 | DOI | MR | Zbl

[17] Nishikawa H., Kitamura K., “Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers”, J. Comput. Phys, 227 (2008), 2560–2581 | DOI | MR | Zbl

[18] Townsend J., Koenoezsy L., Jenkins K., “On the development of a rotated-hybrid HLL/HLLC approximate Riemann solver for relativistic hydrodynamics”, Mon. Not. R. Astron. Soc., 496:2 (2020), 2493–2505 | DOI

[19] Jaisankar S., Raghurama Rao S., “Diffusion regulation for Euler solvers”, J. Comput. Phys., 221:2 (2007), 577–599 | DOI | MR | Zbl

[20] Mazaheri A., Nishikawa H., “High-order shock-capturing hyperbolic residual-distribution schemes on irregular triangular grids”, Comput. Fluids, 131 (2016), 29–44 | DOI | MR | Zbl

[21] Rossi G., Dumbser M., Armanini A., “A well-balanced path conservative SPH scheme for nonconservative hyperbolic systems with applications to shallow water and multi-phase flows”, Comput. Fluids, 154 (2017), 102–122 | DOI | MR | Zbl

[22] Balsara D., “Three dimensional HLL Riemann solver for conservation laws on structured meshes; Application to Euler and magnetohydrodynamic flows”, J. Comput. Phys, 295 (2015), 1–23 | DOI | MR | Zbl

[23] Popov M., Ustyugov S., “Piecewise parabolic method on local stencil for gasdynamic simulations”, Comput. Math. Math. Phys., 47:12 (2007), 1970–1989 | DOI | MR

[24] Popov M., Ustyugov S., “Piecewise parabolic method on a local stencil for ideal magnetohydrodynamics”, Comput. Math. Math. Phys., 48:3 (2008), 477–499 | DOI | MR | Zbl

[25] Kulikov I., Vorobyov E., “Using the PPML approach for constructing a low-dissipation, operator-splitting scheme for numerical simulations of hydrodynamic flows”, J. Comput. Phys., 317 (2016), 318–346 | DOI | MR | Zbl

[26] Kulikov I., Chernykh I., Tutukov A., “A New Hydrodynamic Code with Explicit Vectorization Instructions Optimizations that Is Dedicated to the Numerical Simulation of Astrophysical Gas Flow. I. Numerical Method, Tests, and Model Problems”, Astrophys. J. Suppl. Ser, 243 (2019), 4 | DOI

[27] Kriksin Y. A., Tishkin V. F., “Variational Entropic Regularization of the Discontinuous Galerkin Method for Gasdynamic Equations”, Math. Models Comput. Simul., 11 (2019), 1032–1040 | DOI | MR

[28] Kulikov I., Chernykh I., Karavaev D., Prigarin V., Sapetina A., Ulyanichev I., Zavyalov O., “A New Parallel Code Based on a Simple Piecewise Parabolic Method for Numerical Modeling of Colliding Flows in Relativistic Hydrodynamics”, Mathematics, 10:11 (2022), 1865 | DOI

[29] Kulikov I. M., Karavaev D. A., “A Piecewise-Parabolic Reconstruction of the Physical Variables in a Low-Dissipation HLL Method for the Numerical Solution of the Equations of Special Relativistic Hydrodynamics”, Numer. Anal. Appl, 16 (2023), 45–60 | DOI | MR

[30] Kulikov I. M., “A Low-Dissipation Numerical Scheme Based on a Piecewise Parabolic Method on a Local Stencil for Mathematical Modeling of Relativistic Hydrodynamic Flows”, Numer. Anal. Appl, 13 (2020), 117–126 | DOI | MR

[31] Lee D., Faller H., Reyes A., “The Piecewise Cubic Method (PCM) for computational fluid dynamics”, J. Comput. Phys, 341 (2017), 230–257 | DOI | MR | Zbl

[32] Wang P., Abel T., Zhang W., “Relativistic Hydrodynamic Flows Using Spatial and Temporal Adaptive Structured Mesh Refinement”, Astrophys. J. Suppl. Ser, 176 (2008), 467–483 | DOI

[33] Bhoriya D., Kumar H., “Entropy-stable schemes for relativistic hydrodynamics equations”, Z. Angew. Math. Phys, 71 (2020), 29 | DOI | MR | Zbl

[34] Guercilena F., Radice D., Rezzolla L., “Entropy-limited hydrodynamics: a novel approach to relativistic hydrodynamics”, Comput. Astrophys. Cosmol, 4 (2017), 3 | DOI

[35] Zlotnik A. A., “Spatial discretization of the one-dimensional quasi-gasdynamic system of equations and the entropy balance equation”, Comput. Math. Math. Phys, 52 (2012), 1060–1071 | DOI | MR | Zbl

[36] Gavrilin V. A., Zlotnik A. A., “On spatial discretization of the one-dimensional quasi-gasdynamic system of equations with general equations of state and entropy balance”, Comput. Math. Math. Phys, 55 (2015), 264–281 | DOI | MR | Zbl

[37] Zlotnik A. A., “Entropy-conservative spatial discretization of the multidimensional quasi-gasdynamic system of equations”, Comput. Math. Math. Phys, 57 (2017), 706–725 | DOI | MR | Zbl