BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeffrey Winicour

Jeffrey Winicour appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.COhep-th

Catalog footprint

What is connected

16works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2016arXiv

On solving the constraints by integrating a strongly hyperbolic system

It was shown recently that the constraints on the initial data for Einstein's equations may be posed as an evolutionary problem [9]. In one of the proposed two methods the constraints can be replaced by a first order symmetrizable hyperbolic system and a subsidiary algebraic relation. Here, by assuming that the initial data surface is smoothly foliated by a one-parameter family of topological two-spheres, the basic variables are recast in terms of spin-weighted fields. This allows one to replace all the angular derivatives in the evolutionary system by the Newman-Penrose $ð$ and $\barð$ operators which, in turn, opens up a new avenue to solve the constraints by integrating the resulting system using suitable numerical schemes. In particular, by replacing the $ð$ and $\barð$ operators either by a finite difference or by a pseudo-spectral representation or by applying a spectral decomposition in terms of spin-weighted spherical harmonics, the evolutionary equations may be put into the form of a coupled system of non-linear ordinary differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spectral Cauchy Characteristic Extraction of strain, news and gravitational radiation flux

We present a new approach for the Cauchy-characteristic extraction of gravitational radiation strain, news function, and the flux of the energy-momentum, supermomentum and angular momentum associated with the Bondi-Metzner-Sachs asymptotic symmetries. In Cauchy-characteristic extraction, a characteristic evolution code takes numerical data on an inner worldtube supplied by a Cauchy evolution code, and propagates it outwards to obtain the space-time metric in a neighborhood of null infinity. The metric is first determined in a scrambled form in terms of coordinates determined by the Cauchy formalism. In prior treatments, the waveform is first extracted from this metric and then transformed into an asymptotic inertial coordinate system. This procedure provides the physically proper description of the waveform and the radiated energy but it does not generalize to determine the flux of angular momentum or supermomentum. Here we formulate and implement a new approach which transforms the full metric into an asymptotic inertial frame and provides a uniform treatment of all the radiation fluxes associated with the asymptotic symmetries. Computations are performed and calibrated using the Spectral Einstein Code (SpEC).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The sky pattern of the linearized gravitational memory effect

The gravitational memory effect leads to a net displacement in the relative positions of test particles. This memory is related to the change in the strain of the gravitational radiation field between infinite past and infinite future retarded times. There are three known sources of the memory effect: (i) the loss of energy to future null infinity by massless fields or particles, (ii) the ejection of massive particles to infinity from a bound system and (iii) homogeneous, source-free gravitational waves. In the context of linearized theory, we show that asymptotic conditions controlling these known sources of the gravitational memory effect rule out any other possible sources with physically reasonable stress-energy tensors. Except for the source-free gravitational waves, the two other known sources produce gravitational memory with E-mode radiation strain, characterized by a certain curl-free sky pattern of their polarization. Thus our results show that the only known source of B-mode gravitational memory is of primordial origin, corresponding in the linearized theory to a homogeneous wave entering from past null infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole initial data without elliptic equations

We explore whether a new method to solve the constraints of Einstein's equations, which does not involve elliptic equations, can be applied to provide initial data for black holes. We show that this method can be successfully applied to a nonlinear perturbation of a Schwarzschild black hole by establishing the well-posedness of the resulting constraint problem. We discuss its possible generalization to the boosted, spinning multiple black hole problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gauge Invariant Spectral Cauchy Characteristic Extraction

We present gauge invariant spectral Cauchy characteristic extraction. We compare gravitational waveforms extracted from a head-on black hole merger simulated in two different gauges by two different codes. We show rapid convergence, demonstrating both gauge invariance of the extraction algorithm and consistency between the legacy Pitt null code and the much faster Spectral Einstein Code (SpEC).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Boundary Conditions for the Gravitational Field

A review of the treatment of boundaries in general relativity is presented with the emphasis on application to the formulations of Einstein's equations used in numerical relativity. At present, it is known how to treat boundaries in the harmonic formulation of Einstein's equations and a tetrad formulation of the Einstein-Bianchi system. However, a universal approach valid for other formulations is not in hand. In particular, there is no satisfactory boundary theory for the 3+1 formulations which have been highly successful in binary black hole simulation. I discuss the underlying problems that make the initial-boundary value problem much more complicated than the Cauchy problem. I review the progress that has been made and the important open questions that remain.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Worldtube conservation laws for the null-timelike evolution problem

I treat the worldtube constraints which arise in the null-timelike initial-boundary value problem for the Bondi-Sachs formulation of Einstein's equations. Boundary data on a worldtube and initial data on an outgoing null hypersurface determine the exterior spacetime by integration along the outgoing null geodsics. The worldtube constraints are a set of conservation laws which impose conditions on the integration constants. I show how these constraints lead to a well-posed initial value problem governing the extrinsic curvature of the worldtube, whose components are related to the integration constants. Possible applications to gravitational waveform extraction and to the well-posedness of the null-timelike initial-boundary value problem are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Disembodied boundary data for Einstein's equations

A strongly well-posed initial boundary value problem based upon constraint-preserving boundary conditions of the Sommerfeld type has been established for the harmonic formulation of the vacuum Einstein's equations. These Sommerfeld conditions have been previously presented in a 4-dimensional geometric form. Here we recast the associated boundary data as 3-dimensional tensor fields intrinsic to the boundary. This provides a geometric presentation of the boundary data analogous to the 3-dimensional presentation of Cauchy data in terms of 3-metric and extrinsic curvature. In particular, diffeomorphisms of the boundary data lead to vacuum spacetimes with isometric geometries. The proof of well-posedness is valid for the harmonic formulation and its generalizations. The Sommerfeld conditions can be directly applied to existing harmonic codes which have been used in simulating binary black holes, thus ensuring boundary stability of the underlying analytic system. The geometric form of the boundary conditions also allows them to be formally applied to any metric formulation of Einstein's equations, although well-posedness of the boundary problem is no longer ensured. We discuss to what extent such a formal application might be implemented in a constraint preserving manner to 3+1 formulations, such as the Baumgarte-Shapiro-Shibata-Nakamura system which has been highly successful in binary black hole simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black spacetime. A prime application of characteristic evolution is to compute waveforms via Cauchy-characteristic matching, which is also reviewed.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Some mathematical problems in numerical relativity

The main goal of numerical relativity is the long time simulation of highly nonlinear spacetimes that cannot be treated by perturbation theory. This involves analytic, computational and physical issues. At present, the major impasses to achieving global simulations of physical usefulness are of an analytic/computational nature. We present here some examples of how analytic insight can lend useful guidance for the improvement of numerical approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Well-Posed Initial-Boundary Evolution in General Relativity

Maximally dissipative boundary conditions are applied to the initial-boundary value problem for Einstein's equations in harmonic coordinates to show that it is well-posed for homogeneous boundary data and for boundary data that is small in a linearized sense. The method is implemented as a nonlinear evolution code which satisfies convergence tests in the nonlinear regime and is robustly stable in the weak field regime. A linearized version has been stably matched to a characteristic code to compute the gravitational waveform radiated to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Boundary conditions in linearized harmonic gravity

We investigate the initial-boundary value problem for linearized gravitational theory in harmonic coordinates. Rigorous techniques for hyperbolic systems are applied to establish well-posedness for various reductions of the system into a set of six wave equations. The results are used to formulate computational algorithms for Cauchy evolution in a 3-dimensional bounded domain. Numerical codes based upon these algorithms are shown to satisfy tests of robust stability for random constraint violating initial data and random boundary data; and shown to give excellent performance for the evolution of typical physical data. The results are obtained for plane boundaries as well as piecewise cubic spherical boundaries cut out of a Cartesian grid.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Exact Solutions for the Intrinsic Geometry of Black Hole Coalescence

We describe the null geometry of a multiple black hole event horizon in terms of a conformal rescaling of a flat space null hypersurface. For the prolate spheroidal case, we show that the method reproduces the pair-of-pants shaped horizon found in the numerical simulation of the head-on-collision of black holes. For the oblate case, it reproduces the initially toroidal event horizon found in the numerical simulation of collapse of a rotating cluster. The analytic nature of the approach makes further conclusions possible, such as a bearing on the hoop conjecture. From a time reversed point of view, the approach yields a description of the past event horizon of a fissioning white hole, which can be used as null data for the characteristic evolution of the exterior space-time.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The Asymmetric Merger of Black Holes

We study event horizons of non-axisymmetric black holes and show how features found in axisymmetric studies of colliding black holes and of toroidal black holes are non-generic and how new features emerge. Most of the details of black hole formation and black hole merger are known only in the axisymmetric case, in which numerical evolution has successfully produced dynamical space-times. The work that is presented here uses a new approach to construct the geometry of the event horizon, not by locating it in a given spacetime, but by direct construction. In the axisymmetric case, our method produces the familiar pair-of-pants structure found in previous numerical simulations of black hole mergers, as well as event horizons that go through a toroidal epoch as discovered in the collapse of rotating matter. The main purpose of this paper is to show how new - substantially different - features emerge in the non-axisymmetric case. In particular, we show how black holes generically go through a toroidal phase before they become spherical, and how this fits together with the merger of black holes.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

High-powered Gravitational News

We describe the computation of the Bondi news for gravitational radiation. We have implemented a computer code for this problem. We discuss the theory behind it as well as the results of validation tests. Our approach uses the compactified null cone formalism, with the computational domain extending to future null infinity and with a worldtube as inner boundary. We calculate the appropriate full Einstein equations in computational eth form in (a) the interior of the computational domain and (b) on the inner boundary. At future null infinity, we transform the computed data into standard Bondi coordinates and so are able to express the news in terms of its standard $N_{+}$ and $N_{\times}$ polarization components. The resulting code is stable and second-order convergent. It runs successfully even in the highly nonlinear case, and has been tested with the news as high as 400, which represents a gravitational radiation power of about $10^{13}M_{\odot}/sec$.

0 citations0 reviewsNo code linkedOpen access