Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPuncture black hole initial data: a single domain Galerkin-Collocation method for trumpet and wormhole data sets
119
0
0.0
(
0
)
Added by
Henrique de Oliveira
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present a single domain Galerkin-Collocation method to calculate puncture initial data sets for single and binary, either in the trumpet or wormhole geometries. The combination of aspects belonging to the Galerkin and the Collocation methods together with the adoption of spherical coordinates in all cases show to be very effective. We have proposed a unified expression for the conformal factor to describe trumpet and spinning black holes. In particular, for the spinning trumpet black holes, we have exhibited the deformation of the limit surface due to the spin from a sphere to an oblate spheroid. We have also revisited the energy content in the trumpet and wormhole puncture data sets. The algorithm can be extended to describe binary black holes.
rate research
Read More
We present a new computational framework for the Galerkin-collocation method for double domain in the context of ADM 3+1 approach in numerical relativity. This work enables us to perform high resolution calculations for initial sets of two arbitrary black holes. We use the Bowen-York method for binary systems and the puncture method to solve the Hamiltonian constraint. The nonlinear numerical code solves the set of equations for the spectral modes using the standard Newton-Raphson method, LU decomposition and Gaussian quadratures. We show convergence of our code for the conformal factor and the ADM mass. Thus, we display features of the conformal factor for different masses, spins and linear momenta.
We present improvements to construction of binary black hole initial data used in SpEC (the Spectral Einstein Code). We introduce new boundary conditions for the extended conformal thin sandwich elliptic equations that enforce the excision surfaces to be slightly inside rather than on the apparent horizons, thus avoiding extrapolation into the black holes at the last stage of initial data construction. We find that this improves initial data constraint violations near and inside the apparent horizons by about 3 orders of magnitude. We construct several initial data sets that are intended to be astrophysically equivalent but use different free data, boundary conditions, and initial gauge conditions. These include free data chosen as a superposition of two black holes in time-independent horizon-penetrating harmonic and damped harmonic coordinates. We also implement initial data for which the initial gauge satisfies the harmonic and damped harmonic gauge conditions; this can be done independently of the free data, since this amounts to a choice of the time derivatives of the lapse and shift. We compare these initial data sets by evolving them. We show that the gravitational waveforms extracted during the evolution of these different initial data sets agree very well after excluding initial transients. However, we do find small differences between these waveforms, which we attribute to small differences in initial orbital eccentricity, and in initial BH masses and spins, resulting from the different choices of free data. Among the cases considered, we find that superposed harmonic initial data leads to significantly smaller transients, smaller variation in BH spins and masses during these transients, smaller constraint violations, and more computationally efficient evolutions. Finally, we study the impact of initial data choices on the construction of zero-eccentricity initial data.
We present an implementation of the Galerkin-Collocation method to determine the initial data for non-rotating distorted three dimensional black holes in the inversion and puncture schemes. The numerical method combines the key features of the Galerkin and Collocation methods which produces accurate initial data. We evaluated the ADM mass of the initial data sets, and we have provided the angular structure of the gravitational wave distribution at the initial hypersurface by evaluating the scalar $Psi_4$ for asymptotic observers.
We present a Galerkin-Collocation domain decomposition algorithm applied to the evolution of cylindrical unpolarized gravitational waves. We show the effectiveness of the algorithm in reproducing initial data with high localized gradients and in providing highly accurate dynamics. We characterize the gravitational radiation with the standard Newman-Penrose Weyl scalar $Psi_4$. We generate wave templates for both polarization modes, $times$ and $+$, outgoing and ingoing, to show how they exchange energy nonlinearly. In particular, considering an initially ingoing $times$ wave, we were able to trace a possible imprint of the gravitational analog of the Faraday effect in the generated templates.
We propose and explore a stationary 1+log slicing condition for the construction of solutions to Einsteins constraint equations. For stationary spacetimes, these initial data will give a stationary foliation when evolved with moving puncture gauge conditions that are often used in black hole evolutions. The resulting slicing is time-independent and agrees with the slicing generated by being dragged along a time-like Killing vector of the spacetime. When these initial data are evolved with moving puncture gauge conditions, numerical errors arising from coordinate evolution are minimized. In the construction of initial data for binary black holes it is often assumed that there exists an approximate helical Killing vector that generates the binarys orbit. We show that, unfortunately, 1+log slices that are stationary with respect to such a helical Killing vector cannot be asymptotically flat, unless the spacetime possesses an additional axial Killing vector.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
