1. A well balanced finite volume scheme for general relativity
- Author
-
Gaburro, Elena, Castro, Manuel, Dumbser, Michael, Certified Adaptive discRete moDels for robust simulAtions of CoMplex flOws with Moving fronts (CARDAMOM), Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux Segalen - Bordeaux 2-Université Sciences et Technologies - Bordeaux 1-Université de Bordeaux (UB)-Institut Polytechnique de Bordeaux (Bordeaux INP)-Centre National de la Recherche Scientifique (CNRS)-Université Bordeaux Segalen - Bordeaux 2-Université Sciences et Technologies - Bordeaux 1-Université de Bordeaux (UB)-Institut Polytechnique de Bordeaux (Bordeaux INP)-Centre National de la Recherche Scientifique (CNRS)-Inria Bordeaux - Sud-Ouest, Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria), Universidad de Málaga [Málaga] = University of Málaga [Málaga], University of Trento [Trento], E. G gratefully acknowledges the support received from the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Curie Individual Fellowship SuPerMan, grant agreement No. 101025563. E. G. has been also supported by a national mobility grant for young researchers in Italy funded by GNCS-INdAM, and she received funding from the University of Trento via the Strategic Initiative Starting Grant Giovani Ricercatori 2019. Moreover, the research presented in this paper has been partially funded by the European Union’s Horizon 2020 Research and Innovation Programme under the project ExaHyPE, grant No. 671698 (call FETHPC-1-2014).M. D. also acknowledges the financial support received from the Italian Ministry of Education, University and Research (MIUR) in the frame of the Departments of Excellence Initiative 2018–2022 attributed to DICAM of the University of Trento (grant 605 L. 232/2016) and in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications. Furthermore, M. D. has also received funding from the University of Trento via the Strategic Initiative Modeling and Simulation. The research of M.J. Castro was partially supported by the Spanish Government (SG), the European Regional Development Fund (ERDF), the Regional Government of Andalusia (RGA), and the University of Malaga(UMA) through the projects of references RTI2018-096064-B-C21 (SG-ERDF), UMA18-Federja-161 (RGA-ERDF-613 UMA), and P18-RT-3163 (RGA-ERDF)., and Université Bordeaux Segalen - Bordeaux 2-Université Sciences et Technologies - Bordeaux 1 (UB)-Université de Bordeaux (UB)-Institut Polytechnique de Bordeaux (Bordeaux INP)-Centre National de la Recherche Scientifique (CNRS)-Université Bordeaux Segalen - Bordeaux 2-Université Sciences et Technologies - Bordeaux 1 (UB)-Université de Bordeaux (UB)-Institut Polytechnique de Bordeaux (Bordeaux INP)-Centre National de la Recherche Scientifique (CNRS)-Inria Bordeaux - Sud-Ouest
- Subjects
010308 nuclear & particles physics ,Applied Mathematics ,FOS: Physical sciences ,Well balanced schemes 23 (WB) ,General Relativity and Quantum Cosmology (gr-qc) ,Numerical Analysis (math.NA) ,01 natural sciences ,General Relativity and Quantum Cosmology ,General relativistic magnetohydrodynamics (GRMHD) ,TOV neutron star ,AMS subject classifications. 35L40, 65M08, 83C05, 83C10, 85-08, 85-10, 85A30 ,Computational Mathematics ,0103 physical sciences ,FOS: Mathematics ,Michel accretion disk ,Mathematics - Numerical Analysis ,[MATH]Mathematics [math] ,010306 general physics ,Frst order conformal and covariant 24 reformulation of the Einstein field equations (FO-CCZ4) ,Fnite volume schemes (FV) ,First order hyperbolic systems - Abstract
International audience; In this work we present a novel second order accurate well balanced (WB) finite volume (FV) scheme for the solution of the general relativistic magnetohydrodynamics (GRMHD) equations and the first order CCZ4 formulation (FO-CCZ4) of the Einstein field equations of general relativity, as well as the fully coupled FO-CCZ4 + GRMHD system. These systems of first order hyperbolic PDEs allow to study the dynamics of the matter and the dynamics of the space-time according to the theory of general relativity. The new well balanced finite volume scheme presented here exploits the knowledge of an equilibrium solution of interest when integrating the conservative fluxes, the nonconservative products and the algebraic source terms, and also when performing the piecewise linear data reconstruction. This results in a rather simple modification of the underlying second order FV scheme, which, however, being able to cancel numerical errors committed with respect to the equilibrium component of the numerical solution, substantially improves the accuracy and long-time stability of the numerical scheme when simulating small perturbations of stationary equilibria. In particular, the need for well balanced techniques appears to be more and more crucial as the applications increase their complexity. We close the paper with a series of numerical tests of increasing difficulty, where we study the evolution of small perturbations of accretion problems and stable TOV neutron stars. Our results show that the well balancing significantly improves the long-time stability of the finite volume scheme compared to a standard one.
- Published
- 2021
- Full Text
- View/download PDF