اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدNumerical relativity with characteristic evolution, using six angular patches
56
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular coordinate system has proved effective. Here we investigate the use of a six-patch system in characteristic numerical relativity, by comparing an existing two-patch implementation (using second-order finite differencing throughout) with a new six-patch implementation (using either second- or fourth-order finite differencing for the angular derivatives). We compare these different codes by monitoring the Einstein constraint equations, numerically evaluated independently from the evolution. We find that, compared to the (second-order) two-patch code at equivalent resolutions, the errors of the second-order six-patch code are smaller by a factor of about 2, and the errors of the fourth-order six-patch code are smaller by a factor of nearly 50.
قيم البحث
اقرأ أيضاً
We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at $mathcal{I}^+$ from numerical relativity simulations. Cauchy-characteristic evolution combines an interior soluti
on of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to $mathcal{I}^+$. The foundation of our improved algorithm is a comprehensive method of handling the gauge transformations between the arbitrarily specified coordinates of the interior Cauchy evolution and the unique (up to BMS transformations) Bondi-Sachs coordinate system of the exterior characteristic evolution. We present a reformulated set of characteristic evolution equations better adapted to numerical implementation. In addition, we develop a method to ensure that the angular coordinates used in the volume during the characteristic evolution are asymptotically inertial. This provides a direct route to an expanded set of waveform outputs and is guaranteed to avoid pure-gauge logarithmic dependence that has caused trouble for previous spectral implementations of the characteristic evolution equations. We construct a set of Weyl scalars compatible with the Bondi-like coordinate systems used in characteristic evolution, and determine simple, easily implemented forms for the asymptotic Weyl scalars in our suggested set of coordinates.
We describe a numerical code that solves Einsteins equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used to evolve a
numerical spacetime containing a black hole. We excise the hole from the computational grid in order to avoid the central singularity. We describe in detail a causal differencing method that should allow one to stably evolve a hyperbolic system of equations in three spatial dimensions with an arbitrary shift vector, to second-order accuracy in both space and time. We demonstrate the success of this method in the spherically symmetric case.
The accurate modeling of gravitational radiation is a key issue for gravitational wave astronomy. As simulation codes reach higher accuracy, systematic errors inherent in current numerical relativity wave-extraction methods become evident, and may le
ad to a wrong astrophysical interpretation of the data. In this paper, we give a detailed description of the Cauchy-characteristic extraction technique applied to binary black hole inspiral and merger evolutions to obtain gravitational waveforms that are defined unambiguously, that is, at future null infinity. By this method we remove finite-radius approximations and the need to extrapolate data from the near zone. Further, we demonstrate that the method is free of gauge effects and thus is affected only by numerical error. Various consistency checks reveal that energy and angular momentum are conserved to high precision and agree very well with extrapolated data. In addition, we revisit the computation of the gravitational recoil and find that finite radius extrapolation very well approximates the result at $scri$. However, the (non-convergent) systematic differences to extrapolated data are of the same order of magnitude as the (convergent) discretisation error of the Cauchy evolution hence highlighting the need for correct wave-extraction.
In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial many-boxes-in-many-boxes mesh hierarchi
es and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation using the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. The AMR capability permits the study of a range of new physics which has previously been computationally infeasible in a full 3+1 setting, whilst also significantly simplifying the process of setting up the mesh for these problems. We show that GRChombo can stably and accurately evolve standard spacetimes such as binary black hole mergers and scalar collapses into black holes, demonstrate the performance characteristics of our code, and discuss various physics problems which stand to benefit from the AMR technique.
We present a new expression for the Weyl scalar Psi_4 that can be used in numerical relativity to extract the gravitational wave content of a spacetime. The formula relies upon the identification of transverse tetrads, namely the ones in which Psi_1=
Psi_3=0. It is well known that tetrads with this property always exist in a general Petrov type I spacetime. A sub-class of these tetrads naturally converges to the Kinnersley tetrad in the limit of Petrov type D spacetime. However, the transverse condition fixes only four of the six parameters coming from the Lorentz group of transformations applied to tetrads. Here we fix the tetrad completely, in particular by giving the expression for the spin-boost transformation that was still unclear. The value of Psi_4 in this optimal tetrad is given as a function of the two curvature invariants I and J.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
