No Arabic abstract
We present a symmetric hyperbolic formulation of the Einstein equations in affine-null coordinates. Giannakopoulos et. al. (arXiv:2007.06419) recently showed that the most commonly numerically implemented formulations of the Einstein equations in affine-null coordinates (and other single-null coordinate systems) are only weakly-but not strongly-hyperbolic. By making use of the tetrad-based Newman-Penrose formalism, our formulation avoids the hyperbolicity problems of the formulations investigated by Giannakopoulos et. al. We discuss a potential application of our formulation for studying gravitational wave scattering.
A new evolution algorithm for the characteristic initial value problem based upon an affine parameter rather than the areal radial coordinate used in the Bondi-Sachs formulation is applied in the spherically symmetric case to the gravitational collapse of a massless scalar field. The advantages over the Bondi-Sachs version are discussed, with particular emphasis on the application to critical collapse. Unexpected quadratures lead to a simple evolution algorithm based upon ordinary differential equations which can be integrated along the null rays. For collapse to a black hole in a Penrose compactified spacetime, these equations are regularized throughout the exterior and interior of the horizon up to the final singularity. They are implemented as a global numerical evolution code based upon the Galerkin method. New results regarding the global properties of critical collapse are presented.
We show that the spherically symmetric Einstein-scalar-field equations for wave-like decaying initial data at null infinity have unique global solutions in (0, infty) and unique generalized solutions on [0, infty) in the sense of Christodoulou. We emphasize that this decaying condition is sharp.
In this work we give a complete picture of how to in a direct simple way define the mass at null infinity in harmonic coordinates in three different ways that we show satisfy the Bondi mass loss law. The first and second way involve only the limit of metric (Trautman mass) respectively the null second fundamental forms along asymptotically characteristic surfaces (asymptotic Hawking mass) that only depend on the ADM mass. The last in an original way involves construction of special characteristic coordinates at null infinity (Bondi mass). The results here rely on asymptotics of the metric derived in [24].
A number of astrophysical scenarios possess and preserve an overall cylindrical symmetry also when undergoing a catastrophic and nonlinear evolution. Exploiting such a symmetry, these processes can be studied through numerical-relativity simulations at smaller computational costs and at considerably larger spatial resolutions. We here present a new flux-conservative formulation of the relativistic hydrodynamics equations in cylindrical coordinates. By rearranging those terms in the equations which are the sources of the largest numerical errors, the new formulation yields a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. We illustrate this through a series of numerical tests involving the evolution of oscillating spherical and rotating stars, as well as shock-tube tests.
Bonazzola, Gourgoulhon, Grandclement, and Novak [Phys. Rev. D {bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for the deviation of a conformal metric from a flat fiducial one. The choice of a Diracs gauge for the spatial coordinates guarantees the mathematical characterization of that system as a (strongly) hyperbolic system of conservation laws. In the presence of boundaries, this characterization also depends on the boundary conditions for the shift vector in the elliptic subsystem. This interplay between the hyperbolic and elliptic parts of the complete evolution system is used to assess the prescription of inner boundary conditions for the hyperbolic part when using an excision approach to black hole spacetime evolutions.