Exploring the relativistic regime with Newtonian hydrodynamics: An improved effective gravitational potential for supernova simulations

A. Marek,H. Dimmelmeier,H.-Th. Janka,E. Mueller,R. Buras
DOI: https://doi.org/10.1051/0004-6361%3A20052840
2005-07-13
Abstract:We investigate the possibility to approximate relativistic effects in hydrodynamical simulations of stellar core collapse and post-bounce evolution by using a modified gravitational potential in an otherwise standard Newtonian hydrodynamic code. Different modifications of the effective relativistic potential introduced by Rampp & Janka (2002) are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hydrodynamic simulations with two different codes are compared with fully relativistic results. One code is applied for one- and two-dimensional calculations with a simple equation of state, and employs either the modified effective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale supernova runs including a microphysical equation of state and neutrino transport based on the solution of the Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully relativistic solutions in spherical symmetry, leading to significant improvements compared to previously published approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation.
Astrophysics,General Relativity and Quantum Cosmology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to approximate relativistic effects with an improved effective gravitational potential in the simulation of stellar core collapse and supernova explosions. Specifically: 1. **Importance of Relativistic Effects**: During the process of supernova explosions, especially when the proto - neutron star is formed, the density of matter is extremely high, causing the effects of general relativity to become very important and not negligible. Therefore, it is crucial to correctly handle or approximate these relativistic effects in numerical simulations. 2. **Existing Challenges**: Previous studies (such as the study by Liebendörfer et al. in 2005) have shown that using Newtonian hydrodynamics codes combined with modified gravitational potentials can partially simulate relativistic effects, but in some cases, especially in the post - bounce phase, the results have quantitative differences from fully relativistic simulations. 3. **Objectives of Improvement**: In order to improve the simulation accuracy within the Newtonian framework, the authors have tried several different modification schemes of the effective gravitational potential and evaluated their performance under different conditions to find a method that can better approximate relativistic effects. ### Specific Problem Description The main objective of the paper is to explore how to more accurately simulate relativistic effects by modifying the gravitational potential in Newtonian hydrodynamics codes. Specifically, the authors hope to find a method that can more precisely reproduce the results of fully relativistic simulations while maintaining the simplicity of Newtonian dynamics. To this end, they have carried out the following tasks: - **Developing an Improved Effective Gravitational Potential**: Based on the Tolman - Oppenheimer - Volkoff (TOV) equation, several different forms of effective gravitational potential are proposed, and detailed mathematical derivations are carried out for each form. \[ \Phi_{\text{TOV}}(r)=- 4\pi\int_{r}^{\infty}dr'\frac{r'^{2}(m_{\text{TOV}}+r'^{3}(P + p_{\nu}))}{\Gamma^{2}(\rho + e+P)} \] - **Application in Multidimensional Flows**: For multidimensional flows, the authors replace the "spherically symmetric contribution" with the TOV potential energy to adapt to non - spherically symmetric situations. - **Testing of Multiple Modification Schemes**: In order to reduce the problem of overly strong TOV potential energy, the authors have tested multiple modification schemes (such as Case A to Case G) and compared their effects. - **Numerical Simulation Verification**: A large number of numerical experiments have been carried out using two different simulation codes (CoCoNuT and Vertex) to verify the performance of different effective gravitational potentials and compare them with the results of fully relativistic simulations. ### Conclusion Through these improvements, the authors have found that some modification schemes (such as Case A) can greatly improve the simulation results, making them closer to fully relativistic simulations. This not only improves the accuracy of static solutions but also performs well during the dynamic evolution process. Therefore, this study provides a more reliable and effective tool for future supernova simulations.