An implementation of the microphysics in full general relativity: a general relativistic neutrino leakage scheme
Abstract
Performing fully general relativistic simulations taking account of microphysical processes (e.g. weak interactions and neutrino cooling) is one of the long-standing problems in numerical relativity. One of main difficulties in implementation of weak interactions in the general relativistic framework lies in the fact that the characteristic timescale of weak interaction processes (the WP timescale, t_wp \sim \vert Y_{e}/\dot{Y}_{e} \vert) in hot dense matters is much shorter than the dynamical timescale (tdyn). Numerically this means that stiff source terms appear in the equations so that an implicit scheme is in general necessary to stably solve the relevant equations. Otherwise a very short timestep (Δt < twp Lt tdyn) will be required to solve them explicitly, which is unrealistic in the present computational resources. Furthermore, in the relativistic framework, the Lorentz factor is coupled with the rest mass density and the energy density. The specific enthalpy is also coupled with the momentum. Due to these couplings, it is very complicated to recover the primitive variables and the Lorentz factor from conserved quantities. Consequently, it is very difficult to solve the equations implicitly in the fully general relativistic framework. At the current status, no implicit procedure has been proposed except for the case of the spherical symmetry. Therefore, an approximate explicit procedure is developed in the fully general relativistic framework in this paper as a first implementation of the microphysics toward a more realistic sophisticated model. The procedure is based on the so-called neutrino leakage schemes which are based on the property that the characteristic timescale in which neutrinos leak out of the system (the leakage timescale, tleak) is much longer than the WP timescale. In the previous leakage schemes, however, the problems of the stiff source terms are avoided in an artificial manner. In this paper, I present a detailed neutrino leakage scheme and a simple and stable method for solving the equations explicitly in the fully general relativistic framework. The drawback of the artificial treatment of the stiff source terms is improved. I also perform a test simulation to check the validity of the present method, showing that it works fairly well.
- Publication:
-
Classical and Quantum Gravity
- Pub Date:
- June 2010
- DOI:
- 10.1088/0264-9381/27/11/114107
- arXiv:
- arXiv:1009.3358
- Bibcode:
- 2010CQGra..27k4107S
- Keywords:
-
- Astrophysics - High Energy Astrophysical Phenomena;
- General Relativity and Quantum Cosmology
- E-Print:
- An invited paper from Numerical Relativity and Data Analysis (NRDA) 2009, Albert Einstein Institute, Potsdam Finer-resolution-figure version will be found in the journal reference